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PREFACE 

The  scientific-technical  collection  submitted  to  the  readers' 
attention  unites  thematically  a  group  of  articles  dedicated  to  the 
different  aspects  of  an  increase  in  efficiency  of  contemporary 
[DTRD  *  turbofan  engine]  (ATP/O  with  a  high  bypass  ratio. 

They  include  theoretical  studies  of  different  methods  of 
improvement  of  the  thermodynamic  cycle  of  DTRD  and  of  the  optimum 
selection  of  operating  conditions  parameters,  as  a  result  of  which 
specific  fuel  consumption  decreases.  To  these  methods  on  one  hand 
one  ought  to  relate  the  use  of  high  bypass  ratios  on  the  basis  of 
the  forced  operating  conditions  according  to  parameters,  (temperature 
of  the  gas  before  the  turbine,  the  compression  ratio  of  the 
compressor)  of  the  cycle  of  engine;  on  the  other  hand  -  the 
application  of  heat  recovery  in  DTRD  with  high  bypass  ratio  on  the 
basis  of  a  high-temperature  cycle  with  a  low  compression  ratio. 

A  comparative  efficiency  of  these  two  alternative  methods  to 
a  considerable  degree  will  be  determined  by  design  perfection  and 
specific  weight,  and  operational  reliability  of  the  accomplished 
power  plants. 

A  series  of  articles  is  dedicated  to  the  methods  of  the 
experimental  and  theoretical  investigation  of  DTRD  with  the  booster¬ 
mixing  chambers  and  to  the  creation,  for  this  purpose,  of  original 
experimental  installations. 
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Finally,  a  series  of  articles  Is  dedicated  to  the  extremely 
acute  and  actual  problem  of  contemporary  civil  aviation  -  to  the 
problem  of  a  reduction  DTRD  noise  at  large  bypass  ratio  values. 

In  these  articles  the  basic  physical  characteristics  of  the 
noise  of  aviation  [GTE  *  gas  turbine  engine]  (CTA)  and  the  mechanism 
of  noise  formation  at  the  basic  sources  (reactive  Jet,  compressor, 
fan)  are  examined;  the  range  of  bypass  ratios  is  determined,  in 
which  the  noise  of  fan  turns  out  to  be  predominant;  the  methods 
of  calculation  of  the  noise  of  reactive  jets  and  compressor  are 
examined. 

Furthermore,  in  these  articles  are  examined:  the  role  and 
place  of  full-scale  experiment  in  acoustic  studies,  possible 
schemes  of  the  open  stands  created  for  this  purpose,  the  enumeration 
of  studies  which  can  be  executed  on  such  stands. 

The  published  articles  are  a  summarization  of  several 
scientific  Investigations  made  in  the  Department  for  the  Theory 
of  Aircraft  Engines  [RKIIGA]  (PHMMrA)  during  1968-69. 

The  present  scientific-technical  collection  is  of  interest 
to  specialists  who  work  in  the  area  of  aircraft  engine  construction, 
and  also  to  the  graduates  and  students  of  aviation  [VTUZ's]  (BTy3). 


From  editor 
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AN  ENERGY  STUDY  OF  THE  EFFICIENCY  OF  TURBOFAN 
JET  ENGINES  AT  SUBSONIC  SPEEDS  OF  FLIGHT 


A.  L.  Klyachkin,  A.  Ya.  Dantsygs 


In  this  article  the  results  of  the  calcu¬ 
lated  study  of  the  effect  of  operating  condi¬ 
tion  parameters  and  the  bypass  ratio  on  the 
efficiencies  of  useful  action  of  prospective 
non-boosted  DTRD  at  the  subsonic  speeds  of 
flight  are  examined. 

The  physical  Interpretation  of  the 
obtained  laws  is  given,  and  recommendations 
regarding  the  development  of  DTRD  parameters 
which  provide  its  best  cost-effectiveness  are 
alS'..  given. 


It  is  known  that  the  efficiency  of  a  DTRD  as  any  Jet  engine, 
it  is  estimated  at  three  efficiencies:  effective  (ne),  thrust  (nR) 
and  total  (nQ),  where 

=  (1) 

At  a  fixed  flying  speed  total  efficiency  [KPD]  (Hn/O  uniquely 
determines  the  specific  fuel  consumption  that  follows  from  the 
relationship: 

nr-  (2) 

The  advantage  of  non-boosted  DTRD  over  the  turbojet  engine 
according  to  cost-effectiveness  is  explained  by  the  fact  that 
the  introduction  of  the  secondary  circuit  during  the  optimum  or 
close  to  optimum  energy  distribution  between  contour  retards  the 
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outflow  of  gas  from  the  jet  nozzles  of  the  contours.  Thus,  the 
thrust  engine  efficiency  substantially  Increases.  However,  In 
this  case  the  added  losses  of  energy  are  introduced  during  its 
transfer  from  the  first  contour  to  the  second,  as  a  result  of 
which  the  effective  efficiency  is  lowered. 


Thus,  total  efficiency  Increases  (and  consequently,  the  cost- 
effectiveness  of  a  DTRD  is  increased)  when  an  increase  in  the 
thrust  efficiency  substantially  predominates  a  reduction  in  the 
effective  efficiency,  or  when  as  a  result  of  the  development  of 
the  DTRD  parameters  both  efficiencies  increase:  effective  and 
thrust . 


The  evaluation  of  the  effect  of  the  development  of  the 
operating  condition  and  bypass  ratio  parameters  on  the  efficiency 
of  prospective  DTRD  is  of  considerable  interest.  This  evaluation 
makes  it  possible  to  objectively  establish  the  prospects  for  an 
improvement  in  the  cost-effectiveness  of  DTRD  in  the  coming  years. 


Effect  of  the  Parameters  of  Operating 


Conditions  (tt*  and  TS)  and  Bypass 

H  J 


Ratio  (y)  on  the  Efficiency  of  (n 


nD  and  n  )  DTRD 
n  o 


(ttJJ,  TS)  and  bypass  ratio  (y)  on  the  specific  parameters  (R  and 
no  y,q 


C._),  and  also  on  efficiency  (ne>  nQ)  of  turbofan  engine  for 


standard  subsonic  design  conditions  (MQ  =  0.8  and  H  =  11  km). 


In  these  calculations  the  parameters  of  the  DTRD  changed  over 
a  wide  range  of  the  values: 


f,*  =  1200— I800CK; 
n*  =  20-  70; 
y  «  Q-j-20. 
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Figures  1-4  depict  the  results  of  the  calculations  according 
to  the  study  of  the  effect  of  the  parameters  of  operating  conditions 


- - - - *■-< 


Fig.  3*  The  effect 

specific  parameters 
KEY:  (1)  DTRD;  (2) 


ir*  and  Tt  on  the 
k  3 

of  the  DTRD  (y  =  20). 
kg/kg  h;  (3)  kg/kg  s. 


ip,  ?•* °*°  ***  Q*  J*  ia  ***  av 


4»  JV 


«?*>• 


Fig.  *1.  The  effect  of  the  parameters  of  the  operating  conditions 
with  y  =  0.5  and  y  =  10  on  the  efficiency  of  the  DTRD. 


The  following  improved  values  of  the  particular  efficiencies 
and  loss  factors  of  the  main  centers  of  prospective1  triple  shaft 
DTRD  were  accepted. 

Ti»*~  0,89;  =  r„  =  0,87: 

rir«(*,D*  =  =  =  0,92; 

=  *«c*  ~  0,97;  ^=0,99;  €,,“0.99. 

In  the  calculation  the  air  bleed  from  the  high  pressure 
compressor  for  the  cooling  of  turbine  blades  [VD]  (B/J)  was  also 
considered,  which  comprised,  dependent  on  the  level  of  the 
temperature  of  gases  T^,  from  4  to  10% . 

In  Fig.  4  in  coordinates  ne  and  nR  (with  the  applied  hyperbolic 

grid  of  the  isolines  p  =  const)  represented  the  total  effect  of 

o 

ir  *  and  Tf  at  fixed  values  of  y  (y  =  0.5  and  10)  on  the  efficiency 

J 

of  DTRD.  Such  a  coordinate  system  makes  it  possible  to  estimate 
in  demonstrative  form  how  development  (change)  of  the  parameters 
of  engine  simultaneously  affects  the  three  efficiencies. 

We  see  that  an  increase  of  the  bypass  ratio  from  0  [TJE]  (TPfl)2 
to  10  and  more  in  the  entire  range  of  parameters  of  operating 
conditions  sharply  increases  the  total  efficiency  because  of  the 
predominant  increase  in  the  thrust  efficiency  with  an  insubstantial 
reduction  in  the  effective  efficiency.  A  change  of  T§  and  it*  with 
y  *  const  in  dependence  on  the  numerical  value  of  the  initial 
parameters  of  engine  affects  its  efficiency  differently. 


lTaking  into  account  the  possibility  of  an  improvement  in  the 
hydraulic  and  gas-dynamic  perfection  of  the  air-gas  flow  area  of 
the  engine. 

2Ir  Fig.  4  y  =  0.5. 
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Prom  Pig.  4  it  is  possible  to  make  an  important  conclusion. 

It  is  advantangeous- to  produce  increase  tt*  only  up  to  a  value  of  n 

which  differs  little  from  the  maximum  (n  ) .  If  with  low 

e  MfiK  c 

bypass  ratio  the  increase  of  it*  over  tt#,  s  still  can  somewhat 

K  K  \  X\  ) 

e«Manc 

increase  n  »  then  with  large  y  a  similar  increase  of  u*  leads  to 

O  K 

a  drop  of  nQ*  Thus,  an  approach  to  the  selection  of  the  computed 
value  of  ttJ  at  large  values  of  y  should  be  analogous  to  that 
made  at  high  pressure  turbines. 

Another  conclusion  is  that  increase  of  T^,  as  a  rule,  makes 
the  cost-effectiveness  of  DTRD  worse  at  all  values  of  y  in  question, 
with  the  exception  of  very  high  values  of  tt*  ( ~ 7 0 ) . 

H 


Fig.  5.  The  effect  of  bypass  ratio  on 
the  efficiency  of  DTRD. 
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Figure  5  shows  the  effect  of  bypass  ratio  (for  fixed  values 
of  and  it*)  on  the  efficiency  of  DTRD.  We  see  that  an  increase 
of  y  from  C  to  20  sharply  increases  thrust  efficiency  from  0.34-0.44 
to  0.83-0.90,  i.e.,  more  than  2  times.  Effective  efficiency  in 
this  case  initially  decreases  somewhat  from  0.51-0.62  to  0.43-0.52, 
i.e.,  to  approximately  20$;  it  is  characteristic  that  in  the  range 
of  values  y  =  4-20  effective  efficiencies  virtually  retain  an 
invariable  value. 

The  total  efficiency  of  a  DTRD  with  an  increase. of  y 
continuously  grows  from  0.21-0.23  to  0.37-0.43,  i.e.,  almost  to 
80-100$,  whereupon  it  is  especially  intense  using  the  high 
parameters  of  operating  conditions  (Tl  =  l800°E  and  it*  =  70). 

J  K 


Fig.  6.  Fig.  7. 

Fig.  ,6.  The  temperature  effect  of  the  gas  before  the  turbine 
on  the  efficiency  of  DTRD. 

Fig.  7.  The  effect  of  the  compression  ratio  of  the  compressor 
on  the  efficiency  of  DTRD. 
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Figure  6  shows  the  temperature  effect  of  the  gas  before  the 
turbine  (for  a  series  of  fixed  values  y  and  tt* )  on  engine 

H 

efficiency.  We  will  see,  that  an  increase  continuously  lowers 
thrust  efficiency  (losses  with  outlet  velocity  grow).  It  is 
characteristic  that  a  drop  in  the  thrust  efficiency,  as  a  rule, 
compensates  for  an  increase  in  the  effective  efficiency.  In 
summation,  with  an  Increase  of  T{^  the  total  efficiency  falls  with 
the  exception  of  high  values  of  y  (y  ■  10)  and  tt*  (ir*  *  70°)  j 
however,  even  in  this  case  the  economic  value  of  TJJ  (to  which 
corresponds  the  maximum  of  total  efficiency  and  minimum  Cy^)  does 
not  exceed  a  value  of  l400°K. 

Finally,  Fig.  7  depicts  the  effect  of  the  compression  ratio  of 
compressor  n*  (for  a  series  of  fixed  values  y  and  T^)  on  the 
efficiency  of  DTRD.  From  this  figure  it  follows  that  with  an 
increase  of  tt*,  dependent  on  the  level  of  parameters  y  and  Tt, 
effective  efficiency  can  be  changed  differently:  with  low  y  and 
it  continuously  grows;  with  large  y  and  low  it  continuously 
falls.  Curves  n  can  have  a  maximum  on  it*. 

0  K 

The  total  efficiency  of  the  engine  with  an  increase  of  tt* 

H 

increases  at  high  values  of  and  has  a  maximum  with  tt*  j  40, 
when  Tjj  »  1200°K. 

Conclusions 

During  an  observed  high  level  of  the  efficiency  of  units  of 
prospective  DTRD,  calculated  for  the  large  subsonic  speeds  of 
flight,  to  provide  for  the  best  cost-effectiveness  of  the  engine: 

1.  It  is  advantageous  to  increase  the  bypass  ratio  of  engine 
to  values  of  y  *  10-15,  whereupon  the  larger,  the  higher  the 
accepted  level  of  T^  and  ttJ. 
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2.  It  is  inadvisable  to  increase  the  temperature  of  the  gas 
before  the  turbine  in  cruising  flight  over  values  of  1400°K. 

This  is  profitable,  however,  from  the  viewpoint  of  an  increase 

in  the  specific  thrust  of  a  DTRD,  i.e.,  for  reduction  in  the 
overall  dimensions  and  weight  of  engine. 

3.  It  is  advisable  to  continuously  increase  tt*  (up  to 

rr*  *  70-100)  only  at  high  values  of  gas  temperature  before  the 
turbine  (when  T»  >  1400-1500°K) .  When  y  =  10  and  =  1200°K  the 
economic  value  of  tt#  does  not  exceed  a  value  of  35-40. 

K 

4.  The  realization  of  a  DTRD  with  parameters  of  operating 
conditions  in  cruising  flight  (MQ  *  0.8;  H  =  11  km);  T^  *  1400°K; 
tt*  =  40  and  y  =  10  makes  it  possible  to  create  a  third  generation 

n 

DTRD  with  specific  fuel  consumption  C  ~  0.50-0.52  kg/kg*h  at 
values  of  efficiency  nQ  =  0.38-0.40;  nR  =  0.82-0.84;  ng  ~  0.46-0.48. 

Bibliography 

I.  K.isiiii'H  A  ,1.  Tcojmih  e<>).niiiHo-|>c..MitBiihi\  jmiraic'ieft.  «M; 

.'V.  1969. 

i  Hpo6.icMu  npiiMciicHiiH  ,ib\  \KOHiypiiux  TPH  b  rpjSKaaH.  koA  .ibiibuiii 
ll.u  pea.  iipcxj).  A.  .1.  K.ibmkiimb.  Tp\  iu  PKIIIII.X.  man.  Ilj.  1967. 


9 


APPROXIMATION  Or*  THE  BASIC  THERMODYNAMIC 
FUNCTIONS  USED  FOR  THE  CALCULATION  OF 
GAS  TURBINE  ENGINES 


M.  P.  Budzinauskas ,  V.  P.  Labendik 


In  this  article  equations  are  presented 
which  approximate  the  thermodynamic  functions 
used  for  the  calculation  of  gas  turbine  engines 
and  the  procedure  of  their  application  is 
presented.  Tables  are  presented  which  make  it 
possible  to  estimate  the  accuracy  of  the  approx 
imation  of  the  basic  thermodynamic  functions. 


Principal  Notations 


T  - 
AT  - 


i  - 

S  - 

AS  - 


TT  - 


5 


m  - 

T 

HC  .” 

Hu  ' 
a  - 


temperature 

change  in  temperature 

relative  temperature 

enthalpy 

entropy 

change  in  entropy 
compression  ratio  (expansion) 
the  relative  consumption  of  fuel 
combustion  completeness  coefficient 
fuel  heating  value 
the  excess  air  ratio 
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At  present  gas  turbine  engines  hold  the  ruling  position  in 
aviation.  In  accordance  with  this  the  front  of  the  scientific 
work  in  the  region  of  study  of  these  engines  was  considerably 
expanded.  The  complexity  of  GTE  of  latter  generations  makes  their 
calculations  more  laborious.  It  is  logical,  that  the  acute 
necessity  of  the  rapid  and  precise  thermodynamic  calculation  of 
these  engines  arises.  For  the  acceleration  of  calculations, 
computers  are  being  ever  more  widely  applied.  The  necessity  of 
more  precise  methods  of  thermodynamic  calculation  is  especially 
strongly  perceived  in  the  case  of  investigation  high-temperature 
DTRD  with  high  compression  ratios  and  bypass  configuration.  This 
problem  is  no  less  acute  during  studies  of  engines  for  supersonic 
aircraft . 

One  ought  to  consider  the  method  of  calculation  of  gas  turbine 
engines,  based  on  the  application  of  thermodynamic  tables,  most 
adequate  for  accuracy  and  simplicity.  Such  tables  on  the  basis  of 
the  newest  data  on  the  specific  heat  of  gc ses  were  compiled  in 
1956  by  V.  M.  Dorofeyev.  The  tables  give  the  thermodynamic 
functions  for  air  and  combustion  products  of  standard  kerosene. 

The  method  of  calculation  of  DTRD  by  application  of  tables  of 
thermodynamic  functions  can  be  easily  used  also  for  machine 
calculation,  however  in  this  case  it  is  necessary  to  encumber  the 
memory  unit  of  the  computer  by  a  large  quantity  of  tabulated  data. 
The  subsequent  selection  of  these  data  and  their  unavoidable 
interpolation  substantially  complicate  the  program  and  increase 
calculation  time.  It  is  completely  understandable  that  such  an 
approach  is  acceptable  in  the  case  of  applying  larg^  computers  of 
the  "Minsk-22"  type. 

In  the  case  of  the  application  of  small-scalp  omputers  of 
the  type  of  "Promin'"  etc.,  it  is  necessar'  >  '  er  to  simplify 
procedure  with  a  loss  of  accuracy,  or  to  itcuth  for  adequate 
methods  of  approximation  of  thermodynamic  functions  by  equations. 


The  proposed  method  for  the  calculation  of  QTE  is  based  on 
approximation  by  polynomials  of  the  basic  thermodynamic  functions. 
Such  an  approach  maximally  decreases  the  necessity  to  have  large 
immediate  access  memory,  since  the  relatively  small  quantity  of  the 
initial  values  and  coefficients  is  easily  introduced  directly  from 
a  console. 

During  the  calculation  of  the  processes  of  compression  and 
expansion  for  a  GTE  It  is  necessary  and  sufficient  to  have  the 
following  analytical  expressions  of  thermodynamic  functions  for 
air  (a  *  °°)  and  pure  combustion  products  (a  =  1) 

/=/(D;  9(0  and  15=  *(A7;  T), 

and  for  the  calculation  of  combustion  chambers  (basic  and  booster)  - 
cofunction  f (T*) . 

The  analytical  expressions  of  thermodynamic  functions  i  =  f(T) 
and  T  =  <p(i)  were  btained  by  means  of  the  approximation  of  the 
tabulated  dat  i  polynomials  of  the  form: 

i  =  aT*  +  bT*  +  cT*+dT  (1) 

and 

r=,/«  +  ?/"-pY,i  +  ;/.  (2) 

Approximation  of  thermodynamic  functions  by  polynomials  of  the 
fourth  degree  makes  it  possible  to  attach  equations  to  five  points 
of  the  tabular  values  of  these  functions. 

The  coefficients  of  equations  (1)  and  (2)  were  determined  at 
reference  points  by  the  solution  of  systems  of  linear  algebraic 
equations  of  the  fourth  order  relative  to  these  coefficients.  The 
solution  of  the  systems  of  equations  was  carried  out  on  the 
"Promin"'  computer.  The  values  of  the  coefficients  of  equations 
(1)  and  (2)  respectively  for  air  and  combustion  products  (a  =  1) 
are  given  in  the  following  table: 


Table  1. 


Nv  Coeffl- 
^S^cients 

a 

b 

c 

d 

Air 
a  ■  00 

-0.10353-10"10 

0.35002-10"7 

-0 .15931 • 10*^ 

0.24089 

Combustion 
products 
a  *  1.0 

-0.71294* 10”11 

0.21554-10-7 

0.12027*10*^ 

0.24163 

N.  Coeffi- 
^\Cients 

a 

8 

Y 

Air 
a  =  <*> 

0.90875-10*8 

-0.7331-10*5 

0. 5910-10*3 

4.1765 

Combustion 
products 
a  =  1.0 

0 . 25872 • 10~8 

-0.12449-10*5 

-0.1451*10~2 

4.1613 

Equation  of  the  adiabatic  curve  is 

AS 

„=/*  <3) 

in  coordinates  of  tt,  and  S  gives  the  simple  dependence  of  the 
compression  ratio  (expansion)  on  a  change  of  entropy,  and  in  the 
last  analysis  on  the  temperatures  of  the  beginning  and  end  of  the 
processes  of  adiabatic  compression  or  expansion. 


Let  us  find  dependence  of  AS  on  the  temperatures  of  the 
beginning  and  end  of  the  adiabatic  processes 


1S  = 


(4) 
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In  accordance  with  the  form  of  the  equation  1  ■  f(T),  the 
dependence  of  the  true  values  of  specific  heats  on  the  temperatures 
Is  written  by  the  equation 


C'^AT'  +  ar  +  cT+D 


(5) 


As  a  result  of  the  Integration  of  expression  (4)  taking  into 
account  equation  (5)  we  obtain,  that 


AS  =  5  (r„»  -  7i») + *-  (V  -r1*)+c(rt-  rt)  +  Din  £ . 


(6) 


It  is  not  difficult  to  demonstrate  that  the  coefficients  A,  B,  C 
and  D  unambiguously  depend  on  the  coefficients  of  a,  b,  c  and  d 
of  equation  ( 1) . 


Actually  , 

T 

i~aTi  +  hr  +  cr+,1T=  j  cfdT. 

u 

After  replacement  of  c^  in  integrand  by  its  value  from  equation  (5) 
we  obtain  the  equality, 

a  p  +  ar»  +  c  r  +  <tr=  £  p  - f- 1  p  +  £  r*  +  or, 

from  which  it  follows  that 

.4  =  4a;  5  =  3ft;  C^=2c  and  D  —  d. 

Then  expression  (6)  is  finally  written: 

AS  =  J  a  (7  V  -  r,»)  +  |ft  ( r,a  -  r,») + 2c(  ra  -  r.)  +  ./in  £».  ( 7 ) 

i.e.,  they  obtained  thermodynamic  function  AS  =  ^(AT;  T),  which 
makes  it  possible  to  determine  adiabatic  compression  or  expansion 
ratio  by  equation  (3)  when  the  temperatures  at  the  beginning  and 
end  of  these  processes  are  known. 

Equations  (3)  and  (7)  can  be  utilized  directly  for  determining 
the  dynamic  compression  ratio  ir  at  the  input  device  of  the  engine 
and  the  expansion  ratio  of  gas  in  the  turbine,  tt*,  since  the 
initial  and  final  temperatures  for  these  processes  are  unambiguously 
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determined.  It  i  s  somewhat  more  complex  to  determine  the 
temperature  at  the  end  of  the  process  of  compression  in  compressor 
and  the  thermodynamic  temperature  in  a  jet  nozzle,  since  in 
these  cases,  as  a  rule,  they  are  known  as  tt*  and  tt  .  Let  us  show 

H  p  C 

how  to  use  equations  (3)  and  (7)  to  determine  the  temperature  at 
the  end  of  the  adiabatic  process  of  compression  (T!$a^)  in  the 
example  of  the  calculation  of  a  compressor.  If  tt*  is  given,  then 
according  to  equation  (3)  we  determine  a  change  in  entropy  AS  . 
Further  for  determining  temperature  T!jafl  with  the  assigned  in 
advance  accuracy  equation  (7)  we  solve  by  the  method  of  successive 
approximations.  In  the  first  approximation  of  T2afl*  we 
according  to  the  formula 


where 


/  • 

T  *»  —  T#  i  Lam 
—  'l 


AM,*=102.5V(r.*M*-|). 


That  found  in  the  first  approximation  of  T we  substitute  in 
modified  equation  (7),  which  we  will  rewrite: 

f= ~  « ( -  /',*»)  + 1  b  (Tux“  -  Tn + 2c  (r,./-  r*),  < 


where 


F  =\S~d  In  rQ  - 
•\ 


i.e.  , 


From  expression  (9)  we  find  the  new  value  of  temperature  T£a^, 


lSa-F 


r* _ T  *  0  * 

214  ~  ‘  l  “ 


The  found  value  of  Tt  we  again  substitute  in  (8).  The  process  of 

<-&A 

approximation  is  continued,  when  the  difference  TJ|”  -TJj^-  in 
absolute  value  does  not  become  equal  to  or  less  than  the  preassigned 
value  (for  example  0.1°). 
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The  thermodynamic  temperature  in  Jet  nozzle  is  analogously 


determined,  if  tt  is  known. 
*  pc 


For  determining  the  relative  consumption  of  fuel  mT  both  in 
basic  and  afterburners  it  is  most  convenient  to  use  the  formulas 
given  in  [1]: 

a)  basic  combustion  chamber 


"V 


k* 


*«c 


do) 


b)  afterburner 


-Jil±  mrl /(V)~/(V)) 
5*.c^-/(V)  +  V  ' 


(11) 


where  the  enthalpies  i}$,  i*,  i|J  and  i*  are  determined  for  the 
approximate  temperatures  on  air. 


The  values  of  functions  f(T*)  in  expressions  (10)  and  (11)  are 
given  in  the  table  located  in  [1].  The  tabular  values  of  functions 
of  f(T*)  are  approximated  with  a  sufficient  degree  of  accuracy  by 
the  equation 

f{T)^AT-\-  flP-f-cr-f  D,  (12) 

where  A  =  -0. 33025’ 10"7 
B  =  0.25084-10”3 
C  =  0.35186 
D  =  -17.533. 


To  evaluate  the  accuracy  of  the  approximation  of  the  basic 
thermodynamic  functions  by  equations  (1),  (2),  (7)  and  (12),  in  the 
following  tables  precise  (tabular)  and  computed  values  of  these 
functions  are  given. 


Table  2. 


I  Pete  rained  by  formulas  fll,  12 


203.16 
mi)i 

603.16 

803.16 
1003, H> 

1203.16 

1403.16 

1603.16 

1813.16 


W) 

r- 

rm 

I. 

‘r 

T. 

Tt 

48,557 

203,16 

96,518 

100,60 

402,87 

402,81 

145,81 

153,91 

602,73 

6022*9 

197,01 

210,03 

803,80 

803,98 

250,18 

269,04 

1004,1 

1004.6 

306,05 

330,79 

1203,2 

12036 

363,22 

394,65 

1400,3 

1401,7 

421.08 

460,01 

1600,3 

1600,6 

478,34 

525.82 

1810,5 

1804,2 

The  comparison  of  precise  values  of  the  basic  thermodynamic 
functions  with  the  values  of  these  functions,  calculated  from 
formulas  (1),  (2),  (7)  and  (12),  shows  that  in  the  temperature 
range  in  question  (see  Tables  2,  3  and  4)  the  deviation  of  the 
calculated  values  from  their  precise  values  does  not  exceed  0.5?. 
Keeping  in  mind,  that  tabular  (precise)  values  of  the  basic 
thermodynamic  functions  in  question  are  determined  with  an.  accuracy 
of  0.5%  (see  V.  M.  Dorofeyev's  work  "Thermodynamic  Tables  for  the 
Calculation  of  Gas  Turbine  Engines"),  then  it  is  evident  that  the 
error  obtained  during  the  calculation  of  these  functions  according 
to  formulas  (1),  (2),  (7)  and  (12)  does  not  exceed  the  limits  of 
accuracy  of  the  available  tables  of  the  basic  thermodynamic 
functions . 

Good  convergence  of  results  of  the  approximation  of  thermo¬ 
dynamic  functions  by  the  equations  presented  above  makes  it  possible 
to  apply  formulas  (1),  (2),  (7)  and  (12)  with  a  sufficient  degree 
of  accuracy  for  the  thermodynamic  calculation  of  gas  turbine  engines. 


•*> (wwBartwwtow, _ 


Table  3.  AS  -  ♦(AT;  T2) ,  where  AT  -  T0-T^ ,  and  T*  -  203.16°K: 
Tf  «  403. 16°K.  2  1  1b 


ir 

200 

400 

600 

800 

1000 

1200 

!  1400 

1600 

1800 

— 

2000 

8 

II 

M 

-=f 

II 

* 

f2^ 

=  ft 

52 

►"  R 

0,16447 

— 

0,10710 

0,26352 

0,33699 

0,39648 

0,44670 

0,35811 

0,49042 

0,40155 

0.52W& 

0.44056 

0,56346 

049470 

0,82325 

0,18764 

0.25333 

040928 

0.47603 

040848 

A.S’„  .  s/>ui> 

?. 

II 

*1 

a. 

** . 

3N 

«e-w 

0 

r— 

0,16420 

0,26340 

0.33650 

0,30620 

0,44660 

0.49050 

042890 

046260 

0401 50 

O.6IM0 

_ _  __ 

~  TIT'"  '  ' 

Avlr  hKtU/KC  .  CfHit) 

». 

II 

n 

0,10720 

0,18750 1  0,25320 

1 

0,30900 

_ 

0,3^2o!(),40HH) 

1 

0.44040 

0,1745 

04042 

__  _ 

(3) 


KEY:  (1)  kcal/kg*deg;  (2)  Precise  values;  (3)  According  to 

formula  (7). 


Table  4. 


T*°  K 

40<) 

6U0 

800 

10.  0. 

uoo 

140) 

WkiO 

1800 

2000 

f  (T*)  J  3iia'ifiniM 

161,23 

2764« 

4U/.60 

552,02 

70840 

876.11 

1052J6-1 

12.3642 

1425.4* 

r  rr*\  i  Piicci.  no  4>i>pwy. 
J  '  ’  j^2)  .ic  (12) 

161,23 

276.76 

40749 

_ 

552.14 

_ 

708,84 

876,10 

_ _ 

1052,4 

1236.0 

1425,4 

KEY:  (1)  Precise  values;  (2)  Calculation  according  to  formula  (12). 
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METHOD  FOR  THE  CALCULATION  OF  DTRD 
WITH  HEAT  RECOVERY 

M.  P.  Budzinauskas 

The  work  gives  the  procedure  for  approximate 
computation  of  the  parameters  of  efficiency  of 
DTRD  with  heat  recovery  in  the  first  duct. 

Principal  Notations 

T  -  temperature 
p  -  pressure 
i  -  enthalpy 
Ai  -  change  in  enthalpy 
q  -  heat 

c  -  rate  of  flow  of  gas  (air)  stream 
v  -  flight  speed 
y  -  bypass  ratio 

tt  -  compression  (expansion)  ratio 
G  -  flow  rate  per  second 
9  -  velocity  coefficient 
E  -  coefficient  of  completeness  of  combustion 
Hu  -  fuel  heating  value 
Cp  -  specific  heat  at  constant  pressure 

LQ  -  quantity  of  air  theoretically  necessary  for  complete 
combustion  of  1  kg  of  fuel 
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R  -  universal  gas  content 
k  -  adiabatic  index 


Indices 
Superscript : 

*  -  stagnation  parameter 
Subscripts : 

I  -  parameter  in  the  first  contour  of  engine 
II  -  the  parameter  in  the  second  contour  of  engine 
2,  3,  ^  -  the  parameters  in  characteristic  cross  sections  of 
engine  (see  Pig.  1) 
b  -  parameter  for  air 

r  -  parameter  for  gases  or  pure  combustion  products 
p  -  regeneration 
t  -  fuel 
pc  -  jet  nozzle 

Abbreviations 

[DTRD]  (ATPH)  -  turbofar  engine 
[TVD]  (TBfl)  -  turboprop  engine 
[KPD]  (HPA)  -  efficiency 

Introduction 

The  successful  solution  of  the  proolem  of  creation  of  an 
aircraft  of  large  passenger  capacity  and  flying  range  within  the 
framework  of  acceptable  profitableness  requires  an  essential 
reduction  in  the  specific  fuel  consumption  of  the  power  plants. 


The  selection  of  engine  type  first  of  all  affects  the  level 
of  specific  fuel  consumption,  however,  for  a  passenger  aircraft 
the  possibilities  of  such  a  selection  are  limited.  In  practice, 
selection  falls  exclusively  to  DTRD. 

A  reduction  in  the  specific  fuel  consumption  of  DTRD  may  be 
attained  by  an  increase  in  the  compression  ratio  of  air  in  the 
first  contour  of  an  engine  with  a  simultaneous  increase  in  its 
bypass  ratio  and  gas  temperature  in  front  of  the  turbine.  The 
transfer  to  higher  gas  temperatures  in  front  of  the  turbine  makes 
it  possible,  to  a  great  extent,  to  maintain  sufficiently  high 
values  of  the  specific  thrust  of  the  engine. 

The  possibilities  of  reduced  specific  fuel  consumption  by 
increased  efficiency  of  the  basic  engine  components  are  very 
limited,  and,  at  the  present  time,  exhausted  to  e  great  extent. 

Besides  the  mentioned  factors  which  facilitate  an  increase 
in  the  cost-effectiveness  of  DTRD,  the  so-called  regenerative 
method  of  reduction  in  the  specific  fuel  consumption,  which  is 
widely  applied  for  stationary  gas-turbine  power  plants,  is  of 
great  interest.  Successful  attempts  at  the  application  of  a 
regenerative  cycle  for  TVD  are  known,  [1],  [2]  and  [3]. 

The  essence  of  regenerative  cycle  consists  in  the  fact  that 
the  preheating  of  air  entering  the  combustion  chamber,  the  heat 
of  the  gases  which  escape  from  the  first  contour  of  the  DTRD  is 
partially  utilized.  Thus,  the  degree  of  the  preheating  of  air 
in  the  combustion  chamber,  and,  accordingly,  also  fuel  consumption 
is  reduced. 

The  schematic  diagram  of  DTRD  with  heat  recovery  is  depicted 
on  Fig.  1.  The  reaction  gases  of  the  basic  contour  are  guided  to 
the  heat  exchanger  arranged  directly  after  the  turbine  where  part 
of  the  heat  is  given  up  to  the  compressed  air  which  enters  the 
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combustion  chamber  from  the  compressor.  Reaction  gases,  in  their 
turn,  passing  through  the  heat  exchanger,  are  cooled. 


a) 


o 

Fig.  1.  Schematic  diagram  of  a  turbofan 
engine  with  heat  recovery. 

KEY:  (1)  Heat  exchange. 


Method  of  Calculation 


The  gas-dynamic  and  thermal  design  of  DTRD,  which  work  on  a 
regenerative  cycle,  is  performed  into  two  stages.  During  the 
first  stage  the  parameters  of  gas  in  the  basic  sections  of  the 


turbo-compressor  part  of  the  engine  and  the  relative  consumption 
of  fuel  (m  =  G  /G  , )  are  calculated  according  to  the  usual 

T  T  B  j. 

procedure  with  the  use  of  the  well-known  equations  and  formulas 


of  the  theory  of  turbofan  engines  [ 4 ] .  The  findings  are  accepted 
as  initial  and  utilized  in  the  second  stage  for  determining  the 
basic  parameters  of  DTRD  with  heat  recovery. 

The  regenerative  cycle  is  considered  given,  if,  besides  the 
basic  initial  parameters  and  the  coefficients  of  the  particular 
losses  of  the  usual  cycle  of  a  DTRD  the  degree  of  regeneration  of 
heat  Hp  is  known  and  total  pressure  losses  o*  along  air  and  gas 
circuits  of  the  heat  exchanger  are  known. 


By  the  degree  of  regeneration  of  heat  is  meant  the  relationship 
of  the  actual  preheating  of  air  T2p”T2  in  the  reSenerator  to  the 
theoretically  possible  preheating  Tjj-T*,  i.e., 


R 

v 


(i) 
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The  total  pressure  losses  along  the  air  and  gas  circuits  of 
the  regenerator  are  estimated  by  the  total  pressure  recovery 
coefficient . 

V  =  (2) 

where  o*  tota^  pressure  recovery  coefficient  in 

the  air  circuit  of  the  regenerator  system;  o*  ^  =  p  Hp/pjj  -  the 
total  pressure  recovery  coefficient  in  the  gas  circuit  of  the 
regenerator. 


On  the  basis  of  the  findings  of  calculating  the  turbo¬ 
compressor  portion  of  the  Initial  DTRD,  the  further  calculation 
of  the  parameters  of  efficiency  of  an  engine  with  heat  recovery 
is  conducted  from  the  determination  of  the  relative  fuel 

consumption  m  (m  =  G  /G, ,  where  G  is  conducted  the  fuel 
rp  rp  Tp  1*  Tp 

consumption  per  second  of  a  DTRD  with  heat  recovery). 

For  determining  the  relative  fuel  consumption  let  us  write 
the  equation  of  energy  for  two  cross  sections  (at  compressor  and 
the  combustion  chamber  outlet)  of  the  first  contour  of  a  DTRD. 


(1  +  mrp  -  mu„)t*  =  (1  -  mu,  y)  »,*  4  (1  -  m0„)  Aip*  + 

+  0  -  )<7,  +  n„i0, 


where 


The  latter  equation  is  solved  in  conjunction  with  the  equation 
of  heat  balance  in  the  combustion  chamber  of  the  engine 

He  "u  ^ip  ~~  ( 1 

Bearing  in  mind,  that 


, * _  (i  +/.„)/*, P  .  ,  r .  0  4- /,.) Wjji .  . 

*  1  +  «,p  - *  +|l  1  8*  * 


where  i^r  and  i^B  -  correspondingly  the  total  enthalpies  of  air 
and  pure  combustion  products  at  the  entry  to  turbine;  after  simple, 
but  bulky  transforms  we  will  obtain 

From  expression  (1)  it  follows  that 

v-  V  =  (/■«• -V)v 


Then  finally 


«rp  =  (l  -  Tfh  Zff 


(3) 


For  determining  the  specific  thrust  of  the  first  contour  of  a 
DTRD  with  heat  recovery  we  find  the  exhaust  gas  velocity  from  the 
nozzle  of  this  contour: 


—  7pd 


/ 


k 

f'--M 

2*i±  i'?,v 

(4) 


The  temperature  Tjjp  can  be  found  from  the  equation  of  heat 
flow  through  the  heat-transmitting  surfaces  of  the  regenerator. 


Since 


then 


C,r(l  =  T,*). 

v-  r*  - =\{Ti*  -V). 
VaV-VPi'-W 


o-Cpi(l-f  m,  )V 


25  • 


where 


(5) 


Finally,  the  degree  of  expansion  of  gases  in  the  jet  nozzle 
of  the  first  of  a  DTRD  with  heat  recovery  we  find  according  to  the 
formula: 

sp  '  Rp<i**  ( 6 ) 

In  formula  (6)  is  the  degree  of  expansion  of  gases  in  the 

pc  1 

Jet  nozzle  of  the  first  contour  of  the  initial  DTRD. 

The  specific  thrust  of  the  first  contour  of  a  DTRD  with  heat 
recovery  will  be 


Then  the  specific  thrust  of  a  DTRD  with  heat  recovery  is 
written : 

O  _  .Ip  +  .V#,all . 

- —  (7) 

Finally,  the  specific  fuel  consumption  of  a  DTRD  with  heat 
recovery  is  determined  from  the  formula: 

C,;f  =  3600  (8) 

^yji  I* 

Thus,  formulas  (3-8)  at  given  values  of  n  and  o*  make  it 

P  P 

possible  to  accomplish  the  successive  calculation  of  the  basic 
parameters  of  the  efficiency  of  a  turbofan  engine  with  heat 
recovery  according  to  the  available  data  of  initial  DTRD. 
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STUDY  OF  THE  EFFICIENCY  OF  A  DTRD 

WITH  HEAT  RECOVERY 

M.  P.  Budzinauskas 

In  the  article  are  given  some  results  of 
the  investigation  of  the  basic  technical  and 
economical  indices  of  DTRD  with  heat  recovery 
and  long-distance  main-line  aircraft  equipped 
by  engines  of  this  type. 

Principal  Notations 

-  gas  temperature  before  the  turbine 
7T#  -  the  compression  ratio  of  air  in  compressor 
tt*  -  the  compression  ratio  of  air  in  ventilator 
y  -  a  bypass  ratio 
n*  -  compressor  efficiency 
n*  -  the  efficiency  of  ventilator 
n*  -  the  efficiency  of  turbine 
np  -  regeneration  effectiveness1 

a*  -  total  pressure  recovery  factor  in  air  and  gas  circuits 
of  the  regenerative  system 

o|J  -  total  pressure  recovery  factor  in  input  device 

B  X 

“See  page  23  of  this  collection. 


o*  -  total  pressure  recovery  factor  in  combustion  chamber 

H  C 

£  -  completeness  of  combustion  coefficient 

9  -  the  velocity  coefficient  of  Jet  nozzle 

-  specific  fuel  consumption 
Ry^  -  specific  thrust  of  the  engine 

Gp8r  "  “  relative  weight  of  the  regenerative  system, 

where  Gper  -  weight  of  regenerative  system; 

G  j  -  the  flow  rate  per  second  of  the  air 
through  the  first  contour  of  the  engine 
H  -  flight  altitude 
M  -  Mach  number  of  flight 

Indices 

Superscript : 

*  -  stagnation  parameter 
Subscript : 

I  -  parameter  in  the  first  contour  of  engine 
II  -  parameter  in  the  second  contour  of  engine 
a k  -  economic 

Contractions 

[DTRD]  UTPA)  -  turbofan  engine 
[GTE]  (n\n,)  -  gas  turbine  engine 

Introduction 

The  creation  of  subsonic  main-line  aircraft  with  a  flying 
range  of  up  to  15,000  km  is  possible  in  the  presence  of  highly 
economical  engines.  The  specific  fuel  consumption  of  such  engines 
at  a  cruise  setting  (H  =  11  km  and  M  =  0.8)  should  be  found  within 
the  limits  from  0.5  to  0.53.  Such  a  high  cost-effectiveness  of  a 
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DTRD  can  be  achieved  either  by  an  increase  in  the  compression  ratio 
of  air  in  the  first,  contour  of  the  engine  with  a  high  bypass  ratio 
and  gas  high  temperature  before  the  turbine  or  by  use  of  a  DTRD  with 
heat  recovery. 

In  this  article  are  given  some  results  of  the  investigation 
both  of  DTRD  with  the  heat  recovery  and  the  long-distance  main-line 
aircraft  equipped  with  these  engines. 

The  regenerative  system  of  a  DTRD  with  heat  recovery  always 
includes  a  heat  exchanger  (see  Pig.  1  in  the  foregoing  article).1 
For  a  transport  GTE  heat  exchangers  both  of  the  recuperative  type 
and  periodic  action  with  a  rotating  matrix  ("rotating")  are  used. 

The  "rotating"  heat  exchangers  differ  by  the  great  compactness  of 
the  heat-transmitting  elements,  however  the  large  weight  of 
auxiliary  parts  Impedes  a  creation  of  sufficiently  light  (engines) 
for  application  in  aviation  constructions.  Furthermore,  the 
presence  of  seals  in  the  heat  exchangers  of  this  type  limits  the 
possibilities  of  use  of  elevated  temperatures  and  compression 
ratios  of  the  cycle,  and  the  unavoidable  overflow  of  air  into 
engine  cavity  after  the  turbine  lowers  the  effect  of  heat  recovery. 
Thus,  during  selection  of  the  fundamental  characteristics  of  the 
regenerator  np  and  ap,  the  use  of  tubular  construction  of  the 
recuperative  heat  exchanger  was  proposed  with  the  cold  air  duct 
within  the  tubes.  The  studies  of  recent  years  of  regenerators  of 
this  type  [1],  [2]  and  [3]  show  that  they  can  be  made  sufficiently 
lightly  and  compactly.  From  work  [3]  it  follows  that  it  is 
virtually  possible  to  obtain  a  degree  of  heat  recovery  of  np  =  0.75,' 

in  this  case  the  coefficient  of  retention  of  total  pressure  a*  of 

P 

approximately  0.91^  in  the  air  and  gas  circuits  is  obtained. 

The  values  n  and  a#  were  stated  as  the  basis  of  the  calculation 
P  P 

of  a  DTRD  with  heat  recovery. 


l3ee  page  23  of  this  collection. 


In  this  article  the  results  of  the  investigation  of  DTRD  with 
*  1600°K  over  a  wide  range  of  change  in  the  remaining  operating 
conditions  parameters  (tt*j  *  5-20,  y  ■  2-10)  are  represented,  with 
flight  conditions  H  =  11  km  and  M  *  0.8.  Thus  were  accepted  the 
following  values  of  the  coefficients  of  partial  losses: 

a,/ =  0,98;  i.t*  =  0,96;  ;\c  =  0,99; 

=  0,85;  1,,*  =  0,89;  Tlf*_0,9!;  ?pc,  =  -rPc„  =  0,99. 

The  determination  of  the  parameter  of  DTRD  efficiency  of  the 
regenerative  cycle  was  carried  out  according  to  the  methods  given 
in  the  foregoing  article.1 

Results  of  the  Study  of  DTRD  with 
Heat  Recovery 

Figures  1-4  depicts  the  calculated  dependences  of  the  specific 
fuel  consumption  of  a  DTRD  on  the  compression  ratio  of  air  in  the 
second  contour  tt*  and  on  the  bypass  ratio  y  at  different  values  of 
tt#  .  The  presented  graphic  dependences  make  it  possible  to 
visually  estimate  the  possibilities  of  a  reduction  in  the  fuel 
consumption  during  heat  recovery  and  to  optimize  engines  according 
to  compression  ratios  In  contours. 

The  comparison  of  those  dependences  given  in  Figs.  1-4  for  the 
simple  and  regenerative  cycles  of  a  DTRD  (using  identical  operating 
condition  parameters)  shows  that,  in  the  first  place,  the  cost- 
effectiveness  of  the  engine  substantially  improves  from  heat 
recovery;  in  the  second  place,  in  the  case  of  applying  heat  recovery 
u  ■  '.ler  values  of  economic  compression  ratios  in  the  second  contour 
ar  obtained. 

^See  page  23  of  this  collection. 
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Fig.  1.  Dependence  of  C  on  tt*  and 
y  with  tt*t  =  5.0.  y,q  B 

KEY:  (1)  DTRD  without  heat  recovery 

(2)  DTRD  with  heat  recovery. 
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Fig.  2.  Dependence  of  C  on  tt*  and 
y  with  tt*  =10.  yfl  B 

K  X 

KEY:  (1)  DTRD  without  heat  recovery 

(2)  DTRD  with  heat  recovery. 
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Fig.  3.  Dependence  of  C  on  tt*  and 
y  with  7T * T  =  15.  B 

K  X 

KEY:  (1)  DTRD  without  heat  recovery; 

(2)  DTRD  with  heat  recovery. 


KEY:  (1)  DTRD  without  heat  recovery; 

(2)  DTRD  with  heat  recovery. 
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For  example  from  the  graphs  of  Fig.  2  It  Is  obvious,  that 
with  y  ■  2,  the  reduction  of  because  of  heat  recovery  is 

31.51,  and  with  y  -  10,  C  drops  to  25.7*. 

J 


(1) 


Fig.  5.  Dependence  of  C  tt  and 

yA  B  3H 

on  (Results  of  the  opti¬ 

mization  of  a  DTRD  with  heat  recovery 
f or  C  ) . 

ya 

KEY:  (1)  DTRD  with  heat  recovery. 


33 


The  results  cf  the  optimization  of  the  DTRD  regenerative  cycle 

according  to  specific  fuel  consumption  depending  on  the  compression 

ratio  in  the  first  contour  and  on  the  bypass  ratio  are  given  In 

Pig.  5a.  These  graphic  dependences  were  constructed  according  to 

the  data  presented  in  Figs.  1-4  by  selection  of  minimum  values  of 

dependent  on  y  and  lrjjjj.  The  economic  values  of  compression 

ratio  in  the  second  contour  (Pig.  5c)  correspond  to  these  minimum 

values  of  C  .  Furthermore,  Pig.  5  depicts  the  corresponding 
y  a 

values  of  R  (c) . 
y 

The  analysis  of  the  results  of  the  optimization  of  engines 

for  Cy^  shows  that  the  economic  values  of  dependent  on  y  are 

located  in  the  interval  of  10-13.  The  exclusively  flat  minimum 

of  these  dependences  makes  it  possible  to  consider  that  for  all 

values  of  y  (from  2  to  10)  tr*  =  10.  To  select  lower  values 

is  inexpedient,  since  this  leads  to  a  sharper  reduction  of 

R  (see  Fig.  5c). 
y  ^ 

The  obtained  minimum  values  of  C  for  engines  with  heat 

ya  b 

recovery  with  tt*  =  10  (see  Fig.  5a)  in  cruising  flight  and  also 

H  X 

possible  to  achieve  for  DTRD,  that  work  on  a  simple  cycle,  by  a 
very  considerable  increase  in  the  compression  ratio  of  the  first 
contour  compressor. 

Figure  6  depicts  a  dependence  of  Cy^  and  on  ttJJj  for  a 
simple  DTRD.  The  comparison  of  the  graphs  presented  in  Fig.  5a  and 
b  shows  that,  for  example,  with  y  =  10,  the  level  of  C  of  a 

y  4 

simple  DTRD  becomes  the  same  as  of  the  engine  which  works  on 
regenerative  cycle,  only  when  the  compression  ratio  tt*  =  58,  but 
with  y  =  8,  the  identical  specific  fuel  consumptions  will  be  when 
tt*t  =  65.  The  fact  that  the  compared  engines  have  virtually 
identical  specific  thrusts  (see  Figs.  5c  and  6)  is  remarkable. 


Fig.  6.  Dependence  of  C  and  R  on  it* 
for  a  simple  DTRD .  yA  yA  Kl 


The  reduction  of  C  both  by  heat  recovery  and  by  an  increase 

J 

in  the  compression  ratio  tt*  unavoidably  leads  to  an  increase  in 
the  specific  weight  of  the  engine.  In  the  first  case  the  specific 
weight  of  the  engine  increases  due  to  the  weight  of  regeneration 
system,  and  in  the  second  -  due  to  a  gain  in  weight  of  the  turbo¬ 
supercharger  and  engine  block. 

The  approximate  estimate  of  the  specific  weight  of  the  compared 

engines  (having  identical  C  ),  shows  that  a  simple  DTRD,  with 
ivi  y  **' 

iTkI  =  °  and  y  =  has  a  specific  weight  up  to  16%  greater  than 

a  DTRD  with  heat  recovery  (it#  =  10;  y  =  10  and  G  =  G  /Gx  =  26). 

Hi  per  perl 

Results  of  the  Study  of  DTRD  Equipped 
Aircraft  with  Heat  Recovery 

The  study  of  the  efficiency  of  long-distance  aircraft  equipped 
with  DTRD  with  heat  recovery  is  of  great  interest.  These  studies 
were  carried  out  on  the  assumption  that  the  operating  conditions 
parameters  of  the  engines  were  selected  from  the  condition  of  ensuring 
minimum  fuel  consumption  in  cruising  flight.  DTRD  with  heat 
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recovery,  that  have  =  1600°K,  tt*^.  ■  10,  n*  *  1.59  and  y  «  10 
(see  Pig.  2)  correspond  to  the  condition  of  a  minimum  C  . 

J  ^ 

The  studies  of  the  parameters  of  aircraft  efficiency  were 
carried  out  according  to  the  methods  given  in  work  [4]. 

Since,  at  present,  reliable  data  on  the  structural 
efficiency  factor  of  heat  exchangers  with  the  increased  degree  of 
heat  regeneration  are  lacking,  during  the  calculations  the  relative 
weight  of  regenerator  (G  r)  was  varied  within  the  limits  from 
4. 0  to  26 . 

As  one  would  expect  with  an  increase  of  the  relative  weight 
of  the  regenerator  the  basic  parameters  of  the  efficiency  of  the 
aircraft  studied  deteriorate  somewhat.  For  example,  with  increase 
Cper  from  4.0  to  26,  the  takeoff  weight  and  the  prime  cost  of 
transport  with  a  given  flying  distance  respectively  increase  by 
0. 5  and  4 . 0  percent . 

In  the  following  table  are  presented  the  results  of  the 
study  of  the  basic  parameters  of  the  efficiency  of  aircraft  equipped 
both  simple  DTRD  and  those  with  heat  recovery.  In  this  case  the 
relative  weight  of  the  regenerator  was  taken  as  Gper  =  26. 

From  the  comparison  of  the  data  presented  in  the  table,  it  is 
apparent  that  according  to  the  basic  indices  (G  and  a),  aircraft 
which  have  DTRD  with  heat  recovery  are  indisputably  more  profitable, 
than  aircraft  equipped  with  simple  DTRD.  For  examole,  because  of 
the  application  of  heat  recovery,  GB3J1  for  aircraft  Nos.  1,  2,  3 
and  4  descends  to  2.1,  10,  3.7  and  10.7  percent.  In  this  case  the 
prime  cost  of  transport  correspondingly  drops  to  1.44;  9-75;  3.34 
and  10.4  percent.  A  more  noticeable  reduction  both  of  takeoff 
weight  and  prime  cost  of  transport  Is  observed  for  long-distance 
aircraft  No.  2  and  No.  4. 
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KEY:  (1)  Aircraft;  (2)  Indices;  (3)  Turbofan 
without  regen;  (4)  Turbofan  with  regen; 

(5)  Distance  L  (km);  (6)  Number  of  passengers 
(n);  (7)  Aircraft  takeoff  weight  G  ; 

03J1 

(8)  Freight  load  G  /£  ;  (9)  Power  plant 

G  /f  ;  (10)  Prime  cost  of  transport  A, 

cy  ■  t  cy 

Kop/t  km. 


In  conclusion  it  should  be  noted  that  at  present  on  .  the  way 
to  practical  lmpleme itatlon  of  DTRD,  that  work  with  a 
regeneration  cycle,  the  main  obstructions  are  the  unsatisfactory 
performing  characteristics,  reliability,  and  structural  efficiency 
factor  of  existing  small-scale  heat  exchangers. 

Conclusions 


1.  The  studies  carried  out  show  the  possibility  in  principle 

and  the  advisability  of  the  creation  of  DTRD  with  heat  recovery. 

When  T^  =  l600°K;  =  10,  =  0.75,  and  o#  =  0.91*1,  such  engines 

will  have  C  ~  0.5  kg/kg’hr. 

Y 
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2.  In  comparison  with  simple  DTRD,  the  minimum  specific  fuel 
consumption  for  DTRD  with  heat  recovery  is  reached  with 
substantially  lower  compression  ratios  of  air  in  the  first  contour 

s  10).  For  an  optimum  DTRD  with  heat  recovery  for  lower 
compression  ratios  are  also  required  in  the  second  contour  (see 
Figs.  1,  2,  3  and  4) . 

3.  The  application  of  DTRD  with  a  regenerative  cycle,  makes 
it  possible  to  raise  the  efficiency  of  long-distance  main-line 
aircraft  even  with  the  low  structural  efficiency  factor  of  the 
regeneration  system.  The  prime  cost  of  transport  in  the  case  of 
applying  engines  with  heat  recovery,  depending  on  the  calculated 
flying  range,  drops  from  1.44  to  10.4  percent  (see  table). 
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EFFECT  OF  THE  SM'.LL  DEVIATIONS  OF  OPERATING 
CONDITIONS  AND  LIGHT  CONDITIONS  PARAMETERS 
ON  THE  SPECIFIC  PARAMETERS  OF  DTRD 


V.  P.  Labendlk 


In  the  article  the  application  of  a  method 
of  low  deviations  ir  examined  for  the  analysis 
of  the  efficiency  and  cost-effectiveness  of  a 
twin- shaft  DTRD. 

Working  formulas  for  the  calculation  of 
coefficients  of  the  influence  of  operating 
conditions  and  flight  conditions  parameters 
on  specific  engine  characteristics  are  given. 

As  an  example  the  values  of  coefficients  of 
influence  are  given  with  concrete  data  in 
connection  with  twin-shaft  DTRD  and  the  results 
of  the  calculation  are  analyzed. 


Introduction 


When  selecting  the  basic  operating  conditions  parameters  of 
the  projected  engines  it  requires  finding  their  optimum  relationship 
for  the  assigned  flight  mode,  and  for  this  it  is  necessary  to  know 
how  and  to  what  degree  the  given  parameters  in  the  range  in  question 
affect  the  efficiency  and  the  cost-effectiveness  of  the  engine;  a 
change  in  which  of  them  causes  the  greatest  effect.  In  other  tasks 
it  is  required  to  obtain  specific  engine  characteristics  with  the 
same  initial  parameters,  but  for  several  distinct  flight  conditions; 


at  other  values  of  the  loss  factor,  efficiency  and  so  on.  In  this 
case  it  is  necessary  to  calculate  a  series  of  variants  which  differ 
by  a  small  change  in  the  initial  values. 


These  results  can  be  obtained,  and  without  resorting  to 
repeated,  sufficiently  bulky  calculations  of  the  engine  by  the 
usual  methods.  As  a  simpler  and  more  demonstrative  method  of 
solution  it  Is  convenient  to  use  the  method  of  small  deviations. 

When  using  this  method  it  is  necessary  to  conduct  the  gas- 
dynamic  calculation  and  to  determine  the  operating  conditions 
parameters  only  for  the  initial  system  in  question.  The  use  of 
relationships  of  the  low  deviations  method  actually  means  the 
linearization  of  the  initial  equations  of  the  process  which 
considerably  simplifies  the  analysis  of  the  dependences  between 
the  increments  of  the  connected  parameters. 

The  system  of  fundamental  equations  in  small  deviations  is 
obtained  by  differentiation  of  the  usual  equations  of  [GTE]  (TT/l,) 
operating  conditions.  The  solution  of  the  quite  complex  system 
of  the  latter  is  reduced  to  the  solution  of  a  system  of  linear 
equations,  which  substantially  reduces  the  volume  of  calculation 
work.  Independent  from  the  complexity  of  the  problem,  the  number 
of  variables,  and  the  nature  of  the  connections  between  them,  a 
solution  can  be  obtained  in  the  form  of  an  explicit  analytical 
dependence,  which  allows  study  in  general  form. 

The  limiting  value  of  parameters  changes,  with  which  the  use 
of  a  method  of  small  deviations  is  possible,  is  determined  by  the 
value  of  permissible  error  of  the  final  result.  The  relative 
changes  of  specific  thrust  and  specific  fuel  consumption,  connected 
with  the  deviation  of  any  parameter  by  10—15%  >  are  almost  exactly 
defined  as  the  product  of  a  relative  change  of  the  given  parameter 
by  the  numerical  value  of  the  coefficient  of  influence  [1],  The 
total  values  6R  and  during  a  change  in  several  parameters 

are  defined  as  the  algebraic  sums  of  the  corresponding  partial 
increases. 
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The  calculation  of  the  numerical  coefficients  of  influence 
in  general  solutions  turns  out  to  be  sufficiently  simple.  In  all 
cases  the  range  of  applicability  of  the  method  of  small  deviations 
is  expanded  with  an  increase  in  the  initial  values  of  tt#,  tt#,  Tf,  etc. 

H  T  j 


§  1.  Equations  of  Operating  Conditions 
of  the  Basic  Elements  of  DTRD  in  Small 
Deviations 


Let  us  compose  the  fundamental  equations  of  the  process  in 
the  elements  of  a  twin-shaft  DTRD  with  separate  exhaust  from  jet 
nozzles  with  complete  gas  expansion,  and  with  air  bleed  due  to 
compressor  on  the  cooling  of  high-pressure  turbine. 


Let  us  assume 

-  *  —  -  a 

*11  *%mj|  • 

We  will  designate  coefficients  of  influence  for  the  entire 
engine  in  terms  of  K,  for  the  low-pressure  [LP]  ( HU)  stage  and 
the  ventilator  -  A,  and  for  the  [VD]  (B/0  stage  -  B  with 
corresponding  indexes. 


The  specific  thrust  of  a  DTRD  is  equal  to 

n  __  4*  y  ^V4ll 

1+j' 

Logarithmizing  and  differentiating  this  expression,  we  obtain 
the  connection  between  the  relative  increases: 


_7*  __  “  ”  _  ,l  Rj*n 

Ryi  Ry.\l  +  y  ^?y.UI  ’  /?yJ,  4*  y  /?VJ„ 


(1) 
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— . — - 


^  -  - . . . 


rr«y»i-v:Ti  if 


w*it 


Let  us  Introduce  designations  AR  /R  *  dR  /R  “  ^  n*  etc 

Jr  n  y  in  jr  n  j  n  y  rM 

Then  equation  (1)  is  written  in  the  following  form: 


where 


?  Ryx  =  Kx  RnX  +  (I  -  AT,) «  RyM  +  (1  -  /Qty, 


(la) 


^nl  “(■  y  ^)T4ll 


*t= 


.y 

1+y 


Specific  thrust  of  contour  I  is: 


C;1— •  =  ’  < v  I  23ior4«v-\'). 


Correspondingly,  after  logarithmic  operation  and  differentiation, 
this  expression  in  small  deviations  will  take  the  form: 


* =  <2*?pc«  +  5 T**  +  *4 « «V,)  -  (Kt-  1 ) IV. 


(2) 


where 


v — C‘\ _ ij _ y  . 


v  °.25  . 

-  e„  «.»  _  |  • 


|Kl 


3rpc,  =?n,  +  Sr.,,*  +  ii„*  -  in,*; 
*r4  =  3.5/C,  ?,W  +  $«;,*. 


(3) 

C^O 


0.4AP 

*v_  !  +  0,2/VJ*’ 


The  temperature  of  gas  after  the  turbines  is 


t*  —  t  *  L™*  —  t  *n  -  s  ** 

«  4  —  ><u - pjg —  '4»A  l*  '»*  tIihj  ;• 


By  simlltude  we  obtain  the  expression 


®  —  /4U  (r4|  -f- 


(5) 
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Tr™.  “"“"i 


In  small  deviations  this  expression  will  be  written: 


where 


«  V  =87-,*  +  Kt. (K9 8r„*  -  5 t,.,*). 


Uui  _  0,286k.,"'** 
/'!l  5  n.,*  ji.,*0"*-!  ’ 


IT I 

j/  _  k|  _  1 

l'|i» —  7*  * 


i+K 


(9) 


where 


5  7,*  =  o7'„  +  =  +  Ky  '‘M, 


K»~  =  when  H  <  11  km; 

'  H  *  H 

K  =  0  when  H  >  11  km. 

H 


(10) 


Specific  thrust  of  contour  II  is 

Ryx  —  j  (Vpdi  I 

The  expression  for  R  in  the  low  deviations  takes  the  form: 

/?y4ii  =  /i,  (2'»‘ppt,| -j- <.  7a||*  + /l4  5  rpc„)  (/4S-  (11) 

Coefficients  and  A^  are  analogous  to  and  : 


'•■pell  -  —  ‘Wv 8Af  +  8i„’+ 

3  T.m:  =  >.  Ti*  +  Ai0' (A,  i  s,„*  —  Th„*); 

5l/  =  dAf+*  5  7.. 


(12) 

(13) 

(1*0 


The  specific  fuel  consumption  of  a  DTRD 


C’*-m'Ryil+yRylll' 
hCyA  =  imt  -5/?  -  Ktly. 


' 


(15) 


§  2.  The  Derivation  of  the  Generalized 
Coefficients  of  Influence  of  the 
Parameters  on  the  Design 
Conditions  for  R  and 
C  of  a  DTRD 

y* 

Let  us  find  dependences  5R  , ,  6R  nTT,  6R  and  5C  „  on  6 tt* ; 

F  y,ql  yflll*  y,q  ya  * 

6y ,  6T^,  6M,  6H  on  6n*,  6a*  and  69pC« 

We  begin  the  solution  of  the  system  of  equations  with 
compilation  of  the  equations  of  the  balances  of  power  on  the 
shafts . 


For  the  low-pressure  stage 


^IHJ  ( I  w,)  ■{•  y  Lin* 

=  + (I  -  ,  +  (l 


(16) 


where 


K>  i  = 


£,,,(>  +*»".) 


For  the  high  pressure  stage  the  equation  of  the  power  balance 


takes  the  form: 


£..us£.(U(l+«,-  m«.>, 


4  = 4  +  [ +  ^-^7, 8  m,; 


(17) 


The  coefficient  during  6mT  is  sufficiently  low,  therefore  we 
disregard  the  second  term. 

After  substituting  in  (16)  the  expressions  for  SLTHfl, 
and  6L  TT  and  in  (17)  for  6L  and  5L  ,  it  is  possible  to 

Mil  HBfl  T  Bfl 

determine  the  required  changes  in  values  tt*b^  and  from  these 

equations,  dependent  on  the  change  in  the  remaining  parameters  of 
the  engines. 


iiwwmM UtoMM  SKf,3S«S Pi 


are 


The  further  obtained  expressions,  6tt*  and  6 tt*  , 

T  Bfl  T H,Q 

substituted  in  (3),  (5)  and  (7)  with  the  replacement  6tt* 

K  Hfl 

and  6ir#  ■ 

H  Bfl  Hi  KHfl 


s"!ii 


After  substitution  of  (10)  in  (9),  (9)  in  (8);  (7),  (8),  and 
(9)  in  (6);  (6)  in  (5);  (3),  (*0,  and  (5)  in  (2)  and  some  trans¬ 
forms,  let  us  find  the  dependence  on  a  change  in  the  initial 

parameters : 

i  =  Kt  8  Vl  - I  AT*(1  -  Klt )  K„  8 y  - 
-\K, U#f„  +  (A.,Ai0 -  B,J +  BtiK, (1  +  *,,)]} * ».«*  + 

+  2 Kt (1  —  A’n) +  j  A*|/f4 —  jj  |  + 

-M 1  +  K„)  lA'.X.od  -  *;)  -  M^]|  8  - 

\  K\A"  ('  -  *«  >  <  1  +  -  KuKn  -  ~ 

"  7  ^  {l  (  I  “  K,)  -  B„K;  |  ( I  +  *„)  -  |J  *  W  + 
f  .7  K »  «  n +  J  K*  w,» + ?..*  +  \  K,Kt  5a„f*+ 

+  j  *n|n  +  *„>*„+  *  t%*  +  j(i -k,)  - 

-  2  K  [(»  +  Ku) Kx%  +  #f4  (--  -  3,5)  -  1  jj  0  Af  + 

+  \  J 1  ~ ■**  [< 1  +  K»)  K»  +  ^Jj  K,t  8  H, 


where 

Kvl  —  Ah  -f-  ; 

Substituting  (12)  and  (13)  in  (11),  we  obtain  the  expression 
for  the  change  in  R  : 


’  /?yJ„  —  j  (^«  "I"  ^  rBii*  ~  g  ^*^1®'  ^  t<»ii*  + 

i-  \  a,  a,  s  +  a,  8  Vll +  J-  *.  8/y + 

+  [l-^+j/r##fv(l+3,544)j5AI. 

For  the  relative  fuel  consumption  per  second 


8  m,  =  *„  8  V  +  ( I  -  /Q  [ff„ «  //  +  Afv  *  .**  + 

■f"  K..Kln  '•  K*!41  Atni  ft  HI  8  T,kM  |. 

With  the  aid  of  the  obtained  dependences  it  is  possible  to 
find  the  relative  changes  in  6R  and  6C  in  terms  of  the  change 

jr  iM  j  M 

of  the  independent  parameters  of  the  process  in  both  contours  and 

of  the  flight  conditions,  after  substituting  the  relationships 

found  above  for  6R  T,  6R  TT  and  6m  into  equations  (la)  and  (15). 

yfll  yAli  t 

For  the  determination  of  the  numerical  values  of  the  coefficient 

of  influence  of  the  parameters,  it  is  first  necessary  to  determine 

the  values  of  the  initial  Influence  coefficients  by  the  values  of 

IT?,  *pc>  T2 ’  T*’  T$’  Ry,  etc.  It  is  further  advantageous 

to  calculate  the  values  of  the  combined  coefficients.  Using  the 

general  expressions  of  coefficients  of  influence  and  the  found 

values  ,  K2,  ,  etc.,  the  numerical  values  of  coefficients  of 

influence  are  found.  It  is  more  convenient  to  calculate  these 

coefficients  first  for  6R  T»  5R  „Tt  and  6m,,  and  then  to 

yfli  yflii  t 

substitute  them  into  the  expressions  6R  and  6C  .  The  relative 

yA  yA 

changes  in  the  specific  engine  characteristics  can  be  determined 
after  multiplying  the  calculated  coefficients  of  influence  of  the 
changed  parameters  times  the  relative  changes  in  these  parameters. 
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S  3.  The  Analysis  of  the  Numerical  Values 
of  the  Coefficients  of  Influence  of 
Selected  Parameters ' on  the  Specific 
Characteristics  of  DTRD 

Let  us  compile  a  numerical  table  cf  the  coefficients  of 
influence  for  DTRD  with  initial  values  of  the  parameters: 


y  =  8;  7*,*  =  M00°K;  ^,*  =  35; 

Ti«m  —  —  0,86; 

9. 

II 

o 

00 

J-Nj 

=  W  =  0,92; 

=  0,96:  cKt 

=  0,99;  =  1,0; 

fp«i  ‘  ?pcii  “  0.99; 

mot  =  0,08. 

Rated  flight  conditions: 

M  =  0,8;  //  »  ll  km. 

Gas-dynamic  and  thermal  design  are  determined  by  the  following 
values  of  the  parameters  of  the  process: 


rUI*  =*  278.5°K; 

Wpti  =*  2,285; 

r„*  -  751. 9°K; 

/.[Mi  =  32400  kgm/kg; 

Lull  -  .15-10  kgm/km;  Tiu*  --=  945,2°K; 

rt*  67I,0:K; 

-  5,27; 

"i«i  ~  4.5.3: 

:i 

u 

w 

in 

mT  --  0.0174; 

-  28.2  kg- s/kg: 

=  10.45  kg- s/kg; 
Cyjt  ■-=  0,562  kg/kg -h. 

Ryx-  12,41  kg- s/kg; 

Further  according  to  formulas  we  find  the  numerical  values 
of  the  initial  coefficients: 


J,  5.«  'iyw> 


A,  -  0.227 

A „  -  -  0,33; 

Ai  -  0,252. 

A,  =  0.888 

A,  =  1,858; 

Kt  -  1,19; 

At  =  0,207 

Kj  =  0,935; 

A*  -  2,16; 

A9  =>  0,448 

A10  =  0,676; 

A|,  -  0,1105; 

A„  =  2.588 

Au  =  1.345; 

■h  =  3,31; 

At  =  1,088; 

As  =  0,546; 

'h  =  0,405; 

•1,  -  2,61; 

A,  0=  0,125; 

B,  -  0,472; 

Bs  =  0,438; 

«»  =  0.482; 

B, 0  =  0.629; 

A#A!0  =*  0,303; 

M,0™  0,326; 

B,B. 0  =  0,303; 

Ms  .  0.221; 

BaBf  =  0,207. 

Substituting  these  values  in  the  general  expressions  of  the 
coefficients  of  influence,  let  us  compile  a  numerical  table  of 
coefficients . 

From  the  obtained  table  it  is  possible  to  directly  define 
a  series  of  the  features  of  DTRD  with  the  given  values  of  the 
initial  parameters. 


So,  according  to  the  value  of  the  coefficients  of  influence 

the  advisability  of  increase  T§  for  Increase  in  R  is  apparent, 

i  y  A 

but  with  some  increase  in  An  increase  of  ir*j  and  y  somewhat 

decreases  C  and  R  .  The  changes  in  a*  and  <p  TT  exert  a  strong 
yA  yA  bx  TpcII  c 

influence  on  R  and  C 

yfl  yfl 

With  the  aid  of  the  table  it  is  possible  to  estimate  the  values 

of  specific  parameters  of  an  engine  with  other  initial  data.  Let 

us  suppose  it  is  necessary  to  calculate  R  and  C  of  an  DTRD  with 

y  a  y  a 

the  parameters:  T^  =  150Q°K;  tt =  40;  y  =  9  during  flight  at 
H  =  10  km  with  M  =  0.9. 

H 
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mm 


The  values  of  given  magnitudes  differ  from  the  corresponding 

Initial  values  by  less  than  15$.  After  multiplying  relative  changes 

in  these  parameters  to  the  coefficients  for  them,  and  after  totaling 

the  results,  it  is  possible  to  determine  that  R  decreases  by  9.56$ 

y  A 

and  increases  by  1.62$,  which  v'ill  comprise  correspondingly: 

R  ■  11.23  kg  s/kg  and  C  =  0.571  kg/kg  h. 

j  j  iM 


Table.  Numerical  values  of  the  coefficients 

of  influence  of  the  selected  parameters 

during  design  on  R  and  C  of  DTRD. 

yA  y  A 


»/?ya 

ic,k 

c,  n.,* 

—0,612 

0 

-0.154 

-0,352 

-0,198 

;  _  * 

*  "«ll 

-5.16 

2,34 

0.45 

0 

-0,45 

b' 

-2,14 

0 

-0.539 

0 

—0,349 

ira 

6,86 

0 

1.74 

2,16 

0,42 

*  7ixU  ' 

2,14 

— 0,<07 

0,384 

0 

— 02184 

4  T<X../ 

0.703 

0 

0.177 

0.125 

-0.052 

*  w 

3.88 

0 

0.978 

0,73 

-0,248 

2.02a 

0 

0.51 

0 

-031 

'■ 

2.343 

0 

0,592 

0 

-0,592 

?  3  * 

s 

1,106 

1.8 

1,623 

0 

-1,623 

1.106 

» 

0,278 

0 

-0.278 

1.838 

0 

0,468 

0 

—0.468 

> 

'•  ~  Wtl 

0 

3.31 

2.47 

0 

-2.47 

i.  if 

-1,326 

-0,605 

-  0,713 

-0.2G4 

0.4 ‘9 

IH 

2.1 

— 0,165 

0,406 

0.383 

-0,023 

The  method  of  the  low  deviations  makes  it  possible,  with 
considerably  less  expenditures  of  time,  to  compare  a  multitude  of 
diverse  variants  and  in  the  narrow  area  of  the  combinations  of 
parameters  to  find  their  optimum  relationship  to  provide  for  the 
required  specific  fuel  consumption  and  an  acceptable  value  of 
specific  thrust. 
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METHOD  OF  CALCULATION  OF  THE  MIXING  OF  GAS 
FLOWS  IN  THE  DIFFUSERS  OF  THE  AFTERBURNER 
MIXING  CHAMBERS  (AMC)  OF  DTRD 


Ye.  V.  Barabanov 


Accepted  Designations 


c  =  c  ,/c  , 

pnl  pal 

-  VFi 
R3  = 


6  =  tg  (0np/2) 


-  longitudinal  component  of  flow  rate; 

-  transverse  component  of  flow  rate; 

-  the  average  (by  area)  flow  rate  in  a  given 
cross  section; 

-  specific  heat  ratio  of  passive  gas  to  active; 

-  expansion  ratio  of  diffuser; 

-  equivalent  radius  of  annular  diffuser; 

-  given  angle  of  flare  of  annular  diffuser; 

-  expansion  coefficient  of  rectilinear  diffuser. 


Indices 


parameters  of  active  flow  at  inlet  to  the  diffuser  of 
an  afterburner  mixing  chamber  [AMC]  ($CH); 
parameters  of  passive  flow  at  inlet  to  diffuser; 
total  (averaged)  parameters  of  the  inlet  to  diffuser; 
parameters  of  mixed  flow  at  output  from  diffuser; 
parameters  on  the  x  axis. 


The  present  work  examines  the  overall  task  of  calculating  the 
mixing  of  gases,  differing  in  chemical  composition,  '(c  and  x)  in 
the  annular  channel  of  a  variable  cross  section,  which  has  a  direct 
relationship  to  the  design  and  study  of  the  mixing  chambers  of 
[DTRDF]  (^TP^tti) ,  and  also  the  different  compound  propulsion  systems 
( [TJE]  (Tfm),  [TJEW  AB]  (TPfl*)  +  [RJE]  (riBPA)  »  DTRDF  f  RJE ,  etc.) 
in  the  diffuser  of  their  area  behind  the  turbine. 

Unfortunately,  the  predominant  majority  of  known  methods  of 
calculation  [6]  is  related  to  cylindrical  mixing  chambers  and  cannot 
be  completely  used  by  us  in  the  present  work.  Thus  this  article 
is  dedicated  to  the  construction  of  the  analytical  method  of 
calculation  of  the  mixing  of  heterogeneous  flows  in  the  diffuser 
with  the  purpose  of  calculation  of  J  p  dP,  which,  unlike  cylindrical 
channels,  in  a  diffuser  is  not  equal  to  zero  and  introduces  a 
substantial  change  in  the  entire  procedure  of  calculation. 

§  1,  Basic  Assumptions 

1.  At  the  inlet  to  the  diffuser  of  AMC  the  parameters  of  each 
of  the  flows  are  distributed  on  their  cross  sections  evenly. 

2.  On  the  circuit  of  the  diffuser  of  AMC  the  nonuniformity 
of  the  dimensionless  speed  is  described  by  Schlichting' s  law  [3]: 


where 


T  u  -  !Lr  ?  _  r  (-:)  -  A  (1 ) 
ii,  —  w  1  -  A  (T)  * 


f  (;)  =  (! 


(1) 

(2) 


3.  Static  pressure  in  every  cross  section  is  assumed  constant. 
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^ .  We  consider  the  mixed  gases  incompressible.  The  latter 
two  assumptions  are  introduced  on  the  ground,  that  in  the  diffuser 
of  AMC  essentially  the  subsonic  flows  are  examined,  since  to 
provide  for  a  stable  combustion  of  mixture  in  afterburner  it  is 
necessary  to  have  X^  <  0. 2-0.3*  and  M^  =  0.55-0.7  [^]. 

5.  The  hydraulic  losses  of  gas  friction  against  the  walls  of 
the  diffuser  we  disregard  in  view  of  its  low  extent,  which  is 
dictated  by  the  very  idea  of  AMC  as  a  means  of  economy  of  the 
overall  sizes  of  the  output  part  of  a  DTRDF. 

6.  The  processes  of  flow  and  mixing  of  streams  in  the  diffuser 
of  a  AMC  are  stationary. 

§  2.  The  Fundamental  Equations  of  Mixing 

As  it  follows  fron  assumptions  3  and  4  adopted  above,  at  the 

inlet  to  the  diffuser  of  AMC  pA,  »  p  , ,  or: 

A1  nl 

/.a,)  =—?,,/,,l)  •  (3) 

Under  such  a  condition  of  Inflow  the  fundamental  equations  of  the 
mixing  of  heterogeneous  flows  in  the  diffuser  of  AMC  are: 

1.  The  equation  of  continuity 


—  CA|  ( I  -f  .y). 


2.  The  equation  of  energy 


fri4',,*  r»*  —  7\i*  -f-  (7n,  Cp„  V, 


whence  is  determined 


,>  *_  _  l  -{-yc  w* 

'*  ~''M  7ai*“  14  7" 


and  stagnation  temperature 


r  *  _  /•  *  i  y e 

•»  —  1  A»  —  ,  —  • 

i  +y<- 


(5) 

(6) 
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3.  The  equation  of  momentum 


which  after  using  a  Kiselev  [1]  formula  of  impulse,  and  also 
formulas  (4)  and  (5),  is  converted  to  the  form: 

*  <'*>  1  (l+jO(T+j>f  »*)  =  BAtz  (XAI)  + 

A» 

+  BHtzOl„)ylrc'y  +  Cj  pdF, 

where  1 


The  last  member  of  (7),  which  differs  it  from  similar  momentum 
equations  for  cylindrical  chambers,  we  convert  to  the  more 
convenient  form: 


=  |  2(Xr  — »)  fc 


mKpr  *.V(^Al) 


z. 


jp,“l 


(8) 


Considering  that  the  constant  before  the  sign  of  the  integral  of 
(8)  is  a  function  of  the  parameters  of  the  1st  contour  of  a  DTRD 
(i.e.,  combustion  products),  and  by  substituting  (8)  in  (7),  we 

will  obtain  in  final  form  the  equation  of  momentum  for  the  diffuser 
of  AMC: 


l  (1  +  !>*)=;  + 


Bnl2(Ki)y^rr 


1,045(1  +  q) 
.V(*Ai)a 


(9) 
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It  is  completely  obvious  that  the  system  given  above  of 
equations  of  mixing  can  have  a  determinate  unique  solution  only 
in  the  search  of  the  integral  of  static  pressure  in  the  last  equation. 


3.  The  Calculation  of  the  Integral 
of  Static  Pressure 


For  the  calculation  of  the  integral  of  static  pressure  in 
general  form  (for  isolated  special  cases  successful  solutions  [5,  11] 
are  known)  it  is  necessary  to  know  the  analytical  dependence  of 
static  pressure  distribution  along  the  diffuser,  of  AMC.  For  the 
quasi-stationary  turbulent  flow  of  a  Newtonian  fluid  such  a 
dependence  can  be  derived  on  the  basis  of  the  Navier-Stokes 
differential  equations  of  motion  and  continuity  which  for  the  case 
interesting  to  us  are  converted,  as  is  known  from  [8],  to  the  form: 


dp 

dx 


Id,  .  du  ,  du 


where  tension  of  eddy  viscosity,  according  to  Prandtl  [2]: 


(10) 

(ID 


„du  du 
'■'  =  *r  Ty  Ty 


(12) 


The  joint  transforms  of  equations  (10)  and  (11)  reduce  to 
the  form: 

y‘*dp  _ d  .  d 


yap  up  dp 

j  f “ydy' 


Utilizing  in  (13)  an  equation  of  the  constancy  of  consumption 


f>uydy  =  0 


(13) 


and  having  satisfied  the  boundary  condition  on  wall  y  ■  Rg ,  tt  *  0, 
we  obtain  the  fundamental  integrodifferential  equation  for  static 
pressure  distribution: 


dp  _ 

d  X 


2n  J 
F  dx 


u*ydy. 


(14) 


which  can  be  solved  only  in  conjunction  with  equation  (13). 


For  convenience  and  the  universality  of  the  solution,  we  will 
introduce  the  dimensionless  coordinates  and  coefficients: 


where 


x  =x 


a  a 


K  P  U 
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(15) 


F  _  V:J  z.  +  i . 

I  -  |  * 


V 
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We  solve  equation  (13)  relative  to  the  complex  coefficient  of 
diversity  X,  after  making  the  replacement  of  variables  in 
correspondence  with  (15)  in  it  for  this,  utilizing  Schlichting* s 
law  (1)  and  a  procedure  of  transforms  from  work  [10],  whereupon 
we  obtain: 


where 


,  - --  -</}.  _  d\'  f  dV  f 

>  /(I  -t-o >•)--  =  -  *(>.  2  >.(<+>. d)  l-L  +  2e~\A 

dx  dx  dx 

— t§=  ‘  »-rf>  *«>“»;= 


(16) 


/<;> 

/(f) 


~  /.(;) 
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l  -  mixing  length  where,  in  turn: 

^  (€)  =  /» <0M  V  f(5)l  +  ^(D/Mf)  *  2|fl(l){  —  2B(0|; 

.V({)=18{*?(*>; 

/(? ) «  2/i  (;)  1 1  +  A  ( 1 1 1  +  2 1 P  0  ( 1)  -  B  (?)); 

H(()  =  /}(f)-;*B0)-2p(<).  ln  which> 

£ 

/»(;)  =  >! (5) I);  .i(i>  =  2 

o 

;  *  • 

^ ( 1 )  =  2  jnDldl- 0,2571;  fl($>  =  2  / 

o  o 

1 

/i(l)  =  2  y  ?*(;)?</?  =  0,1334, 
o 

and  <p(£)  is  in  accordance  with  (2). 


) 

i 


In  the  case  interesting  to  us  of  the  mixing  of  flows  in  a 
diffuser  with  rectilinear  generatrices  VF  =  1  +  8x,  the  obtained 
Abel  equation  (16)  of  the  2nd  order  [9]  is  expressed  ln  quadratures 
and  is  reduced  to  the  following  differential  equation  with  separable 
variables . 

_ (1  -f  <7 \)  d*. _ _  dx 

(I?d  --  2(pr-f  *)X  +  (2?r-’*)~f+5i’  (l8) 

The  relative  coefficients  entering  this  equation  were 
calculated  by  computer  as  the  mathematical  expectations  of  the 
random  variable  £  by  Simpson's  method  and  the  following  mean- 
integral  values  obtained: 

ac|»  =  0.9.3;  ctp  =  2.714:  rfc„  =  0.573:  i*fp=~  3,287;  *tp  =  39,6. 
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With  consideration  of  empirical  data  <(£)  of  different  authors 
[7,  8] 


Equation  ( 1 8 )  was  solved  by  the  method  of  the  undetermined 
coefficients  of  the  initial  conditions  x  *  0,  \  =  T-j.  : 


where 


-'•I  -  Yi 


M 

=  I  +  '?*> 


_  M1  f  «Y|) 

"  ~  (2?  d  -  k)[ Y,  -7j) 


m  = 


Ml  ±  qyJ 

(2?d —  *)  (7,  —  7,)’ 


(19) 


while  ,  y2  are  the  roots  of  the  trinomial  of  the  denominator  in 
equation  (18). 


After  the  transforms  of  equation  (19),  we  obtain  the  final 


solution  in  the 


where 


form : 


a  =  _L±?_l 

9-2(\+a) 


1 


2?  (c+ d)  +  k 

k~2'?d 


(20) 


From  this  equation  it  follows  that  A  depends  on  the  geometry  of 
the  diffuser  of  AMC  (B,  f)  and  the  gas-dynamic  coefficient  of 
diversity  at  the  inlet  to  it 

■f-_  *A, 

'cpll'  KEm 

since,  according  to  assumption  1  (see  S  1) 


=  V  =  V;  E,  =  £m. 


Now  already  according  to  equation  (14)  we  can  find  the 
distribution  law  of  static  pressure,  for  which  we  will  carry  out 

p 

both  of  its  parts  to  a  constant  value  of  ram  pressure  Pwcpj/2  In 

the  reference  section  and,  remembering  that  w  /w  T  =  F,/F  =  1/f. 

op  cpi  1 

after  the  transforms  with  consideration  of  (1),  (15)  and  (17),  we 
obtain 


58 


'IMW  MIIUII  lllimillWIlWll.MPIl  III. 


WJ-'P  BWBPW^W.  W(Plji'i  WIP  'WF"W  JHWW 
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•>*  'rwywv 


*■7  ■*,- 


d(--£_\  ..  ,  .  _  , 

\ pwfn)  _dp  _  s  d  »  /■•  —  2h  >■  4-  / 

dx  dx  /  dx  I  7"  J 


(21) 


where  in  accordance  with  (17):  h  =  B(l)  -  A^(l)  =  0.067^; 


1  +  B(  1)  -  2A(  1)  =  0.6192;  s  =  2/Cl  -  Ad)]*  =  3.62. 


We  Integrate  equation  (21),  having  preliminary  replaced  its 
right  side  with  the  differential  of  the  polynomial  to  the  same 
degree,  but  with  different  coefficients,  and  f  -  by  formula  (20), 
whereupon,  after  conducting  the  transforms,  we  obtain  the  resultant 
expression  for  static  pressure  distribution: 


where 


p—pt-ip  =  (<t>V-  -j.  ).,  +  «•»)_ 

<t>  -  ?^i4- 1 ?(a-  2d)  +  *);  -i  =  0,488 |l  +  ; 


(22) 


=  1.81  Jo, 


•  =  i.ai  10,6192  -  0.0674  /=-.2'i(a  —  d)  +  k. 


As  we  see,  formula  (22)  is  given  as  a  function  of  X  and  static 
pressure  distribution  along  the  diffuser  of  AMC  can  be  constructed 
only  during  the  joint  solution  of  equations  (22)  and  (20).  Since 
the  latter  is  a  complex  dependence  of  X  on  6,  f,  X^,  which  cannot 
be  expressed  in  an  explicit  form,  then,  for  the  imparting  of  an 
engineering  nature  to  the  solution  of  the  system  of  these  equations 
by  the  computer  of  "MIR-1,"  equation  (20)  for  the  entire  reasonable 
range  of  parameters  of  the  diffusers  of  AMC  DTRDP  was  computed: 


0_  _/2  =  3-16°  with  pitch  of  1°  and  A  =  l-\  =  y  with  pitch 

np  Max  ' 

AA  =  0.05.  As  a  result  of  these  calculations  a  universal  nomogram 


Fig.  1,  was  constructed  for  any  value  of  A^,  with  which  it  is 
possible  to  find  both  =  (X^  -  l)1^^/(y  -  X^)q  and  X  in  any 


(dependent  on  f)  cross  section  of  the  diffuser  (see  Fig.  1)  for  a 


given  ,  which  significantly  simplifies  practical  calculations. 
After  thus  finding  the  value  of  A  according  to  formula  (22)  from 


Fig.  1,  3t  is  now  unambiguously  possible  to  determine  Ap  in  this 
cross  section. 
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Fig.  1. 


However  with  the  purpose  of  achievement  of  greater  accuracy 
and,  simultaneously,  of  engineering  calculations  we  have  conducted 
the  Joint  solution  of  equations  (20)  and  (22)  for  the  system 
Ap  =)  _  _ 

Y  _>  <p(A)  with  AX  =  0.05  over  a  wide  range  of  the  gas-dynamic  and 

geometric  parameters  of  the  subsonic  diffusers  of  AMC  on  the  "MIR-1" 

computer:  A^  =  1. 0-1.8  and  0  =  10-30°.  According  to  the  results 

of  this  calculation  the  graph,  Fig.  2,  of  Ap  =  qp(f)  was  constructed, 

which  is  also  universal,  since  in  the  indicated  range  of  0  °  it  is 

n  p 

applicable  for  diffusers  of  different  length.  The  presence  of  the 
graph  Fig.  2  makes  it  possible  in  engineering  calculations, 
generally,  to  not  use  formulas  (20)  and  (22),  but  to  directly  take 
the  necessary  value  of  Ap  from  the  graph. 


h-js  . 


JirffT  VwV 


! 

-  i - 


gw  9-KbQ*  -\ 
-g*  9-JO0 
-9*  9-frSO*  f 
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Knowing  the  dependence  of  static  pressure  distribution  on  f , 
we  can  now  solve  the  basic  task  of  calculation  of  the  integral  of 
static  pressure.  As  it  appears  from  equation  (9),  this  integral 

f. 


ought  to  be  solved  in  the  form 


f  E" 

J  Pi 


df.  In  accordance  with  (21)  and 


(22)  let  us  introduce  Into  the  integrand  the  relative  function 
P  =  p/pw^/2: 


where 


u  /, 
/%«-■ 


1 1P*/. 


H  =  ?— *  = _ & _ i»2L _ 

‘  *Pi  *  Z.+  1/,  *■-!. 

\  ~  +  I  A,tp  / 


Let  us  make  "integration  by  parts"  of  the  obtained  integral,  after 
replacing  for  this  the  differential  dp  with  its  expression  (21): 


j>ipdf=z\ip/-n  Jfdp=z\t  p(X)/CX)  _+ 
+ 1 +  '  j = /7,^]  + 
+,  J**£=**&+i  _  (/,  v-  »*;+ q|. 

“*  J  A  J 


Considering  that  u  =  1/p^ ,  f^  =  1,  and  in  accordance  with  (22), 
=  Pj  +  Ap,  the  final  solution  of  the  integral  of  static  pressure 
for  the  diffusers  of  AMC  with  rectilinear  generatrices  obtained 
in  the  following  final  form: 


/  =  </,-  1 )  4-  !‘ A pf \  +  3,62,J '1 -***■+*  _ 

^  /  i 


-  (A  >.  +  *)]>  (23) 

where  f^  -  the  usually  assigned  magnitude  of  the  expansion  ratio 
of  a  diffuser;  Ap  -  the  increase  of  static  pressure  at  the. assigned 
diffuser  length,  computed  according  to  formula  (22)  or  determined 
from  Fig.  2;  X^  -  determined  from  Fig.  1  from  assigned  magnitudes 
of  0Q >  f  and  from  the  computed  value 


where 


i  _  '*A1 

'tpi  1  'V  £.\1 


f  —  Am  ~t~ )  A  ~  l 
'•“/M-I/Xfl 


According  to  formula  (23)  we  conducted  the  calculations  of 
f 

J up  df  for  that  diffuser  length  of  the  space  behind  the  turbine 

1 

which  is  most  acceptable  for  aircraft  engine  construction,  T  =  2.5 

calibers  [4],  in  the  following  range  of  parameters:  0  =  12-28°; 

np  ’ 

\1  =  1-1.8;  and  Mjcp  =  0. 2-1.0.  The  results  of  the  calculations 
are  given  ..  r.  Fig.  3,  from  which  it  is  clearly  evident  that  to  the 
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§  4.  Determination  of  tt  ,  in  the  Diffuser 
of  AMC  3 


The  basic  purpose  of  the  calculation  of  any  mixing  chamber  is 
the  determination  of  its  criterial  parameter  -  the  pressure  ratio 
=  p3^pnl*  always  possible  to  find  the  total  pressure  of 
mixture  at  the  chamber  outlet,  which  is  especially  important  in 
our  case  of  calculating  the  diffuser  or  AMC.  Thus  for  its 
determination  we  utilize  in  relationship  (4)  the  known  consumption 
equation : 

/V  ^*0  (*j.  ^3)^:1  -  Pm*  Dxtq0.xl.*M)FM  (  ^**’0  +!')• 


Utilizing  in  this  the  equality  of  expression  (5)  and  F,  =  f  (F.,  + 

5  fl  A1 

+  F  , ) ,  after  the  transforms  we  obtain: 
nl 


__  ^.M  0  *A\|)  8 

«.*(>  i.'^AO+a) 


(1  +.v)(l  rye  *►*). 


(24) 


where 


(25) 

(26) 


Since  we  are  investigating  the  mixing  of  heterogeneous  flows, 
the  gas-dynamic  functions  in  the  given  formulas  are  each  found  for 
their  own  index  x.  The  adiabatic  index  of  the  gas  mixtures  is 
determined  from  the  known  formula: 


^:t  —  \\  1 


1  +yc 

1  "I"  y  C  M'^nl 


(27) 


The  integral  of  static  pressure  affects  the  value  tt3  in  terms 
of  a  function  of  q(X^,  x^)  in  (24),  determined  in  accordance  with 
equation  (9). 


In  the  method  given  above  the  calculation  of  ir  was  conducted 

depending  on  the  available  pressure  differential  nQ  ■  P*]/P*i  *>or 

a  series  of  diffusers  of  AMC  in  the  range  0  *  12-28°;  and  also 

for  the  cylindrical  mixing  chamber  (0  «  0°),  for  which  equations 

(9)  and  (24)  were  simplified  due  to  the  introduction  of  conditions 

f 

of  cylindricity :  J p  d7  *  0  and  ffl  »  1. 

1 

The  calculation  was  performed  at  constant  values  of  *  0.45 
and  Mlcp  =  0.6,  those  most  characteristic  and  optim  .m  for  space 
behind  the  turbine  of  a  LTJE]  (TPA)  and  DTRD  [4,  5]»  As  can  be 
seen  from  graph  4,  constructed  according  to  the  results  of  the 
aforementioned  calculation,  the  values  in  the  range  actually 
used  in  DTRD,  no  =  1.0-1. 3  for  a  cylindrical  mixing  chamber  do  not 
considerably  exceed  that  for  the  diffusers  of  AMC:  for  example 
in  a  diffuser  with  0°  =  12°,  tt3  is  only  0.5-1. 3%  less  than  in  the 
cylindrical  chamber.  However,  as  the  same  calculations  showed, 
in  the  cylindrical  chamber  X ^  at  the  indicated  values  of  iro  and  tt3 
it  remains  at  a  level  =0.6;  while  in  diffuser  with  0  =  12°  at  the 
same  values  of  vQ  and  tt3  ,  X^  =  0.24,  i.e.,  it  completely 
corresponds  to  the  conditions  of  an  afterburner  [4],  It  we 
consider  moreover  the  savings  of  weight  and  overall  sizes  of  such 
a  DTRDF,  then  the  advantages  of  organization  of  mixing  in  the 
diffuser  of  AMC  are  indisputable  as  compared  with  mixing  ir  the 
cylindrical  chamber  of  DTRDF. 

Conclusions 

1.  The  experimentally  checked  and  sufficiently  substantiated 
initial  assumptions  form  the  basis  of  the  proposed  method  of 
calculation  of  mixing  in  the  diffusers  of  afterburner  mixing 
chambers . 
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Pig.  4. 


2.  On  the  basis  of  continuity  and  the  joint  solution  of  the 
differential  Navier-Stokes  equations,  the  analytical  calculation 
of  the  integral  of  static  pressure  for  diffusers  of  AMC  with 
rectilinear  generatrices  is  given  in  the  final  form,  which  afforded 
the  possibility  to  transform  the  equation  of  momentum  in  connection 
with  the  diagram  in  question. 

3.  On  the  basis  of  the  proposed  dependences  universal 
nomograms  for  the  single-valued  determining  of  the  criterion  of 
diversity  of  T,  Ap,  and  J"p  df  are  calculated  as  a  function  of  the 
Initial  gas-dynamic  and  geometric  parameters  of  the  diffuser  of 
AMC. 

4.  The  calculations,  carried  out  according  to  this  method, 
showed  the  great  advisability  of  organization  of  the  mixing  of 
flows  of  both  contours  of  DTRDP  directly  in  the  diffuser  of  the 
AMC,  rather  than  in  the  preliminary  cylindrical  mixing  chamber. 
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THE  ExpERIMENTAL  STUDY  OF  THE  MIXING 
OF  GAS  FLOWS  IN  THE  DIFFUSER  OF  THE 
AFTERBURNER  MIXING  CHAMBER  OF  A 
DTRDF 


Ye.  V.  Barabanov 


The  present  article  Is  dedicated  to  the 
experimental  research  of  the  mixing  of  con¬ 
centric  heterogeneous  gas  flows  In  the 
annular  diffuser  of  the  afterburner  mixing; 
chamber  of  the  DTRDF  of  a  supersonic  aircraft. 

Investigation  is  conducted  over  a  wide 
range  of  a  change  of  the  gas-dynamic  and 
geometric  parameters  of  diffuser  nixing 
chambers .  The  optimum  parameters  of  such 
chambers  are  revealed  and  the  evaluation  of 
their  effect  on  relative  specific  impulse  is 
produced,  as  is  the  efficiency  of  the  mixing 
process  on  the  most  "integral"  criterion. 

The  concrete  quantitative  and  qualitative 
derivations  and  the  recommendations  regarding 
the  design  of  optimum  diffusers  of  afterburner 
mixing  chambers  are  given,  as  is  the  selection 
of  the  most  logical  parameters  of  mixing  in 
them.  ° 
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Accepted  Designations 


p*.  »,-T» 


P 

G 

F 


n  =  F  /F, 
A  A  1 


a. 


np 

».  ,  T-./T*, 


*  p!i/p 


«  _ 


y  =  Gnl/GA1  " 


IT  =  pS/p*T 

3  F9  Fnl 
X  =  w/a 

Hp 

J  . 
y,q  i 


J  =  J  n/J  , 
yn9  yfll 


p*  -  P^/Pjp! 


pressure  and  temperature  of  stagnation 

static  pressure  in  the  contour  of  device 

flow  rate  per  second  of  air  (gas),  kg/s 

cross-sectional  area 

expansion  ra'io  of  the  diffuser 

given  angle  of  flare  of  the  diffuser 

available  drop  in  total  temperatures 

available  drop  in  total  pressures 

bypass  ratio 

the  pressure  ratio 

velocity  coefficient 

specific  impulse  in  the  appropriate 

cross  section 

relative  specific  impulse  at  outlet  from 
the  diffuser 

total  pressure  recovery  coefficient 


t*„  =  p* 


Max  ^cp9’ 


p9 

r#  =  T*  /T* 
t9  Max'  cp9 


ft  - 


ht/cm 

cen 


variation  factor  of  pressure  and 
temperature  in  the  exhaust  section  of  the 
mixed  flow 

relative  variation  factor  of  flow 
2 


kg/cm 

s 


Indices 


A,  - 


n,  - 


1,  2,  3,  ...9  - 


9  - 


parameters  of  primary  (active)  flow  at 
inlet  to  diffuser 

parameters  of  secondary  (passive)  flow  at 
inlet  to  diffuser 

the  total  parameters  of  inlet  to  the  model 
cross  sections  along  the  length  of  the 
diffuser 

parameters  of  mixed  flow  of  outlet  from 
diffuser 
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Aataaliid  BttafcfcaL&ifc 


. 


. 


Contractions 


[DTRDF]  (ATPAG)  - 
[AMC]  (OCH)  - 


turbofan  engine  with  afterburner 
afterburner  mixing  chamber 


Introduction 


In  connection  with  wide  application  of  forced  DTRD  on  super¬ 
sonic  aircraft,  the  selection  of  the  most  logical  layout  of  a  DTRDF 
with  mixing  chambers  becomes  urgent.  The  advisability  of  the 
setting  of  such  a  task  is  determined  by  the  trend  toward  reduction 
in  weight  and  overall  sizes  of  the  entire  power  plant  of  a  super¬ 
sonic  aircraft,  without  worsening  considerably  its  gas  dynamics. 

The  existing  requirements  in  gas  dynamics  [1-3]  to  provide 
satisfactory  mixing  recommend  construction  of  cylindrical  mixing 
chambers  with  a  length  of  6-7  calibers ,  and  afterburners  with  a 
length  of  2-3  calibers.  Thus,  for  a  satisfactory  afterburning  of 
the  mixed  flows  in  DTRDF  with  a  mixing  chamber  a  very  bulky  and 
heavy  system  behind  the  turbine  will  be  required. 


At  the  same  time  the  numerous  experimental  works  on  mixing 
in  the  cylindrical  chambers,  which  were  carried  out  earlier  [^-6], 
did  not  confirm  the  assumptions  of  researchers  for  the  expected 
thrust  increment  of  DTRD  as  a  result  of  the  increase  of  completeness 
of  mixing.  As  a  result ,  there  arose  the  need  for  study  of  the 
mixing  of  gas  flows  directly  in  the  diffuser  of  an  AMC  of  a  DTRDF 
with  the  goal  of  development  of  the  optimum  mode  and  geometry 
of  such  a  diffuser  located  in  the  space  behind  the  turbine. 


The  study  of  the  mixing  of  heterogeneous  flows  in  the  diffuser 
of  an  AMC  was  conducted  for  the  solution  of  this  problem  on  a 
special  experimental  device  [7]  which  simulates  the  operation  of 
the  exhaust  systems  of  DTRDF  (Fig.  1).  In  every  testing  were  taken 
the  diffuser  performances  of  determined  geometry  assigned  by  the 
setting  of  female  cone  07  necessary  for  this  (Fig.  2)  with 
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determined  X^  and  by  fixed  settings  of  the  corresponding  exit 
nozzle  05.  In  all,  each  of  the  three  investigated  diffusers  (with 
anp  *  12°,  1 6°  and  20°)  was  tested  in  the  range  of  the  velocities 
corresponding  X^  =  0.2-0. 9* 


Fig.  1.  The  diagram  of  experimental  devices 
for  the  study  of  the  mixing  of  gas  flows. 


Fig.  2.  Experimental  afterburner  mixing  chamber. 
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The  parameters  of  spatial  flow  during  treatment  which  were 
measured  in  the  process  of  the  experiment  at  several  points  of 
every  cross  section  were  averaged  both  along  the  radius  and  along 
the  circumference  of  an  annular  axisymmetric  diffuser.  For  our 
mixer  circuit  of  flows  with  different  p#  and  T*,  in  accordance  with 
[8]  the  principle  of  the  simultaneous  averaging  of  the  parameters 
according  to  impulse  and  according  to  consumption  was  selected. 

§  1.  The  Fields  of  Flow  Parameters 
at  the  Outlet  from  a  Diffuser 

Figures  3  and  4  depict  the  fields  of  the  total  pressures  and 
temperatures  along  an  outlet  radius  of  a  diffuser  measured  in  the 
experiment  and  averaged  in  the  zones  of  the  circumference  of  the 
corresponding  radii.  The  rearrangement  of  flow  pattern  and  the 
transformation  of  the  fields  T#  =  f(R^)  during  a  change  of  ttq  and 
the  transformation  of  fields  p#  =  f(R^)  during  a  change  in  the 
drop  in  temperatures  is  distinctly  visible.  Moreover  as  is  evident 
from  Fig.  3,  a  change  of  t7q  affects,  to  consider  .bly  greater  degree 
than  a  change  of  $#,  the  structure  of  the  field  T*  =  f(Rg). 

However,  in  spite  of  the  great  uniformity  of  the  field 
T#  =  f(R^)  with  a  lew  ttq,  these  projections  differ  by  the  greater 
nondiffusion  of  the  flow  core. 

As  the  consequence  of  greater  uniformity  of  T#  on  the  cross 
section  with  ttq  =  1.0-1.25  and  =  0.64,  the  warming  up  of  mixed 
flow  on  the  periphery  of  the  diffuser  proves  to  be  substantially 
higher  than  at  their  greater  values  of  tto  =  1.53. 

Together  with  this  it  follows  to  show  that,  concerning  the 
increase  in  the  available  drop  in  the  temperatures,  the  mixing  in 
the  diffuser  of  flows  with  small  loses  its  advantage  in  the 
warming  up  of  the  peripheral  zones,  which  is  distinctly  evident  in 
Fig.  3,  when  4#  =  0.4,  where  T#  on  the  periphery  both  for  the  flows 
with  ttq  =  1.25  and  with  ttq  =  1.52,  is  identical  and  equal  to  4l0°C, 
which  is  285°  (or  230$)  greater  than  T#^ . 
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Fig.  3.  Temperature  field  at  the  outlet  of  mixed  flow  from  the 
diffuser  of  an  AMC  with  a  =  16°  ,  n  =  3.2  at  X„,  =  0.7-0. 8. 


Fig.  4.  Field  p#  of  a  mixed  flow  at  the  outlet  from  the  diffuser 
of  an  AMC  with  a  =  20°,  n  =  4.0  with  X.,  =  0.7-0. 8. 


However  fields  T#  =  f(R)  at  low  values  of  ttq  also  in  the  given 

case  turn  out  to  be  more  uniform  than  when  it  =  1.52,  but  this  '  - 

0 

achieved  only  due  to  the  drop  of  in  the  flow  core.  This 
explained,  apparently,  by  the  fact  that  during  larger  value? 
the  active  (central)  flow  possesses  considerably  greater  mass  than 
the  passive  (values  of  "y"  are  written  on  every  curve  in  Fig.  3), 
and  the  low-pressure  peripheral  flow  cannot  show  an  essential 
effect  on  its  nucleus,  as  it  is  three  times  lesser  in  mass.  In 
spite  of  large  nonuniformity  when  =  0.4,  the  average  T#  as  a 
whole  on  cross  section  when  ttq  =  1.52  is  greater  than  when  ttq  £  1.25, 
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and  comprifes  780°K  (instead  of  700°K  when  7io  =  1.25),  which  cannot 
be  said  about  the  flows  which  are  mixed  when  If  =  0.6*4.  In  the 
latter  case,  the  mean  temperature  at  all  values  of  ttq  (1.0;  1.215 
and  1.53)  is  approximately  identical  and  ~560-570°K,  but  the  field 
T*  is  considerably  more  uniform  with  small  ttq.  Hence  the 
conclusion  about  the  advisability  of  an  increase  of  tto  (for  the 
purpose  of  obtaining  larger  average  value  of  T^)  with  an  increase 
in  the  drop  in  temperatures  (by  decreased  $#)  asserts  itself  and, 
contrarily,  about  the  advisability  of  a  decrease  of  7tq  with  a 
decrease  in  the  difference  between  the  temperatures  of  the  initial 
flows  (Fig.  3). 

In  the  process  of  the  experiment  numerous  starts  under 
identical  conditions  of  the  initial  flows  at  the  inlet  to  diffusers 
of  different  geometry  (see  page  70)  were  conducted.  In  this  case 
it  was  noted  that  using  identical  available  gas-dynamic  parameters 
with  a  decrease  in  the  slope  of  the  diffuser  (i.e.,  with  a  decrease 
in  the  values  of  n  and  a  )  increases  the  warming  up  of  its 
peripheral  areas,  i.e.,  increases  the  uniformity  of  the  field  T#, 
which  indicates  the  preferability  of  smaller  values  of  n  (in  the 
range  n  =  *4-3). 

In  Fig.  *4  attention  is  drawn  to  the  transformation  of  fields 
p*  both  during  the  change  in  the  values  of  tt  and  during  a  change 
in  the  values  of  $#.  Moreover  it  is  noticeable  that  with  an  increase 
in  the  preheating  of  the  active  flow,  fields  p#  acquire  a  more 
uniform  nature  even  at  invariable  value  tt  ,  which  speaks  about  the 
intensification  of  the  process  of  mixing  in  diffuser  with  a 
preheating  increase. 

Furthermore,  a  determinate  characteristic  law  governing  a 
change  in  the  structure  of  fields  p#  with  an  increase  of  the 
available  pressure  differential  was  established. 


With  the  retention  of  the  constant  value  of  the  available 
drop  in  the  temperatures,  it  was  noted  that  with  an  increase  of  tto 
the  nonuniformity  of  the  field  p*  sharply  grows  initially.  For 
example,  when  **  *  0.8  =  const,  with  an  increase  of  TrQ  from  1.0  to 
a  value  of  1.3  the  nonuniformity  of  the  field  p*  grows  from  2.4 
to  7.7%,  but  with  further  increase  of  tt^,  the  nonuniformity  remains 
approximately  at  the  same  level. 

The  indicated  transformation  of  fields  p*  was  repeated  with 
the  blasting  of  diffusers  differing  from  aforementioned  in  their 
geometry. 

Thus,  the  geometry  of  diffuser  did  not  exert  a  substantial 
influence  on  the  transformation  of  fields  p#  during  a  change  in  the 
available  gas-dynamic  parameters. 

§  2.  A  Change  in  the  Pressure  of  Mixed 
Flows  Along  the  Diffuser 

Graphs  5-8  present  the  results  of  the  measurements  of  static 

pressure  at  the  internal  and  external  cowlings  (07  and  01  in  Fig.  2) 

along  the  investigated  diffusers.  The  flow  of  the  mixed  flows 

5 

passed  with  numbers  Re  =  (10-6.5)*10  ,  i.e.,  in  self-similarity 
according  to  the  coefficient  of  friction  of  the  area. 

For  each  of  the  indicated  charts  it  is  distinctly  evident 
that  an  increase  in  the  static  pressure  (but  this  means  also  a 
decrease  in  the  velocity  of  the  mixed  flows)  occurs  most  intensely 
in  the  initial  part  of  the  diffuser  -  at  a  length  less  than  the 
half  of  the  entire  length  of  the  investigated  diffuser.  Subsequently 
an  increase  in  the  static  pressure  is  sharply  retarded  in  all 
investigated  systems  and  with  any  diffuser  geometry.  At  compara¬ 
tively  low  speeds  of  the  order  of  =  0.3-0. 5  (Fig.  6)  at  the 
inlet  to  the  same  diffuser,  the  static  pressure  grows  only  up  to 
a  determinate  cross  section,  which  corresponds  to  expansion  ratio 
=1. 7-2.0,  and  then  remains  constant.  But  on  one  of  the  systems 
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(  TT. 


1.22;  ##  =  0.49)  at  the  end  of  the  diffuser  even  a  definite 


drop  in  the  static  pressure  was  fixed  along  the  external  cowling. 


pressure  distribution  along  a  diffuser 

of  an  AMC  with  a  =  12°  and  n  =  2.5. 

np 

KEY:  (1)  Distribution  of  P_..  along 


CM 


the  internal  surface  of  the  diffuser; 
(2)  Distribution  of  P _  along  the 


CM 


external  surface  of  the  diffuser. 
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Fig.  6.  The  projection  of  distribution 
of  pct  along  diffuser  of  an  AMC  with 

a  =  12°,  n  =  2.5. 
np  ’ 

KEY:  (1)  Along  the  internal  surface  of 

the  diffuser;  (2)  Along  external  surface 
of  the  diffuser. 


a  =  16°  and  n  =  3.2. 
n  p 

KEY:  (1)  Internal  diffuser  cone; 

(2)  External  diffuser  cone. 


Fig.  8.  The  projection  of  distribution 
of  pcT  along  the  diffuser  of  an  AMC  with 

a  =  20°,  n  =  iJ.O. 
rip  * 

KEY:  (1)  Internal  diffuser  cone; 

(2)  External  diffuser  cone. 
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This  is  explained  by  the  fact  that  beginning  with  a  certain, 
sufficiently  low  value  of  velocity,  the  increase  in  the  static 
pressure,  which  appears  as  a  result  of  deceleration,  is  compensated 
for  by  the  increasing  losses  of  pressure,  but  with  the  further 
drops  in  velocity  proves  to  be  even  less  than  those  losses  in  the 
recently  noted  case. 

It  was  also  established,  that  for  every  divergence  angle  of 
a  diffuser  there  is  its  limiting  value  of  the  expansion  ratio 
which  corresponds  to  a  maximum  increase  in  static  pressure.  This 
limiting  value  of  the  expansion  ratio  depends  also  on  the  value  of 

X1‘ 

It  is  remarkable,  that  the  limiting  value  of  the  expansion 

ratio  depends  on  the  parameters  of  mixing  and  in  the  entire  range 

of  the  investigated  parameters  (tto  =  1.0-1.55;  =  0.65-0.49; 

A., /A  ,  =  1.0-5. 8)  it  remained  constant. 

Al  ni 

Since  an  Increase  in  the  static  pressure  occurs  most  intensely 

in  the  Initial  part  of  the  diffuser,  then,  in  order  to  decrease  Its 

overall  sizes  during  the  design  one  ought  to  select  an  expansion 

ratio  considerably  less,  namely:  for  diffusers  with  a  =  12°  - 

the  value  n  =  1.5;  for  diffusers  with  a  =  16°  -  the  value 

*  np 

n  =  2.0;  for  diffusers  with  a  =  20°  -  the  value  n  =  2.2,  with 

np  ’ 

which  the  pressure  Increment  ~ 90—  95%  of  the  maximum.  Such  values 
of  expansion  ratios  correspond  to  T  =  2. 3-1.6  calibers  (see  Figs. 
5-8). 

Moreover  the  increase  of  static  pressure  at  the  beginning  of 
the  diffuser  occurs  faster,  the  greater  A^,  however,  the  limiting 
mode  with  greater  A^  begins  later. 

Of  special  interest  is  the  flow  of  the  mixed  flows  in  diffusers 
with  a  steeper  opening  (Figs.  7  and  8). 


If  the  flow  in  a  diffuser  with  anp  *  12°  (Figs.  5  and  6)  in 

the  entire  range  of  the  investigated  values  of  X-^  and  of  the 

parameters  of  mixing  (tto,  *#,  X^/X^)  was  distinguished  by  the 

nonseparable  nature  on  the  entire  length  of  the  diffuser  (with 

different  expansion  ratios);  then  this  cannot  be  said  about  the 

flow  of  the  mixing  flows  in  diffusers  of  the  same  length,  but 

with  a  =  16°  and  20°  (Figs.  7  and  8).  On  the  basis  of  the 
np 

analysis  of  the  available  experimental  data  the  conclusion  was 
made  that  with  the  purpose  of  the  prevention  of  separate  systems 
it  is  not  recommended: 

a)  to  use  diffusers  with  a  >  20°  for  the  mixing  in  them  of 

np  — 

heterogeneous  flows  with  X^  >  0.7; 

b)  to  apply  in  diffusers  with  anp  =  16°  available  pressure 
differentials  of  >1.5  and  X^  >  0.8  (see  Fig.  7). 

§  3.  Characteristics  of  the  Diffuser 
According  to  Pressure  Differential 
of  the  Mixed  Flows 


The  most  essential  diffuser  characteristic  is  the  determination 

of  the  path  losses  which  were  estimated  on  the  experiments  carried 

out  on  the  value  of  the  coefficient  of  the  rise  of  total  pressure 

.'■re s sure  -  a*.  They  also  grew  on  all  operating  modes  of  the 

experimental  device  and  with  any  geometry  of  the  investigated 

diffuser  with  an  Increase  tt  of  loss.  Moreover  within  the  limits 

o 

up  to  the  value  ttq  =  1.2,  losses  remained  constant,  a  result, 
mainly,  of  friction  and  diffuser  flow  expansion. 


With  an  Increase  of  tt  more  than  the  value  1.2,  the  increase 


n  losses  (decrease  of  a#)  becnnir  more  essential  due  to  the  greater 

ilxirig,  which,  at  a  =  12°, 
t3’  *  np  5 


■  i 


Influence  of  the  specific  losses 

were  2%  when  ttq  =  1.5  (in  comparison  with  the  process  when 

tt  =  1.0-1. 2)  and  3.5$  when  tt  =  1.75.  The  described  phenomenon 
o  o 

was  explained  by  the  fact  that  with  an  increase  of  tt  an  increase 
in  the  difference  of  the  initial  flow  rates  (see  Figs.  9  and  10) 
was  observed.  The  essential  difference  in  X^  mixed  flows  with 
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large  n  produced  a  considerable  stagnation  of  an  active  flow  and 

as  a  result,  the  added  losses  are  analogous  to  the  losses  of 

collision  of  "inelastic  spheres."  These  losses  increased  with  an 

increase  of  the  expansion  ratio  and  cxnp  of  the  diffuser.  In  the 

experiment  the  degree  of  irregularity  of  the  parameters  on  a 

pressure  differential  was  determined  also.  As  can  be  seen  from 

Figs.  9  and  10,  the  coefficient  of  pressure  irregularity  grows 

more  or  less  considerably  to  the  value  TrQ  =  1.4,  and  with  a 

further  increase  of  it  to  1.8  it  remains  almost  constant.  This 

□ 

derivation  is  in  complete  agreement  with  the  analysis  of  the 
fields  p*  (see  §  1  this  chapter).  During  the  tests  of  diffuers 
of  different  geometry,  the  maximum  value  of  did  not  exceed 
values  of  1.06-1.075  (nonuniformity  =  6-7.5 %)  when  ttq  =  1.5-1. 8. 
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Fig.  9.  Characteristic  of  a  diffuser 

of  an  AMC  with  a  =  12°,  n  =  2.5  on  tt 

np  *  c 

KEY:  (1)  For  mixed  flow;  (2)  For 

separate  flow. 
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Pig.  10.  Characteristic 

of  the  diffuser  of  an 

AMC  with  a  -16°, 
np  * 

n  =  3*2  on  tt  . 

o 

KEY:  (1)  For  a  mixture; 

(2)  For  separate  flows. 


The  greatest  interest  in  characteristic.;  on  a  pressure 

differential  is  the  comparison  of  specific  imoulses  at  the  inlet 

and  outlet  from  the  diffusers  in  question.  As  can  be  seen  from 

charts  9  and  10,  if  the  value  of  output  impulse  falls  sufficiently 

sharply  with  an  increase  of  then  the  value  of  input  impulse 

remains  almost  constant  on  tto,  which  is  completely  understandable, 

since  the  value  of  J  ,  was  affected  neither  by  loss  nor  the 

yfll 

dissipation  of  energy  as  a  result  of  the  exchange  of  the  impulses 
of  the  mixed  flows. 


The  loss  of  the  output  impulse  (in  comparison  with  input),  at 
any  degree  of  preheating  of  the  active  flow  (it  is  equal  as  without 
the  same)  with  an  increase  of  tt o ,  increased  in  diffusers  of 
different  geometry  (Figs.  9  and  10),  reaching,  when  it  =  1.5-1. 6 
and  X m  =  0.85-0.9,  the  significant  magnitude  of  13-18$.  Such  a 
combination  of  the  limiting  values  of  and  X ^  in  the  given 
experiment  clearly  spoke  of  the  Inexpediency  of  their  application. 
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And,  the  lesser  the  given  available  pressure  drop  in  the  experiment 

the  lesser  was  the  loss  of  impulse.  However  this  decrease  was  not 

continuous  for  it  ,  rather  it  had  a  certain  optimum  with  ttq  which 

differed  from  one,  namely  when  tt  =1.1. 

o 

Hence  the  conclusion  is  asserted  that  the  most  advisable  of 
all  considerations  (losses,  the  momentum  conservation,  uniformity 
of  the  parameters,  etc.)  is  a  diffuser  with  the  mixing  of  flows 
during  low  ttq  z  1.1.  Moreover  the  advantage  of  the  application 
of  this  tt  t  is  greater,  the  lesser  i*. 

§  Characteristics  of  a  Diffuser 
of  an  AMC  According  to  Drop  in 
Temperatures 

With  an  increase  of  the  temperatures  drop  of  the  initial 
flows  the  losses  of  pressure,  the  value  of  tt3  ,  t*^,  and  t*  remained 
virtually  constant. 

As  can  be  seen  from  Pigs.  11  and  12  and  the  range  of  temperature 
drops  =  0.6-0. 4  is  most  definite  from  the  viewpoint  of  change 
in  the  basic  gas-dvnamic  parameters  of  mixture. 

The  consequence  of  this,  obviously,  is  also  the  increase  in 

specific  impulse  of  25-^0!?  with  an  increase  of  temperature  drop 

in  the  aforementioned  limits.  Furthermore,  with  an  increase  of 

temperature  drop  the  loss  of  outlet  pulse  noticeably  decreases  in 

comparison  with  the  inlet.  As  can  be  seen  from  Fig.  11,  this 

difference  in  the  values  of  impulses  for  tt^  =  1.0  is  decreased  from 

the  value  of  9.5%  when  =  1.0  to  2.7 %  when  =  0.^7.  This 

comparison  was  conducted  for  the  case  of  high-speed  mixing  which 

occurred  when  X^  z  Xnl  =  0.7  in  a  diffuser  with  geometric 

parameters  anp  =  16°  ,  n  =  3.2.  The  picture  of  a  change  in  the 

relationship  of  inlet  and  outlet  impulses  in  a  temperature  drop 

becomes  even  more  favorable  at  lesser  (than  normally  used)  initial 

velocities  of  the  mixed  flows.  So,  with  XA,  =  X  ,  -  0.35  and  tt  =1.0 

A1  nl  o 
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(see  Fig.  12)  the  loss  of  outlet  impulse  in  the  same  diffuser 

(a  =  16° )  in  comparison  with  the  inlet  for  an  isothermic  mixing 
n  p 

comprises  no  longer  9*5/?»  hut  only  2%.  With  an  increase  of  the 

temperature  drop  this  difference  sharply  decreased,  and  at 

=  0.88  the  values  of  both  Impulses  are  equalized.  The  graph 

of  Fig.  12  is  comprised  of  relative  values  where  the  horizontal 

J  =  1.0  corresponds  to  the  equality  of  inlet  and  outlet  impulses, 

while  the  curve  J  =  J  A/J  ,  represents  a  change  in  the  outlet 

yfl  y^9  yfll 

pulse  relative  to  inlet.  With  a  further  increase  in  the  temperature 
drop  the  value  of  outlet  impulse  begins  to  prevail  over  its  inlet 
value,  in  which  the  advantages  of  the  process  mixing  are 
effectively  realized.  This  increase  of  the  output  pulse, 
continuously  increasing,  reaches  at  =  0.5  values  of 


Fig.  12.  Characteristics  of  a 

diffuser  of  an  AMC  with  a  =  16°, 

np 

n  =  3.2  according  to 

KEY:  (1)  Note:  The  curve  of  the 

change  in  7t3  when  ttq  =  1.0 

virtually  coincides  with  the  curve 

of  when  t  =1.0. 
o 

Such  a  change  in  the  relationship  of  inlet  and  outlet  impulses 
can  be  explained  in  the  following  form.  The  investigated  diffuser 
mixing  chamber  is  a  two-contour  system  in  which  the  energy  exchange 
between  contours  is  realized  with  the  mixing  of  the  gas  flows. 

Thus  the  laws  governing  energy  exchange  in  the  model  in  question 
are  based  on  the  theory  of  energy  exchange  in  DTRD  [33.  As  can  be 
seen  from  Fig.  12,  the  speed  of  mixed  flow  is  always  less  than  the 
speed  of  separate  flows  at  the  mixing  chamber  inlet.  However,  with 
an  increase  of  the  temperature  drop,  this  difference  in  the  speeds 
continuously  grows  on  the  strength  of  the  fact  that,  with  a  change 
of  from  one  to  0.49,  X A1  increases  by  44%,  while  Xnl  and  X^ 
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barely  change.  But  this  means  that  with  an  Increase  0  ■  1/**, 
thrust  efficiency  increases  and,  as  a  result,  so  does  the  output 
pulse  of  the  retarded  mixed  flow. 

Thus,  an  increase  in  the  temperature  drop  between  the  initial 
flows  does  not  worsen  the  operation  of  the  diffuser  mixing  chamber, 
but  for  the  impulse,  conversely,  substantially  improves  it.  Hence, 
it  is  recommended,  other  conditions  being  equal  (especially  with 
ttq  =  1.0),  to  select  an  operating  mode  of  the  diffuser  mixing 
chamber  in  the  region  of  increased  temperature  drops  (low  values  $■*). 

§  5*  The  Effect  of  the  Geometry  of  A  Diffuser 
of  an  AMC  on  the  Efficiency  of  its  Operation 

It  is  easy  to  note  from  Figs.  13  and  14,  that  a  drop  in  the 

total  pressure  in  a  diffuser  with  smoother  expansion  (a  =  12°) 

np 

proves  to  be  substantially  less  (a*  is  large),  than  in  a  steeper 

diffuser  (otnp  =  20°).  This  advantage  with  respect  to  geometry 

turns  out  to  be  more  essential  with  a  smaller  value  of  and 

larger  tt^  (dot-dash  lines  of  a*  on  graphs  13,  14),  but,  this  means 

also  during  low  values  of  y:  thus,  when  tto  =  1.50;  =  0.64  and 

the  different  X^,  ^nl  Pressure  loss,  upon  transition 

to  a  diffuser  with  a  lesser  slope  of  the  channel,  become  less  than 

approximately  on  7%,  while  when  tt  =  =  1.0  and  XA1  z  X  ,  (solid 

o  Ai  m 

lines  of  a*  in  Fig.  13)  -  it  is  less  only  by  3.4$,  i.e., 
approximately  double. 

As  a  consequence  of  the  lower  losses  p*  and  X  with  a  decrease 
of  the  geometric  parameters  of  the  diffuser,  the  impulse  of  outward 
flow  also  increased  somewhat:  by  4-9$.  Moreover  during  lesser 
ir  (it  =  1)  momentum  grows  more  intensive  than  during  larger 
TT  (TT  =  1.5)  • 

O  0 


ft 
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Pig.  13.  The  effect 
of  the  geometry  of  a 
diffuser  of  an  AMC 
on  its  characteristics 
when  pj^  *  2.0  =  const 

and  XA1  *  const. 

KEY:  (1)  Parameters  of 
separate  flows  into 
Section  I-I. 
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Fig.  14.  The  effect 
of  the  geometry  of  a 
diffuser  of  an  AMC 
on  its  characteristics 
when  p#-^  =  2.0  =  const, 

*A1  =  cons1:  and 
=  0.5  =  const. 

KEY:  (1)  Note:  The 

curve  of  the  change 

in  tt  when  tt  =1.0 
3  0 

virtually  coincides 

with  curve  of 

with  tt  =  1.0. 
o 
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The  preferability  of  a  diffuser  with  smaller  values  of  n  and 
<*np  (within  the  limits  of  the  investigated  series  of  geometric 
parameters)  also  ensues  from  the  comparison  of  the  values  of 
specific  impulses  at  the  inlet  (heavy  line)  and  at  the  outlet  (fine 
line  in  Pig.  13)  from  a  diffuser,  depending  on  its  geometry. 
Moreover  in  a  diffuser  of  any  geometry  with  mixing  with  high  speeds 
*A1  “  °* 7-0.9,  the  losses  of  specific  impulse  were  minimal  when 
it  =1.0  and  *•  =  0. 6h  (broken  lines  in  Pig.  13)  and  they  composed 
in  a  diffuser,  with  n  =  2.5  and  anp  =  12°,  altogether  only  3% » 
i.e.,  the  advisability  is  asserted  for  diffusers  of  any  geometry 
of  low  pressure  drops  and  large  temperature  drops. 

As  can  be  seen  from  Pig.  1*J,  this  advisability  becomes  even 

greater  at  low  initial  velocities  which  correspond  to  A^  =  0.3* 

The  value  of  J  at  such  velocities  and  when  tt  =  1.0;  5*  =  0.5, 
ya  o 

in  the  entire  range  of  geometric  parameters  is  substantially  more 
than  one,  i.e.,  in  a  diffuser  of  an  AMC  of  any  geometry  using  the 
aforementioned  gas-dynamic  parameters  outlet  specific  impulse  is 
always  greater  than  the  inlet  by  a  value  of  3-5  to  5-5%  (in 
proportion  to  a  decrease  in  the  given  divergence  angle).  Both  in 
the  case  of  the  greater  mixing  velocities  (A^  >  0.7),  and  in  this 
case  when  ttq  >  1.0;  A^  >  Anl,  the  relative  specific  impulse  is 
considerably  less  than  that  when  ir  =  1.0  and  less  than  one,  i.e., 
the  specific  impulse  at  the  outlet  of  the  diffuser  mixing  chamber 
is  substantially  less  than  the  specific  impulse  at  the  inlet  to  it. 
However  this  difference  is  considerably  less  than  with  large  A^  and 
comprises  7 -  5— 5 % . 

Thus,  as  the  optimum  combination  of  the  geometric  and  gas- 
dynamic  parameters  in  this  study  one  ought  to  consider  the 
organization  of  the  mixing  of  heterogeneous  flows  with  the  maximum 
preheating  of  active  gas  jet  when  tt  =  1.0  in  a  diffuser  with 
n  =  2.5  and  a  =  12° . 
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5  6.  Characteristics  of  a  Diffuser  of  an  AMC 
According  to  the  Speed  of  Plow  at  the  Inlet 
to  the  Model 


As  can  be  seen  from  Fig.  15 »  with  an  increase  in  the  speed 
of  the  initial  flows  at  the  inlet  to  the  diffuser  mixing  chamber, 
the  total  pressure  recovery  coefficient  monotonically  decreases, 
while  when  approaching  the  value  XA^  *  0.8,  it  falls  very  sharply. 
In  the  range  of  numbers  X^  *  0.3-0. 5  pressure  recovery  factor 
retains  an  approximately  constant  and  very  high  value  which  is 
equal  to  0.9-0.98;  with  a  drop  in  total  pressure,  it  is  equal  to 
1.0.  The  further  decrease  of  a*  during  an  increase  of  X^  is 
connected  with  the  increasing  values  of  static-pressure  gradient 
along  the  length  of  the  diffuser,  which  is  possible  to  trace  on 
Figs.  5-8. 
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Fig.  15.  Characteristics 
of  a  diffuser  of  an  AMC 
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An  especially  noticeable  drop  in  o#  was  observed  with  pressure 
differentials  greater  than  one  (dash-line  in  Pig.  15)  and  with 
larger  XAr 

As  has  already  been  indicated  above  (§  3),  with  an  increase 

of  i  a  continuous  increase  of  A.,  and  drop  of  A  ,  was  noticed, 
o  Al  nl  ’ 

so  that  with  large  7tq  *  1.7-1. 8  (Pig.  9)  AA1  reached  the  transonic 

value  A..  =  1.0.  But  even  at  smaller  values  of  n  >  1.0  and, 
ai  o 

respectively,  of  A^  *  0.8-0. 9  the  outflow  was  accompanied  by 
development  in  the  inlet  part  of  the  active  flow  of  the  local 
supersonic  flows  (which  was  established  by  separate  sensors  of 
p*)  and  by  an  increase  of  added  losses  in  connection  with  this. 

Also  the  appearance  of  flow  breakaway  from  the  wall  of  channel 
with  tt  >  1.0  and  large  AA1  (see  §  2)  contributed  to  an  increase 
in  the  losses  in  this  area. 

The  comparison  of  inlet  and  outlet  impulses  using  the 
different  parameters  of  mixing  was  of  the  greatest  interest. 

As  a  result  of  analysis,  the  essential  effect  of  preheating 
the  active  Jet  on  an  increase  in  the  outlet  pulse  with  a  decrease 
in  the  velocity  coefficient  at  the  inlet  to  the  chamber  was 
noticed.  This  effect  is  especially  considerably  pronounced  with 
ttq  =  1.0.  Thus,  if,  with  the  mixing  of  both  "cold”  flows  and 
it  ■  1.0  (solid  line  in  Fig.  15),  the  loss  of  the  outlet  impulse 
in  comparison  with  inlet  in  the  entire  range  of  A-^  retained  its 
value  as  a  constant  and  equal  to  11%,  then,  with  the  mixing  of 
heterogeneous  flows  by  temperature,  its  value,  in  the  first  place, 
became  considerably  less  even  with  large  A^  =  0.7  and  the  same 
ttq  (6. 0-6. 5%  Instead  of  ll£),  and  most  important,  this  loss  of 
momentum  no  longer  remained  constant  with  drop  of  A^,  but  sharply 
decreased,  equalizing  at  a  certain  value  of  A^  of  the  inlet  and 
outlet  impulses  (the  point  of  intersection  with  the  straight  line 


7  »  1.0).  This  A1  alignment  with  the  increase  of  the  degree 

of  preheating  of  the  active  jet  was  displaced  to  the  side  of  tr. ' 
greater  values:  up  to  A ^  0.51  when  **  ■  0.^7. 

Further  decrease  of  A,  of  nonisotnermal  flows  with  tt  ■  1.0 

I  o 

gave  no  further  loss,  but  an  increase  of  outlet  impulse  in  comparison 
with  the  inlet,  which  reached  5-5?  with  *  0.^7  and  A^  «  0.35* 

In  connection  with  this  the  significant  effectiveness  of  the 
preheating  of  active  jet  on  the  process  of  mixing  with  low  A^ 
should  be  noted,  which  changes  the  relationship  of  the  inlet  and 
outlet  impulses  from  11%  loss  of  the  latter  with  ♦*  =  1.0  to  5.55E 
increase  with  **  =  0.^7  and  A^  =  0.35.  With  an  increase  of  A-^, 
as  has  already  been  indicated  above,  the  preheating  of  active  jet 
to  a  lesser  degree  affects  a  change  in  the  relationship  of  the 
inlet  and  outlet  impulses.  It  is  obviously  possible  to  explain 
such  a  change  according  to  A1  by  a  change  In  the  relationship 

of  the  specific  gravity  of  the  losses  caused  by  friction,  by  the 
expansion  of  diffuser  channel,  by  the  mixing  of  flows,  and,  with 
large  A1,  also  by  the  appearance  of  local  supersonic  zones  (i.e., 
determined  by  value  a*),  on  one  hand,  and  losses  at  outlet 
velocity,  on  the  other  hand. 

As  follows  from  Fig.  15  during  the  extrapolation  of  all 

available  curves  up  to  A^  =  0.85-0.9  in  the  area  of  transonic 

speeds,  variation  by  temperature  and  pressure  drops  gives  no  effect 

at  all  obviously  in  view  of  the  prevailing  effect  of  the  wave 

losses  of  mixed  flow.  But  in  the  area  of  the  moderate  numbers 

of  \A1  there  is  a  noticeably  larger  advantage  of  the  mixing  of 

flows  with  tt  =1.0  and  4*  <  1.0. 
o 

Conclusions 

The  experimental  research  carried  out  on  the  mixing  of  gas 
flows  in  diffusers  of  AMC  made  it  possible  to  establish  the. 
following : 
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1.  The  uniformity  of  fields  depends  not  only  on  the  value 

of  the  available  pressure  differential,  but  with  it  ■  const  also 

on  the  value  of  the  available  temperature  drops;  also  exactly  as 

the  uniformity  of  fields  Tj*  depends  not  only  on  value  **,  but  with 

♦  *  =  const  also  on  the  value  of  it  .  Moreover  to  the  structure 

o 

of  field  change  ttq  affects  in  larger  degree,  than  change  ♦#. 

In  this  case  the  geometry  of  a  diffuser  in  the  range  anp  *  12-20° 
does  not  exert  a  substantial  influence  on  the  transformation  of 
fields  p^,  T*  with  the  retention  of  the  constants  of  initial 
gas-dynamic  parameters. 

2.  In  the  investigated  range  of  gas-dynamic  parameters  of 

the  given  diffusers  of  AMC:  rro  =  1.0-1. 9,  *  0.4-2. 5,  y  =  0. 2-5.0, 

Xj  =  0.1-0. 9,  completely  satisfactory  values  of  the  variation 

factors  of  the  pressure  and  temperature  of  the  mixed  flow  were 

obtained:  t*  *  1.00-1.08  and  =  1.03-1.14. 
p9  t9 

3.  For  every  value  of  upp  and  of  a  diffuser  of  an  AMC  the 

limiting  value  of  its  expansion  ratio,  corresponding  to  the  maximum 
Increase  in  static  pressure,  which  does  not  depend  on  the  mixing 
parameters  Is  established.  Taking  into  account  this  distribution 
of  pressure  along  the  length  of  a  diffuser  of  an  AMC,  with  the 
purpose  of  a  decrease  in  its  overall  sizes  during  the  design,  it 
follows  to  limit  the  degree  of  its  expansion  to  values  less  than 
maximum,  namely: 


for 

a  = 
np 

12° 

-  to 

the 

value 

of 

n  = 

1.5; 

for 

a  = 
np 

16° 

-  to 

the 

value 

of 

n  = 

2.0; 

for 

a  = 
np 

O 

o 

OJ 

-  to 

the 

value 

of 

n  = 

2.2. 

4.  The  boundary  of  separation  systems  is  established  according 

to  pressure  differentials,  which  in  a  diffuser  of  an  AMC  with 

a  =  12°  does  not  occur  in  the  entire  range  of  the  studied 

parameters:  in  a  diffuser  of  an  AMC  with  a  =  16°  passed  with 

it  >1.5  and  X.,  >  0.8,  while  in  a  diffuser  of  an  AMC  with 
o  —  A1 

a  =  20°  -  at  all  values  of  available  pressure  differentials  and 
np 

Xx  >  0.7. 
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5.  Prom  the  considerations  of  lesser  losses  and  larger 

uniformity  of  the  gas-dynamic  parameters  and  retention  of  impulse 

it  follows  to  consider  the  effect  of  the  mixing  of  heterogeneous 

flows  ir  a  diffuser  of  an  AMC  when  a  *  12°  with  the  small 

np 

7i  ®  1. 0-1.1  as  the  greatest  and  possibly  the  maximum  preheating 

of  the  active  gas  jet.  Somewhat  better  results  were  obtained  with 

values  of  it  closer  to  1.1,  rather  than  to  1.0. 
o  ’ 

6.  In  a  dif^ser  of  an  AMC  with  otnp  =  12°  and  with  acceptable 
pressure  differentials  for  DTRD  (tto  ~  1.0)  and  with  such 
temperatures  ($*  =  0.45-0.5),  the  mixing  becomes  favorable  for  the 
specific  impulse  (but  it  also  means  for  C  ),  beginning  with  value 

j 

X i  <  0.5.  With  a  decrease  of  X^  using  simple  invariable  parameters, 
the  increasing  effectiveness  of  the  preheating  of  active  Jet  for 
the  intensification  of  mixing  is  noted.  At  the  different  speeds 
of  the  initial  flows,  in  relationship  to  inlet  and  outlet  specific 
impulses,  the  level  of  the  value  of  X  of  high-pressure  flow  is 
definite,  but  not  the  difference  in  X  of  the  active  and  passive 
flows,  which  with  XA1  =  const  up  to  relationship  does  not  affect 
virtually  XA1/Xnl  “  6  the  relative  specific  momentum. 
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BENCH  CHARACTERISTICS  OF  THE  NOISE 
OF  TURBOFAN  ENGINES 


V.  G.  Yenenkov 


In  the  article  the  acoustic  characteristics 
of  the  basic  noise  sources  of  turbofan  engine 
during  engine  operation  on  a  stand  are  examined. 
Concepts  about  acoustic  characteristics  are 
given  and  the  methods  of  their  calculational 
determination  are  presented.  Information  about 
the  mechanisms  of  the  excitation  of  noise  is 
given.  The  data  of  experimental  studies  of  the 
noise  of  DTRD  with  bypass  ratios  of  y  =  0-6  are 
analyzed. 


Introduction 


The  noise  which  accompanies  takeoff,  landing  and  flight  of 
contemporary  aircraft,  becomes  the  object  of  the  fixed  attention 
of  aviation  specialists.  In  the  countries  with  a  highly  developed 
aircraft  industry  and  civil  aviation,  in  practice,  standards,  norms 
and  certificates  are  developed  and  introduced  which  determine  the 
permissible  noise  levels  of  aircraft  during  their  takeoff  and 
landing  [1],  The  concept  evolves  that  aircraft,  whose  noise 
exceeds  standardized  values,  are  considered  underfulfilled  in 
design  and  they  are  not  allowed  in  operation. 
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The  increase  in  the  Intensity  of  air  transport,  connected 
with  the  advent  of  heavy  duty  aircraft  and  an  increase  in  the 
number  of  voyages,  inevitably  leads  to  the  need  for  efficient 
combat  with  aviation  noise.  Many  specialists  consider  the  problem 
of  a  reduction  in  aircraft  noise  one  of  the  most  important  and 
complex,  removing  it  to  a  place  immediately  after  flight  safety 
control.  In  connection  with  this  it  is  understandable,  how 
important  it  is  to  be  able  to  correctly  estimate  the  acoustic 
characteristics  of  contemporary  civil  aviation  aircraft  engines, 
since  precisely  this  power  plant  is  the  main  noise  source  during 
takeoff  and  landing. 

In  the  domestic  and  foreign  press  many  works  are  published 
dedicated  to  the  study  of  the  individual  characteristics  of  the 
gas  turbine  engine  noise.  The  questions  of  the  evaluation  of  the 
noise  of  turbofan  engine,  which  have  different  bypass  ratios  and 
especially  large  y  =  4-6,  are  not  examined  in  them,  as  a  rule. 

In  this  article  primary  attention  is  given  to  the  acoustic 
characteristics  of  DTRD  with  y  =  0-6. 

§  1.  The  Basic  Concepts  About  the  Acoustic 
Characteristics  of  the  Noise  Sources  of 
DTRD 

The  acoustic  characteristics  of  the  noise  sources  of  engines 
are  described  by  [2]: 

1)  complete  acoustic  power; 

2)  sound  intensity  (by  sound  pressure  or  received  noise  levels) 

3)  noise  spectrum; 

4)  directional  characteristic  of  the  noise; 

5)  duration  of  the  effect  of  noise  and  by  the  frequency  of 
its  recurrence. 
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The  acountic  (sonic)  power  W  is  the  complete  sound  energy 

radiated  by  the  noise  source  into  the  surrounding  space  per  unit 

of  time.  Acoustic  energy  is  frequently  expressed  in  the  levels 

of  acoustic  power  L.,  expressed  in  decibels  relative  to  the 

w  _  1 2 

conditional  threshold  value  of  Wc  *  10  W: 

*.  =  H%5r(dB).  (1) 

The  complete  acoustic  power  is  a  fundamental  characteristic 
of  a  noise  source  as  a  generator  of  sound  energy. 

2 

The  intensity  or  force  of  sound  I  (W/m  )  is  the  energy  content 

transferred  in  a  free  field  by  acoustic  wave  for  one  second  across 

2 

an  area  per  1  cm  ,  normal  to  the  direction  of  motion  of  the  wave. 
Sound  intensity  is  expressed  also  in  intensity  levels 

L  —  10lg~  (dB) , 

’  0 

p 

more  frequently  represented  through  sound  pressure  p  (N/m  )  in  the 
form 

/.  =  10lg/  =  201g£  (dB)  (2) 

••  Pu 

since 

Here  p  -  air  density:  a  -  speed  of  sound  in  air:  p  =  2*10-^ 

2  o  o  o 

N/m  -  the  conditional  threshold  value  of  sound  pressure  which 

-12  2 

corresponds  to  threshold  Intensity  I  =  10  W/m  . 

Value  L  is  called  the  sound  pressure  level  and  is  the  relation¬ 
ship  of  the  operating  pressures  or  sound  intensity  to  their 
conditional  threshold  values. 

The  sound  Intensity  determines  the  energy  of  sound  vibrations 
at  different  points  of  the  acoustic  field.  For  nondirectional 
source  I  =  W/F,  for  a  directed  noise  source 
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(3) 

t 

Here  P  -  the  surface  through  which  the  sound  emission  arises,  and 
the  value  0  characterizes  the  directivity  of  the  radiation. 

Noiee  spectrum  -  the  totality  comprising  the  sound  pressure 
levels  Lj  (levels  of  acoustic  power  Lw^)  obtained  during  the 
frequency  response  analysis  of  the  investigated  noise,  which  is 
the  description  of  the  sound  pressure  or  power  levels  ol  the 
individual  components  of  the  noise  signal  in  different  frequency 
bands  ('octave,  1/2  octave,  1/3  octave).  An  example  of  the  noise 
spectrum  of  a  compressor  is  given  in  Fig.  1.  With  the  noise 
spectrum,  by  means  of  logarithmic  averaging,  it  is  possible  to 
determine  total  (overall)  sound  pressure  level  of  the  complex 
noise  signal  in  the  entire  frequency  range 


To 

4=t01g2;i0  (dB). 
/-I 


Fig.  1.  The  noise  spectrum  of  the  compressor 
of  a  DTRD  with  y  =  1.0  at  a  distance  r  =  30  m, 
at  an  angle  <p  =  10°  to  the  axis  of  the  engine 
inlet . 

KEY:  (1)  Total  levels;  (2)  r/min;  (3)  Hz;  (4)  dB. 
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Besides  the  evaluation  of  total  noise  level,  the  spectral 
characteristic  of  the  noise  signal  is  necessary  for  the  determination 
of  the  level  of  received  noises  -  the  subjective  criterion  which 
considers  the  noisiness  of  the  signal  depending  on  frequency 
make-up  of  the  noise  heard. 


The  level  of  reoeived  noiae  is  measured  in  PN  dB  and  is 
determined  by  the  calculation  according  to  noisiness  in  separate 
frequency  bands  expressed  in  "noises"*  (Pig.  2).  The  total 
noisiness  of  complex  sound  is  equal  to 


\i  » 


where  N  -  the  value  of  noisiness  in  the  most  noised  frequency 

M  <3  H  C 


n 


band 


;  ^  -  the  sum  of  noisiness  in  all  bands;  F  -  the  coefficient 


i-1 


which  considers  the  bandwidth  of  the  filter  used  during  measurement. 
For  a  bandwidth,  equal  to  1/3  octave,  F  =  0.15;  1/2  octave  F  =  0.2 
and  of  1  octave  F  =  0.3« 


M/W  dB)  —  40  33.3 Ig  /Vortlll. 


(5) 


Directional  charaoteriatio  -  the  distribution  of  sound  pressure 
levels  to  a  long  acoustic  field,  according  to  different  radial 
directions  from  the  sound  source,  is  constructed  in  polar  coordinates 
both  for  the  total  noise  (Fig.  3,  [3])  and  for  its  components  in 
different  frequency  bands  (Fig.  4  [^]).  The  long  sonic  field  of 
noise  source  is  understood  to  be  the  space,  in  which  the  sound 
pressure  decreases  inversely  proportional  to  the  distance  from  the 


•This  is  a  reference  to  an  unknown/unidentifiable  sound 
measurement  unit  which  is  obviously  taken  from  the  English  word 
noise.  Further  uses  will  be  placed  in  quotes  to  distinguish  them 
as  measurement  units,  rather  than  the  term  for  unwanted  sound. 
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Fig.  2.  The  dependence  of  noisiness  in 
"noise's"  on  sound  pressure  level  per 
octave . 

KEY:  (1)  "Noise's";  (2)  The  sound 
pressure  level,  dB;  (3)  Hz. 


Fig.  3.  The  directional 
characteristic  of  the 
noise  of  an  aircraft  with 
jet  engines. 


center  of  radiation  to  the  measuring  point.  Prom  the  directional 
characteristic  the  factor  of  directivity  of  10  lg  0  (dB)  is 
determined,  that  is,  the  difference  between  sound  pressure  levels 
at  a  certain  point  of  space  from  the  investigated  source  and  from 
a  fictitious  nondirectional  sound  source  of  the  same  power. 
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(  7 )  J)va?pa#Ma  >tarrpa6/rew<ocmt/  tvyta 
cedufoi/  oKmcfot  (2400/4600*4) 


Pig.  4.  The  directional 
characteristic  of  noise 
of  the  seventh  octave 
(2400/4800  Hz). 

KEY:  (1)  dB;  (2)  deg; 

(3)  The  radiation  pattern 
of  the  seventh  noise 
octave  (2400/4800  Hz). 


The  duration  of  the  effect  of  noise  is  characterized  by  the 
time  of  heard  -  by  the  time  of  effect  of  the  upper  10  dB  or  noise 
above  the  level  of  90  dB.  During  repeated  irradiations  by  noise 
the  sum  of  the  time  of  sounding  of  upper  10  dB  or  of  the  level 
90  dB  is  taken. 

§  2.  The  Basic  DTRD  Noise  Sources 

Turbofan  engines  have  several  sources  of  noise:  the  jets 
of  the  first  and  secondary  contours,  the  fan  and  the  compressor, 
the  turbine,  the  combustion  chamber  and  assemblies.  Depending  on 
the  bypass  ratio,  operating  condition  parameters,  and  engine  power 
rating,  some  of  them  are  basic  and  determinate,  while  the  noise 
of  others  is  masked  by  the  noise  of  the  basic  ones. 


Figure  5  gives  the  data  of  the  experimental  research  on  the 
noise  of  DTRD  with  small  bypass  ratio  [5].  The  propagation 
pattern  gives  a  demonstrative  picture  of  the  noise  source  orientation 
relative  to  the  engine.  Moving  left  to  right,  the  observer  first 
encounters  the  noise  of  the  compressor  with  maximum  intensity  of 
radiation  at  about  a  30°  angle  to  the  axis  to  the  inlet  which 
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prevails  In  the  entire  intake  sector  on  the  majority  of  systems, 
except  those  close  to  maximum,  when  the  noise  of  the  Jet  urowns 
everything.  Moving  to  the  rear  sector  and  passing  through  the 
short  section  of  the  middle  part  of  the  engine,  irradiated  by 
the  noise  of  the  assemblies  and  combustion  chamber,  the  observer 
enters  the  field  of  the  superimposed  radiation  of  two  separate 
noise  sources.  This  is  the  turbine  noise  and,  as  is  characteristic 
for  a  bypass  engine,  the  noise  directed  back  from  the  compressor 
(fan),  bearable  through  the  tract  of  the  second  contour.  The 
propagation  angle  of  this  noise  changes  from  60  to  80°  to  the  axis 
of  the  Jet,  depending  on  the  engine  parameters,  and  its  level  is 
basically  determined  by  the  bypass  ratio.  Finally,  observer  enters 
the  area  of  the  Jet  noise  -  the  most  powerful  source  during 
increased  engine  power  ratings,  with  the  maximum  intensity  of 
propagation  at  a  30-60°  angle  to  the  axis  of  the  Jet.  With  an 
increase  in  bypass  ratio  from  0  to  6,  the  relationships  between 
the  noise  intensities  of  different  sources  and  directional 
characteristic  change  substantially.  The  fundamental  picture  of 
this  change  is  shown  in  Fig.  6.  Especially  attention  is  drawn  to 
Jet  noise  reduction  and  a  simultaneous  increase  in  fan  noise. 


Fig.  5-  The  basic  DTRD 
noise  sources. 

KEY:  (1)  Body;  (2)  Com¬ 

pressor;  (3)  Second 
contours;  (4)  Jets; 

(5)  Turbine;  (6)  Angle 
to  inlet  axis;  (7)  Rela¬ 
tive  powers;  (8)  PN  dB  on 
line . 
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Fig.  6.  The  propogation  pattern  of  the 
DTRD  noise  with  y  =  0,  1  and  6. 

KEY:  (1)  Compressor;  (2)  Engine  axis; 

(3)  Jet;  (4)  Turbine;  (5)  Fan. 


Figures  7-10  give  the  experimental  data  of  the  measurements  of 
the  noise  of  the  engine  which  have  bypass  ratios  0,  1,  2,  3  and  5. 
The  analysis  of  the  directional  characteristic  and  noise  spectra 
visually  shows  that: 

a)  with  bypass  ratio  y  =  0  (Fig.  7)  on  all  engine  power  ratings 
from  the  takeoff  to  landing,  the  specific  noise  source  is  the  Jet, 
its  noise  exceeds  the  noise  of  the  remaining  sources  by  10-20  dB; 

b)  with  bypass  ratio  y  =  1  (Fig.  8),  in  takeoff  mode,  the 
noise  of  the  Jet  (maximum  Intensity  of  radiation  at  an  angle  to 
the  axis  of  the  Jet  0  =  60°)  exceeds  the  noise  of  compressor  by 
4  PN  dB,  and  in  landing  mode  the  noise  of  compressor  directed 
into  front  section  prevails  approximately  by  5  PN  dB ,  with  the 
maximum  intensity  of  propagation  at  an  approximately  20°  angle  to 
the  inlet  axis  engine; 

c)  with  bypass  ratio  y  =  2  (Fig.  8),  at  all  engine  power 
ratings,  the:  compressor  noise  directed  forward  exceeds  the  noise 
of  the  Jet:  in  takeoff  mode  by  6  PN  dB,  on  landing  -  by  18  PN  dB; 

d)  with  bypass  ratio  y  =  3  (Fig.  9),  at  takeoff  engine  power 
rating  the  fan  noise  propogated  from  the  shortened  external  contour, 
already  exceeds  the  noise  of  the  Jet  by  more  than  15  PN  dB  with  the 
maximum  intensity  at  an  angle  approximately  70°  to  the  axis  of  the 
Jet.  In  the  front  sector,  it  is  also  the  greatest.  The  noise  of 
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the  Jet  is  overlapped  by  the  noise  of  the  turbine  (almost  by  10  PN 
dB),  this  Is  a  new  noise  sourcf  ,  for  which  Interest  was  not 
displayed  at  low  bypass  ratios  and  which,  apparently,  will  have 
particular  importance  at  large  y  [6]; 

e)  with  bypass  ratio  y  ~  5,  the  fan  drowns  all  other  sources 
and  in  the  rear  and  front  sectors,  the  directional  characteristic 
of  the  noise  of  such  an  engine  is  givtn  in  Fig.  10  [7]. 


eo*  *  *;• 


-  compressor  noise  radiated  forward 

—  -  compressor  noise  radiated  back 

—  -  —  Jet  noise 

-  total  noise  level 


Fig.  7.  The  directional  characteristic  of  the 
noise  of  turbofan  engine  in  takeoff  and  landing 
systems . 
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Fig.  8.  The  directional  characteristic 
of  DTRD  noise  with  y  =  1  and  2  in  the 
takeoff  and  landing  modes. 

KEY:  (1)  Landing  mode;  (2)  Takeoff  mode; 

(3)  deg;  (4)  dB. 


fO  50  50  70  90  no  t30  '*0 

(6)  '’V**3 

Fig.  9*  The  directional 
characteristic  of  the  DTRD 
ncise  with  y  =  3,  in  tak<-^ff 
mode . 

KEY:  (1)  F-  i;  (2)  Total  noise; 

(3)  Turbine;  (4)  Jet;  (5)  High 
pressure  compressor;  (6)  deg; 

(7)  dB. 
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Pig.  10.  The  directional 
characteristic  of  the  DTRD 
noise  with  y  -  5. 

KEY:  (1)  Modified  TP-39 

DTRD.  Radius  t  *  70  m; 

the  peripheral  speed  of 

the  1st  stage  of  the  rotor 

at  the  periphery  'J  *  308 

nep 

m/s;  (2)  deg;  (3)  dB. 


The  results  of  the  analysis  made  can  be  presented  briefly 
in  the  form  of  a  table  of  the  determinate  DTRD  noise  sources 
dependent  on  the  bypass  ratio  and  engine  power  ratings: 


Table . 


Mode 

ratio 

Takeoff  mode 

Landing  mode 

0 

Jet 

Jet  -  compressor 

1 

Jet  -  compressor 

Compressor  -  Jet 

2 

Compressor  -  Jet 

Compressor 

3 

and  greater 

Compress zr  (fan) 

Compressor  (fan) 

The  comparative  calculated  evaluation  of  the  noise  of  the  fan 
and  the  jet  of  a  DTRD  with  T?  -  1300°K,  tt*  =  17  and  bench  thrust 

j  K  0 

R0  =  8  T  during  a  change  in  y  from  0  to  9,  which  was  carried  out  by 
the  author  in  conjuncticn  with  B.  N.  Mel’nikov,  also  shows  that  at 
y  >  2  fan  noise  becomes  determinate. 
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Pig.  11.  The  effect  of 
bypass  ratio  on  compressor 
and  Jet  noise. 

KEY:  (1)  Compressor; 

(2)  Jet;  (3)  PN  dB. 


§  3.  The  Nature  of  Noise  Sources 


A  noise  source  which  radiates  acoustic  waves  in  all  directions, 
can  be  an  oscillating  sphere  or  monopole  (Fig.  12). 


Two  spheres  in  line  and  which  oscillate  at  antiphase  (one  is 
compressed,  another  is  expanded)  are  a  dipole  sound  source. 

Similar  to  this  model  is  the  radiating  sphere  which  oscillates 
along  a  straight  line.  The  important  feature  of  a  dipole  source 
is  the  directional  effect:  a  large  part  of  the  sound  is  propagated 
in  the  plane  of  the  dipole.  For  dipole  excitation  the  harmonic 
force  w:  .ch  acts  on  the  part  of  hard  boundary  is  necessary. 


Monopole  Dipole  Quadrupole 


Fig.  12.  An  Illustration  of  the  nature  of 
noise  sources. 

If  we  place  together  two  dipoles,  that  oscillate  in  antiphase, 
a  quadrupole  sound  source  will  be  obtained.  In  this  case  there  are 
two  mutually  perpendicular  directions,  every  point  of  which  obtains 
a  compensated  "null"  signal.  At  an  angle  of  ^5°  signals  do  not 
compensate  for  each  other,  since  phases  and  amplitudes  are 
different,  and  in  this  direction  the  intensity  of  radiation  is 
greatest . 

Of  the  enumerated  sources,  the  quadrupole  is  the  least 
effective  noise  generator. 

§  ^4.  The  Characteristics  of  the  Noise 
Created  by  the  Jet 

1.  The  Mechanism  of  the  Jet  Noise  Generation 

The  noise  of  the  Jet  arises  as  a  result  of  the  turbulent 
mixing  of  the  gas  flow,  which  escapes  from  the  jet  nozzle,  with 
surrounding  air.  In  the  zone  of  turbulent  mixing,  which  has 
lengths  of  15-25  diameters  of  nozzle,  fluctuations  of  the  speeds 
and  eddies  appear,  in  connection  with  which  occur  the  pulsations 
of  pressures  and  densities.  Lighthill  [8]  proposed  to  consider 
turbulent  flow  as  a  large  quantity  of  centers  of  quadrupole 
radiation,  driving  at  a  certain  convective  speed  in  an  ambient 
medium.  Actually,  the  oscillating  gas  flow  creates  in  the 
atmosphere  those  density  fluctuations,  which  would  arise  in  a 
uniform  acoustic  medium  under  the  action  of  a  system  of  tensions 
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applied  from  without  to  the  elements  of  the  medium.  Each  element 
of  the  gaseous  medium,  during  the  appearance  of  the  tensions, 
experiences  the  effect  of  two  equal.  In  magnitude,  and  opposed 
forces.  But  any  force,  acting  on  an  element  of  medium,  in 
acoustics  Is  the  equivalent  of  a  dipole,  while  a  pair  equal  in 
magnitude  and  opposed  forces  -  to  a  quadrupole. 

That  means  the  jet  can  be  broken  into  many  volumes,  equal  to 
the  dimension  of  an  eddy,  and  each  of  them  considered  as  a 
separate  noise  source  in  the  form  of  an  acoustic  quadrupole. 

2.  The  Acoustic  Power  of  th;  Jet 

At  present,  the  calculation  of  the  acoustic  power  of  a  Jet, 
in  accordance  with  Lighthill's  theory,  is  conducted  from  the  gas- 
dynamic  and  geometric  parameters  of  the  Jet  in  the  jet  nozzle  area. 
Many  investigations  established  that  the  best  agreement  of  results 
of  calculations  and  experimental  data  for  hot  and  cold  Jets  is 
observed  during  the  calculation  of  acoustic  power  according  to 
formula  [2]1 

(W).  (6) 

By  applying  gas-dynamic  functions  for  the  expression  of  density 
and  velocities  through  the  parameters  of  stagnant  flow  and  the  given 
velocity  X,  let  us  write  formula  (6)  in  the  form 

W  =-.  KH 1 0,)p/  Ta'  ’  |  MW),  ( 7 ) 


‘These  investigations  are  related  to  jets  which  have  subsonic 
and  small  supersonic  speed  of  outflow  (1.5  <  tt  <_  2.2).  At  a 

large  supersonic  speed  of  outflow  the  noise  of  the  Jet  is  proportional 
to  the  discharge  velocity  by  the  cube.  At  low  discharge  velocities 
(c,-  <  300  m/s)  the  acoustic  thickness  of  the  jet  is  proportional 

tO  Cr. 
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where 


o-*V‘ ■(»£)' 


(for  air  B  =  13*65;  for  gas  B  *  12.44). 


Since,  for  the  calculation  of  acoustic  power,  the  parameters 

in  the  outlet  section  of  a  subsonic  nozzle  are  taken,  for  the 

subcritical  system  of  outflow  of  gas  (ti  <  1.89,  Pc-  *  pn)  it  is 

pc  5  0 

possible  to  write 


'=o|  i 


P.W/#  F, 


During  a  supersonic  system  of  gas  outflow  from  a  subsonic 
nozzle  (tt  >  1.89;  p^  >  Pq;  Aj.  =  1.0)  formula  (7)  will  take  the 
form 


*r=,»(sTi 


The  calculation  of  acoustic  power  for  supercritical  systems 

of  outflow  according  to  formula  (9)  agrees  well  with  experiment 

up  to  values  of  tt  <2.2. 
v  pc  — 

The  proportionality  factor  K  in  formulas  (6-9)  has  the 
following  values  : 

a)  for  engines  which  work  at  nominal  and  maximum  ratings, 

K  =  1.5*10”^; 

b)  for  engines  which  work  under  conditions  less  than  0.8 
nominal,  K  =  2.5*10~S 

c)  for  gas  jets  which  escape  from  a  convergent  nozzle, 

K  =  0.8*10'\ 


...  *c«.-j** I  -I  -■  i..- >T  iimm ■  “ 


1 

i 

i 

The  calculations  of  the  acoustic  power  of  the  jets  of  DTRD 
without  the  premixing  of  the  flows  of  the  two  contours,  according 
to  presently  available  experimental  materials,  can  be  produced  with 
the  summation  of  the  acoustic  powers  of  the  Jets  of  the  first  and 
second  contours 

VP,tM  =  W'.+  W'ii. 


By  formulas  (7-9)  the  values  of  acoustic  power  are  given  in 
connection  with  the  dimensionality  of  engine,  determined  by  the 
discharge  area  of  the  jet  nozzle  P^.  It  is  possible  uo  eliminate 
the  effect  of  dimensionality,  by  determining  the  acoustic  power 
per  unit  of  thrust.  For  this  Fc  is  represented  in  the  form 


while  the  formula  of  acoustic  power  after  the  transformations 
appear  in  the  f  wing  manner: 


IP 

R 


Knik  f.  n 
/;/  (  2  / 


i 


v  i  /;/,*)■ 


3.  The  Noise  Spectrum  of  the  Jet 


(10) 


The  noisv,  spectrum  of  a  jet  with  a  subsonic  pressure 
differential  in  the  nozzle  is  continuous  with  a  weakly  marked 
maximum.  During  an  increase  in  pressure  differential  and  an 
Increase  in  the  temperature  of  the  Jet,  the  portion  of  the  high- 
frequency  component  in  the  noise  spectrum,  increases  with  an 
increase  in  diameter,  the  low-frequency  portion  [2]  increases. 


The  experimental  research  showed  that  the  relative  spectrum  of 
acoustic  power  may  be  presented  by  one  curve  in  the  form  of  a 
dependence  of  on  the  Strouhal  number  Sh  *  fd^/c^ 

(Fig.  13).  In  this  case,  the  levels  comprising  the  spectrum,  and 
the  level  of  acoustic  power  Lw  are  connected  by  the  function 
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'•  =  10Ik2>10. 

/  1 

Having  calculated  the  total  acoustic  power  W  by  one  of  the  formulas 
(7-10),  depending  on  the  system  of  the  Jet  outflow  from  the  nozzle, 
then  having  determined  Lw  -  10  lg  W  +  120,  it  always  is  possible  to 
find  the  level  of  acoustic  power  in  1/3  octave  frequency  bands 
with  the  aid  of  the  curve  given  in  Fig.  13. 


*  jo  -s~- 

'■J  ,n  o  o 
*3  60 - 


\U-in 
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1  L-  LI- 1 
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( 2  )Vvc/ro  CmpyiaM  S jjp 

Fig.  13.  The  dependence  of  the  level  of 
acoustic  power  of  a  Jet  on  the  Strouhal 
number. 

KEY:  (1)  Models  with  d^,  in  mm; 

(2)  Strouhal  number  Sh  =  f’d^/v,.; 

(3)  Turbojet.  3  ^ 


006  0.06  0080,1  0.2  06  0.6  03  /0  2.0  6,0  6.0 

(1  )Vuc*o  Cmpyxa/ta  Sh  ■ 

/- 20^0^60';  2-70'i  0*180' r 
( 2  )tfodej!u  c  ds=  1QQ**  npupas*M  Tj 
- ,  •  - 525‘K.  k ,  *  -8?3‘U,  ii.  Ik -TPS 

Fig.  14.  The  dependence  of  the  sound  pressure 
level  of  a  Jet  on  the  Strouhal  number. 

KEY:  (1)  Strouhal  number  Sh  »  f’d^/v^;  (2)  Models 

with  d^  =  100  mm  with  different  Tjl. 
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The  dependence  of  the  Jet  noise  spectrum  (Pig.  14)  on  the 
angular  position  of  observation  point  can  be  reduced  to  the 
consideration  of  two  angular  areas  20°  <  0  <  60°  (curve  1)  and 
70°  <  0  <  l80°  (curve  2). 

4.  Sound  Pressure  and  Received  Noise  Levels 


The  sound  pressure  level,  according  to  the  determination  of 
(2)  and  (3),  is  equal  to 


£=10lg^  +  l01g4>  (dB) 


During  sound  propogation  into  hemisphere,  P  =  2nr  ,  therefore 

L  =!01g  IP -j~  10lg—  -201g/'-f-l0lg«l>  (dB) 
or 

/.- 10lgU/-fl09  -  20!g/-  +  101g«I>(dB).  (12) 

The  factor  of  directivity  10  lg  $  is  determined  according  to  the 
results  of  the  processing  of  experimental  research  data.  Figure  15 
depicts  the  change  of  this  factor  depending  on  the  angle  of 
directivity  for  cold  and  hot  Jets  with  gas  temperatures  up  to 
1000°K  [2],  In  the  direction  of  the  greatest  intensity  of 
radiation,  the  factor  of  directivity  composes  a  value  on  the 
order  of  8  dB. 

The  levels  of  received  noise  are  determined  from  the  Jet  noise 
spectrum,  constructed  for  given  d^  and  c^  in  accordance  with  the 
experimentally  established  distribution  of  the  components  of 
sound  pressure  levels  in  1/3  octave  bands  (Fig.  14). 

For  a  less  ,  trict  evaluation  of  the  levels  of  received  noise, 
it  is  possible  to  recommend  quite  a  simple  method.  As  a  result 
of  the  considerable  similarity  of  spectra  of  jets,  which  have 
different  gas-dynamic  and  geometric  parameters,  their  levels  of 
received  noise  from  sound  pressure  levels  by  a  certain  value 
APN  dB.  The  adjustment  of  the  APN  dB  has  the  following  expressions 
[2]: 


1)  for  a  TJE  and  a  turbofan  engine  without  mufflers 

A  MV  dB  «  8,3  -  0,0083  r  (13) 

2)  for  a  TJE  and  a  turbofan  engine  with  noise-abating  nozzles 
of  low  effectiveness  (3-6  dB  of  a  reduction  in  the  noise) 

A/J.V  dB  =  10,:;  -  0,0083  r  ( in ) 


Fig.  15.  The  dependence  of  the 
directivity  factor  on  the  angle  0. 

KEY:  (1)  Factor  of  directivity 

10  lg  $  dB;  (2)  deg. 


Thus , 


L  (MVdB)= /.  (dB)  +  AMV  dB . 


5.  The  Order  of  Calculation  of  the  Acoustic 
Characteristics  of  the  Jet 


The  order  of  calculation  of  the  Jet  noise  can  be  recommended 
as  follows : 


1.  The  total  acoustic  power  of  the  Jet  W  is  determined  from 
one  of  formulas  (7-10). 
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2.  The  level  of  the  acoustic  power 

Lw  =  l01gU/  +  *20  (dB) . 

3.  The  components  of  the  spectrum  of  acoustic  power  are 
determined  in  1/3  octave  frequency  bands 

Lwl  =  Lw  -j-  A£v)l 

where  -  components  of  the  dimensionless  spectrum  of  acoustic 

power  of  the  Jet  given  in  Pig.  13.  In  this  case  the  frequency 
of  the  spectrum 


4.  The  level  of  the  total  sound  pressure  at  the  points 
arranged  at  a  distance  r  from  the  nozzle  edge  of  the  engine  is 
determined  from  formula  (12). 

5.  The  levels  of  the  components  of  the  spectrum  of  sound 
pressure  in  1/3  octave  bands  are  equal  to 

where  AL^  -  the  components  of  the  dimensionless  spectrum  of  sound 
pressure  given  in  Fig.  14. 

6.  The  level  of  received  noise  of  the  jet  is  determined  from 
the  spectrum  of  sound  pressure  in  1/3  octave  bands  or  from  formula 
(15)  with  the  use  of  the  a^ proximate  relationships  (13)  and  (14). 
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§  5*  The  Characteristics  of  the  Noise 
Generated  by  the  Compressor  (Pan)  of 
a  DTRD 


1.  The  Mechanism  of  Noise  Excitation 
in  the  Compressor  Stage  and  the 
Basic  Determining  Factors 


The  noise,  generated  by  the  compressor,  is  broadband  noise  of 
aerodynamic  origin,  against  the  background  of  which  are  marked 
several  discrete  components  which  correspond  to  the  frequencies 
of  blade  passage  (Fig.  16). 
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Fig.  16.  The  compressor 
noise  spectrum  of  a  DTRD 
with  y  ~  5- 

KEY:  (1)  Modified  TF-39 
DTRD.  Data  1/3  octaves; 

R  »  70  m.  Peripheral 
speed  of  the  1st  rotor 
stage.  U  =  308  m/s 

M  6  U 

0  =  70°  (maximum  PN  dB) ; 

(2)  Frequency  of  blade 
passage  of  the  blades; 

(3)  1st  stage;  (4)  2nd 
stage;  (5)  Hz;  (6)  L  [dB]. 


The  fundamental  components  of  the  noise  of  an  axial-flow 
compressor  are  noise  of  rotation,  eddy  noise  and  noise  from  the 
heterogeneity  of  the  stream  which  flows  around  the  blade. 

The  noise  of  rotation  is  formed  due  to  the  forces  of  the 
effect  of  blades  on  the  flow,  has  a  dipole  nature  and  are  discrete 
in  character.  It  is  connected  with  the  motion  of  the  rotor  wheel 
blades,  which  create  a  periodic  disturbance  of  the  medium  in  the 
plane  of  rotation. 
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Eddy  noise  is  caused  by  the  vortices  which  are  formed  as  a 
result  of  air  flow  around  the  inlet  components  and  rotor  wheel 
blades  and  the  stator.  Eddy  noise  has  quadrupole  nature  and  is 
broadband. 

Noise  from  the  heterogeneity  of  flow  is  produced  by  additional 
forces  on  the  blades  which  appear  as  a  result  of  the  nonuniformity 
of  flow  parameters  in  the  plane  of  rotation  of  the  wheel.  The 
heterogeneity  of  flow  in  the  stage  can  be  caused  by  the  features 
of  the  ir.let  to  the  compressor,  and  also  by  distortions  of  the 
field  of  velocities  after  the  stator  blades.  Noise  from  the 
heterogeneity  of  flow  has  a  dipole  nature  and  is  discrete  in 
character . 

Separate  components  of  compressor  noise  are  currently  not 
studied  sufficiently,  and  thus  far,  there  is  no  strict  theory. 
During  engineering  empirical  and  semi-rational  formulas  are 
utilized  which  relate  the  gas-dynamic,  geometric  and  acoustic 
parameters . 

The  purpose  of  many  very  recent  investigations  is  detailed 
research  on  the  mechanisms  of  noise  excitation  in  the  compressor 
stage  and  isolation  of  its  determinate  factors.  One  of  the  most 
thorough  works  of  this  iescription  is  published  by  Smith  and  Howes 
[5].  The  authors  gave  primary  attention  to  eddy  noise  and  noise 
from  the  heterogeneity  of  flow.  The  mechanisms  of  the  excitation 
of  broadband  noise  examined  by  them  and  the  discrete  comprise 
are  visually  illustrated  by  the  diagrams  A-D  (Fig.  17).  Diagram 
A  gives  a  representation  about  the  random  process  of  changes  of 
circulation,  pressure,  and  lift,  occurring  on  the  profile.  The 
initial  flow  turbulence  considerably  amplifies  this  process, 
therefore  diagram  B  reflects  a  stronger  mechanism  of  noise 
excitation.  Diagram  D  shows  routine  and  periodically  repeated 
changes  in  the  angles  of  attack  of  the  rotor  blades  as  a  result 
of  the  deviation  of  flow  lines  at  the  exit  edges  preceding  stator 
vanes,  which  produces  the  appearance  of  discrete  currents.  The 


most  intense  discrete  noise  appears  during  the  passage  of  the  rotor 
blades  through  the  trailing-edge  wakes  of  stator  grid  (diagram  D) , 
when  a  change  of  the  velocity  in  every  trailing-edge  wake  causes 
a  cyclic  variation  in  the  pressure  on  every  vane  of  the  subsequent 
grid. 


Fig.  17.  The  diagram  of  the  mechanisms 
of  the  excitation  of  noise  in  compressor. 

The  more  detailed  classification  of  the  mechanisms  of  noise 
excitation  in  a  compressor,  described  by  M.  Pianko  [9],  is  given 
in  Fig.  18. 

Thus,  the  basic  factor  which  determines  the  noise  of  a 
compressor  is  the  interaction  between  turbulent  perturbations  and 
the  traces  which  are  formed  after  preceding  vanes  and  the 
subsequent  series  of  blades  of  the  rotor  and  stator.  The  intensity 
of  this  interaction,  to  a  significant  extent,  depends  on  the  axial 
clearance  between  blades  and  the  relative  number  of  rotor  and 
stator  blades  in  each  stage.  In  the  process  of  design  it  is 
possible  to  control  the  noise  of  compressor  by  lowering  the 
broadband  noise  and  determinate  discrete  tones  by  means  of  an 
increase  in  the  axial  clearances  and  the  selection  of  the  optimum 
relationship  of  numbers  of  blades  of  the  stator-rotor  unit. 
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2.  The  Basic  Parameters  Which  Determine 
Compressor  Noise 

Many  researchers  were  inclined  to  consider  the  basic  parameter, 
determining  compressor  noise,  to  be  the  absolute  velocity  of  flow 
in  the  peripheral  cross  section  of  the  rotor.  The  dependences  of 
noise  level  on  peripheral  speed  established  from  this,  evidenced 
a  considerable  change  in  complete  acoustic  power  with  a  change  in 
the  engine  power  rating,  which  does  not  agree  with  experience. 

The  careful  analysis  of  the  results  of  experimental  research, 
made  with  consideration  of  the  features  of  the  mechanisms  of  the 
excitation  of  compressor  noise  by  Smith  and  Howes,  showed  that  both 
the  random  and  regulated  fluctuations  of  the  flow  velocity,  pressure, 
and  lift  on  the  profile  depend  on  the  magnitude  of  the  relative 
velocity  vector  of  the  flow  around  the  blade.  The  specially  set  up 
experiments  confirmed  the  accuracy  of  the  analysis. 

In  the  area  of  the  linear  dependence  of  air  consumption 
through  the  compressor  on  revolutions  (at  invariable  angles  of 
attack),  the  noise  level  changes  according  to  a  law,  close  to  the 
sixth  power  of  the  relative  velocity. 

In  practice  a  difference  in  the  intensity  of  sound  emission 
in  front  and  rear  sectors  is  observed.  The  transfer  of  larger 
portions  of  the  jet  noise  energy  downstream  than  upstream  selves 
as  the  reason  for  this,  the  so-called  effect  of  convection  (if 
there  were  no  consumption  of  the  air  through  the  compressor,  then 
uniform  noise  distribution  between  front  and  rear  hemispheres 
could  be  expected.  They  propose  that  the  portion  of  radiated  energy 
which  is  propagated  upstream,  is  equal  to  0.5  (1  -  M),  and  downstream 
0.5  (1  +  M) ,  where  M  is  equal  to  the  Mach  number  of  the  flow  in 
the  axial  direction.  Smith  and  Howes  give  corrections  for  noise 
convection  of  AF  (Figs.  19  and  20). 
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An  important  parameter  which  determines  compressor  noise  is 
the  flow  turbulence  in  the  stage.  In  a  compound  compressor  the 
rotors  of  the  subsequent  stages  are  located  under  conditions  of 
higher  turbulence,  than  the  first  stage  rotor,  as  a  result  of 
which  the  Intensity  of  noise  excitation  can  grow  by  10  lg  N/N-l 
at  the  stage.  The  effect  of  turbulence  can  be  combined  with  the 
effect  of  convection  and  a  new  correction  8^  for  any  stage  with 
number  N  obtained  (Fig.  21). 


Fig.  21.  Total  correction  for  noise 
convection  in  the  front  and  rear  sectors 
of  a  compound  compressor  (fan). 

KEY:  (1)  Number  of  stages;  (2)  Increase 

In  the  i.'iise  In  rear  sector;  (3)  The 
reduction  of  noise  in  front  sector; 

(4)  Total  corrections  8^,  dB; 

(5)  Relative  M  number  on  the  blade 
periphery . 


The  values  of  angle  of  attack  a  exert  a  great  effect  or.  the 
dependence  of  noise  intensity  on  the  relative  velocity,  however, 
this  problem  has  been  studied  insufficiently. 

The  value  of  relative  axial  clearance  (S/C)  between  stator  and 
rotor  blades  substantially  affects  the  levels  of  discrete  a.  d 
broadband  noise  of  the  compressor,  however,  theoretically  the  laws 
governing  this  effect  thus  far  are  not  developed.  The  experimental 
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data  In  this  area  are  quite  limited,  and  the  available  information 
is  related  to  models  of  a  different  dimensionality.  The  tendency 
toward  a  decrease  in  the  intensity  of  the  fundamental  tone  with  an 
increase  in  the  axial  clearance  is  frequently  emphasized,  although 
the  intensity  of  separate  harmonics  can  increase  in  this  case. 

3.  The  Acoustic  Power  of  a  Compressor 

Work  [12]  is  dedicated  to  the  investigation  of  the  dependence 
of  acoustic  power  of  a  compressor  on  basic  determining  factors, 
in  which  analytical  expressions  of  the  levels  of  acoustic  power  for 
the  cases  of  the  nonseparated  flow  and  vortex  flow  of  the  blades 
are  given.  They  take  the  form  (respectively): 

/.,  =:  08  +  47,ti!g  ^  (-  10IK  (1  ±  /Vf)  +  101-  O  ( dB )  , 

£,  =  118  f -34.41R  I0I«(1  ;_/Vf)+10!g(;(dB)  . 

Here:  w  -  relative  velocity  of  flow  around  the  blades  in  the 
peripheral  cross  section  (m/s);  M  -  Mach  number  at  the  entry  to 
first  stage;  G  -  air  consumption  (kg/s). 

Work  [10]  gives  a  similar  expression  for  the  level  of  acoustic 
power,  in  which  the  term  which  considers  the  factor  of  convection 
is  omitted 

/.,=  130.5  +  501gf  +  IOIgO'(dB)  ,  (17) 

*1 

where  c  -  axial  velocity  at  the  inlet  to  the  compressor  (m/s). 

Formulas  (15-17)  give  the  dependence  of  the  level  of  acoustic 
power  on  the  relative  rate  of  flow  arourd  the  blades  on  the 
periphery,  respectively  in  powers  of  5-76,  *1.^  and  6  (consumption 
is  function  of  velocity  to  the  first  power)  .  The  law  w^  is 
characteristic  for  dipole  noise  sources. 


(15) 

(16) 
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4.  The  Compressor  Noise  Spectrum 

In  the  compressor  noise  spectrum  broadband  noise  and  discrete 
component  (Fig.  22)  are  distinctly  distinguished. 

Broadband  noise,  approximately  to  1000  Hz,  corresponds  to  the 
noise  of  the  Jet,  and  therefore  can  be  considered  as  the  overlapping 
noise  from  the  rear  end  of  the  engine.  In  the  range  of  frequencies 
above  1000  Hz,  the  broadband  component  of  the  spectrum  is  formed  by 
the  sonic  radiation  of  compressor.  In  Fig-  23  large  numbers  of 
arbitrarily  selected  spectra,  constructed  relative  to  the  maximum 
level  in  dB  [53  are  presented.  Deviations  at  lower  frequencies  are 
caused  by  the  difference  in  modes  and  by  the  effect  of  the  introduc¬ 
tion  of  the  jet  noise.  The  characteristic  frequency  of  maximum  has 
a  value  on  the  order  of  5000  Hz,  which  is  connected  with  the 
typical  dimension  of  the  vortices. 


Frequency  at  the  middle  of  bands  of  1/5  octave.  Hz 

Fig.  22.  A  typical  noise  spectrum  of  a  single- 
stage  fan. 
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Frequency  in  the  center  of  band  £  > 
of  1/3  octaves ,  Hz 

Fig.  23.  Determined  spectrum  of  broadband 
noise  at  a  distance  of  30.4  m  along  a  line. 

Discrete  noise  tones  have  frequencies,  multiple  to  the 
frequencies  of  the  passage  of  the  blades  of  the  rotor,  and  have 
basic  and  upper  harmonics  of  the  first  compressor  stages.  The 
fundamental  frequency  has  the  value  f  =  nz,  where  z  -  number  of 
rotor  blades  and  n  -  number  of  revolutions  per  second;  the  second 
harmonic  with  double  frequency  has  a  sound  intensity  6  dB  less 
than  the  fundamental  tone,  the  third  on  3*5  more  dB  less  etc., 

(in  accordance  with  the  law  of  reduction  10  lg  k  ,  where  k  -  number 
of  harmonics). 

5.  The  Sound  Pressure  Levels  and  Received  Noise 

Sound  pressure  levels  in  the  zone  of  maximum  intensity  of 
propagation  along  a  line,  at  a  distance  of  30  m  from  the  source, 
can  be  estimated  separately  for  broadband  noise  and  discrete 
components  by  the  semi-empirical  dependences  proposed  by  Smith 
and  Howes  [5]. 
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A.  Single-Stage  Pan 


The  level  of  broadband  noise  at  the  1/3  octave  band,  which  has 
maximum  intensity  of  propagation,  is  estimated  with  respect  to  the 
formula 

L  =s  50lg  -f  tOlg  0  -j-  AF  -f-  a  -f  M  f  dB  )  ,  ( 1 8 ) 

where  a  -  mean  deviation  of  attack  angle  from  that  corresponding 
to  maximum  lift  (in  degrees);  AP  -  correction  factor  which  considers 
the  expenditure  effect  of  flow  on  the  propagation  of  noise  (dB). 

The  form  of  the  spectrum  is  borrowed  from  Pig.  23  with  a 
frequency  shift  proportional  to  the  ratio  of  the  average  width  of 
the  active  and  guide  vanes  on  periphery  to  that  accepted  in 
investigation  [5]  50  mm. 

The  sound  pressure  level  of  the  fundamental  tone  of  discrete 
component  is  determined  from  the  formula: 

A  =  50lg  —  lOlg  (^)*  +  61  (dB),  (19) 

where  S/C  -  ratio  of  the  value  of  the  axial  clearance  between  the 
stator-rotor  unit  to  the  width  of  the  grid  on  periphery. 

The  sound  pressure  levels  of  the  second  and  subsequent 
harmonics  are  lower  than  the  fundamental  tone,  by  20  lg  k.  The 
fundamental  tone  has  frequency  f  =  nz  and  the  harmonic  frequencies 
are  multiples  of  this  value. 

Discrete  components  of  the  fundamental  tone  and  harmonics  are 
included  in  the  frequency  bands  corresponding  to  them  in  the 
spectrum  of  broadband  noise. 
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B.  Compound  Compressor 


The  maximum  level  of  broadband  noise  in  the  frequency  band  of 
1/3  octaves  is  calculated  from  the  formula: 

/-  =  5°lg-j^-f-101g0+a-t-?s  +  5l  (dB)  ,  (20) 

where  0^  -  correction  factor  which  considers  the  expenditure  effect 
of  the  flow  for  N  similar  stages  (dB). 

The  sound  pressure  level  of  discrete  components  of  the 
fundamental  tone  of  the  noise  radiated  forward  is  estimated 
according  to  the  formula 

/.  =  50lgl^-  10lg(§)V‘W  +  61  (dB).  (21) 

It  is  assumed  to  be  that  every  subsequent  stage  lowers  the  level  of 
discrete  components  of  the  noise  which  is  propagated  back,  by 
approximately  5  dB.  The  correction  factor  NAF  in  this  case  is 
replaced  by  [AF  -  5(N'  -  1)];  AF  is  taken  from  Fig.  19,  and  N'  is 

the  number  of  stages  counted  from  the  end  of  the  compressor. 

The  levels  of  received  compressor  noise  can  be  determined  by 
a  noise  spectrum,  constructed  in  accordance  with  the  given  procedure. 

M.  Pianko,  in  work  [9],  gives  semi-rational  formulas  for  the 
calculation  of  the  maximum  level  of  the  received  noise  in  PN  dB, 
determined  at  a  distance  300  m  from  the  source. 

i (f'X  dB)  -  IOIjv  R—  46ijj  s?  -  10lg  4  —  10lg —  15. 


Here  r  -  outside  radius  of  the  compressor  (m) ;  R  -  thrust  (given); 

H 

R  -  specific  thrust  (kg/kg-s). 
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6.  Order  of  Calculation  of  the  Compressor 
Noise 

The  following  can  be  recommended  for  the  order  of  calculation 
of  the  noise  of  compressor: 

1.  The  maximum  level  of  broadband  noise  in  a  1/3-octave 
frequency  band  is  determined  from  formulas  (15)  or  (17). 

2.  The  spectrum  of  broadband  noise  is  constructed  in 
accordance  with  the  results  of  the  processing  of  a  large  number 
of  experimental  data  which  is  presented  in  Fig.  23. 

3.  The  sound  pressure  levels  of  the  discrete  components  are 
determined  from  formulas  (16)  and  (18)  and  will  pertain  to  the 
spectrum  of  broadband  noise  in  the  appropriate  frequency  bands. 

4 .  The  level  of  received  noise  is  determined  from  the 
compressor  noise  spectrum  in  1/3-octave  bands ,  or  from  the  semi- 
rational  formula  (22). 

Conclusions 

1.  The  basic  noise  sources  which  determine  the  acoustic 
characteristics  of  turbofan  engines  are  the  jet  and  the  compressor 
(fan).  Depending  on  the  bypass  ratio,  operating  conditions 
parameters,  dimensionality,  and  engine  power  rating,  either  the 
noise  of  the  Jet  or  the  compressor  turns  out  to  be  prevalent. 

Thus,  in  takeoff  mode  with  y  <  2  the  Jet  noise  predominates,  with 
y  >  2  -  the  compressor  noise.  In  landing  mode,  the  compressor 
noise  'ceeds  the  jet  noise  already  at  y  >  1. 

2.  The  acoustic  responses  of  the  Jet  (in  the  range  of  the 

gas  expansion  ratios  1.5  <  tt  <  2.2)  have  been  investigated  at 

pc 

present  with  sufficient  completeness.  A  strict  theoretical  base 
has  been  created,  on  the  basis  of  which  engineering  methods  of 
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calculation  have  been  developed.  Results  of  vast  experimental 
research  will  confirm  the  reliability  of  the  calculation  methods. 

3.  Compressor  (fan)  noise  has  been  little  studied,  these 
experimental  studies  are  insufficiently  clear,  and  the  obtained 
results  do  not  always  agree  with  each  other.  A  strict  acoustic 
theory  of  the  compressor  thus  far  does  not  exist,  therefore  during 
the  engineering  calculations  semi-empirical  dependences  are 
utilized.  The  primary  attention  of  the  researcher  during  research 
is  now  concentrated  on  the  physical  processes  of  noise  generation 
and  on  their  analytical  description,  which  itself  agrees  with 
experiment . 

Bibliography 


1.  B  .1  a  c  o  ■  E.  B  ,  Kbiitka  BF,  Mc.ibhiikob  6.  H.,  My- 
ii  n  ii  A.  T.  Oipaiiieiciiiie  my  mb  ca.wu.ioT(i«  (npucxr  coBercKoro  CTaimapu). 
«l'p;i/K.iaiuKiin  juii.imia»  AV  2.  1970. 

2.  I’iic.tfTU  n  iiJMipciiiin  .\apKKTepiicniK  iuy.ua,  cn3.iaBaeMoro  ■  .id.ii.hcm 
JB\  KOBIJM  r.O.IC  pP.lxrilBIIUMII  l  iMO.ICi  aMII.  Uo.l  pe.IBKIIIIcn  .1.  II.  CopKIlHB. 
M.,  Mai.uiiiocTpoi-iiiie,  1908. 

3.  RibnerH.  S.  Jcti  and  Noisc.CAS  Journal  vol.  14,  1969,  .Si  It). 

4.  Marin  A.  II.,  Me  Pike  A.  L.  Noise  Level  ol  Turbojet  —  and  Turbo- 
Ian  —  Powered  Aircraft.  Sound.  1963,  2,  >4  5. 

5.  C  vi  n  r.  X  a  \  i  Buy  ipeiiiuu'  ii.  ro'iiuiMi  myna  b  raioTy  pfiiiwiux  .mma- 
1C.IHX.  I |3MCPCIIIIH  II  leopilM.  TpSlIJ  A.MCplIK UlICKWO  OfllllCl  TBd  IIIDKHIcpOB-MO- 
xainlKOB,  Cepun  A  «3iuprelii‘uvKi!c  uaiiiiiiiu  ii  ycTanoiiKti*.  Pycciciifl  nepeno.i. 
At.,  «Mnp»,  .V*  2,  1967. 

6.  HayJ.  A.  The  Eitimatiuii  of  Jet  Engine  Noise.  Presented  ol  Aeron. 
Acoust.  Symp.,  Toulouse,  6-r-8,  Mar.,  1968. 

7.  Progress  ol  NASA  Research  Relating  to  Noise  Alleviation  ol  Large  Sub¬ 
sonic  Jet  Aircraft.  (NASA  SP  —  1-M»l.  Prrseutid  'I  Conference  Held  at  Langley 
Research  Center.  1968. 

h.  .'I  a  ft  t  x  ii.i.i  M.  lUy  .M  CTpytt.  PaKOTHaa  tc.miiikb  h  kocMoiiivBTinta. 
.\?  7.  1963. 

9.  Pi  a  nk  o  M. Fludr  thcorique  et  expcrimeiitalc  du  bruit  de  compiesseur 
Serv  Techn.  Aeron.,  Paris,  avr.  1969. 

10.  Kobryiiski  M.  Recentei  mesures  de  luuil  de  jet*  el  de  compres* 
sours.  ON'ERA,  T.  p.  Ni  6ti2,  l%8. 

11.  Prediction  ol  Tutbine  Engine  Compressor  or  Fan  Noise.  SAI!  Proposed 
air  973. 

12.  Bragg  S.  L.  B r i d  g  «•  R.  Noise  Horn  lurbojet  Compressors.  JRAS, 
1964,  68.  Nt  637. 

13.  <t>  p  a  ii  u  e  s  11.  K.  Oupc.ie.iemie  ypomui  uiy  mb  ii  iiapaMcipoH  uiyMor.iy- 
uiine.iH  ,.iiii..iuioiiiiu\  i  j.ioTyp6iiiiiiLi\  .loiiraTe.ivfi.  Pura,  1964. 


127 


SOME  PROBLEMS  OF  THE  CREATION  OF  AN  OPEN 
STAND  FOR  ACOUSTIC  STUDIES  OF  DTRD 

A.  I.  Balmakov,  V.  G.  Yenenkov 


In  the  article  the  contemporary  methods  of 
the  experimental  research  on  the  acoustic 
characteristics  of  DTRD  are  examined.  Primary 
attention  is  given  to  the  full-scale  experiment 
on  an  open  stand.  The  purposes  of  the  acoustic 
studies  of  DTRD,  and  the  basic  requirements  for 
test  conditions,  facilities  and  equipment  of 
the  open  stand  are  set  forth.  Two  versions  of 
the  structural  solution  of  the  stand  are 
described:  stationary  and  mobile.  Basic 
information  about  necessary  facilities  and 
equipment  of  the  open  stand  is  given. 


Introduction 


Recently  in  civil  aviation,  jet  aircraft  with  the  turbofan 
engines  with  high  bypass  ratios  y  have  been  accepted.  With  the 
use  of  such  DTRD  the  Jet  noise  of  the  exhaust  gases  decreased, 
but  the  noise  generated  by  the  fan  simultaneously  increased. 

To  solve  the  problems  of  turbofan  engine  noise  reduction, 
it  is  first  necessary  to  conduct  detailed  experimental  studies, 
both  model  and  full-scale. 


128 


•*!vv*w«*vri 


F 

i 

In  the  course  of  the  model  experiment  It  Is  also  possible 
to  t^st  a  large  quantity  of  diverse  variants  of  the  Investigated 
object  with  insignificant  expenditures,  to  study  the  mechanism 
of  the  noise  generation  by  different  sources  in  detail,  and  to 
develop  methods  of  its  reduction. 

Model  tests  are  carried  out,  as  a  rule,  in  muffled  (anechoic) 
chambers  (Pig.  1).  The  internal  surface  of  the  chambers  is 
covered  with  the  sound-absorbing  wedges,  which  makes  it  possible 
to  produce  acoustic  measurements  at  frequencies  above  250  Hz 
without  distortions  caused  by  reflections  from  the  walls. 


Fig.  1.  The  preparation  of  the  model 
experiment  in  an  anechoic  chamber. 


The  studies  of  noise  problems  of  individual  and  coaxial  Jets, 
in  connection  with  the  diagrams  of  the  nozzles  of  the  first  and 
second  contours  of  DTRD  the  power  Jets  of  supersonic  transport 
planes,  and  also  DTRD  compressors  (fans)  (Fig.  2),  are  carried 
out  in  anechoic  sound  chambers. 

Full-scale  experimentation  on  an  open  stand  is  necessary  for 
the  purpose  of  final  proof  of  these  studies  of  models,  and  for  the 
generalization  of  them  in  connection  with  turbofan  engines. 


-Vm,  i  >  XAftM&UI 


Fig.  2.  The  interior  of  an  anechoic 
chamber  with  a  model  of  a  muffler. 


Fig.  3.  An  overall  view  of  an 
engine  installation  on  the 
open  stand. 

The  open  stand  for  full-scale  DTRD  (Fig.  3)  makes  it  possible 
to  carry  out  studies  not  only  of  acoustic  characteristics,  but 
also  to  test  the  operational  reliability  of  the  engine. 
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In  the  open  stand  it  is  possible  to  make  a  study  of  the 
intake  and  exhaust  devices  (air  intakes,  engine  nacelles,  nozzles), 
a  highly  reliable  final  adjustment  of  the  elements  and  systems  of 
the  engine  by  applying  methods  of  the  early  flaw  detection,  etc. 

A  study  program  can  also  include  testing  of  thrust  reverse 
and  research  on  the  effect  of  the  rate  of  the  growth  of  load  and 
its  cyclic  recurrence  on  the  engine  lifetime.  The  open  stand  makes 
it  possible  to  study  the  questions  of  the  icing  of  air  intakes  and 
engine  nacelles  and  prevention  of  the  incidence  foreign  objects 
into  the  engine,  and  other  problems. 

From  the  aforesaid,  it  is  evident  that  the  multipurpose  open 
stand  which  has  concrete  purpose  (investigation  of  the  acoustic 
characteristics  of  DTRD)  makes  it  possible  to  study  a  wider  circle 
of  problems  of  theoretical  and  operational  nature.  Thus,  at  present 
the  creation  of  open  stands  is  extremely  necessary  and  real. 
Requirements  for  such  a  stand  very  high,  and  the  stand  itself 
should  be  universal. 

In  this  article  the  problems  of  the  creation  of  the  multipurpose 
open  stand  for  the  study  of  the  acoustic  characteristics  of  turbofan 
engine  are  brought  to  light. 

§  1.  Purposes  of  Experimental  Acoustic 
Studies  on  an  Open  Stand 


1.  The  Study  of  the  Characteristics  of  the 
Ambient  Engine  Noise  in  a  Distant  Acoustic 
Field 

On  the  open  stand  it  is  possible  to  produce  a  study  of  the 
characteristics  of  the  ambient  noise,  both  of  production  and 
experimental  engines,  in  distant  and  near  acoustic  field. 
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In  a  distant  acoustic  field  the  sound-pressure  decreases 
inversely  proportional  to  distance  from  the  center  of  radiation  to 
measuring  point.  During  stand  tests  it  is  desirable  to  produce 
measurements  at  a  distance  of  50-300  m. 

These  problems  enter  into  the  study  of  the  characteristics  of 
the  ambient  engine  noise  in  a  distant  acoustic  field: 

a)  the  determination  of  directional  characteristics  in  the 
levels  of  the  received  noise  and  in  sound  pressure  levels  during 
engine  operation  under  basic  conditions; 

b)  obtaining  spectral  characteristics  in  octave  (or  1/3  octave) 
frequency  bands  at  all  measuring  points; 

c)  establishment  of  the  maximum  level  of  the  received  noise 
created  by  the  engine; 

d)  the  plotting  of  curves  of  noise  reduction  In  the  direction 
of  the  maximum  propagation  at  a  distance  of  50-300  m  and  character¬ 
istic  noise  spectra  in  octave  frequency  bands. 

On  the  basis  of  the  findings  it  is  possible  to  give  a 
comparative  evaluation  of  engines  with  regard  to  ambient  noise. 

2.  The  Study  of  the  Characteristics 
of  Ambient  Noise  in  a  Near  Acoustic 
Field 

In  the  immediate  proximity  of  the  engine  there  is  a  field  of 
hydrodynamic  pressure  which  is  called  the  near  acoustic  field. 

During  the  study  of  the  characteristics  of  DTRD  noise  the  need 
for  producing  measurement  in  the  near  field  frequently  arises. 
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The  following  enters  into  the  study  of  the  characteristics  of 
near  field  noise: 

a)  determination  of  sound  pressure  levels  under  basic 
conditions  of  engine  operations; 

b)  determination  of  directional  characteristics; 

c)  obtaining  spectral  characteristics  in  octave  bands  at  the 
measuring  points  of  the  near  field; 


j 
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d)  establishment  of  the  maximum  noise  level. 


The  numerous  experimental  studies  of  the  near  field  show,  that 
the  greater  the  signal  frequency,  the  less  the  distance  of  the 
sound  pressure  maximum  from  the  source. 

j 

1 

3.  The  Study  of  Methods  of  Noise  Reduction 

3 

The  basic  purpose  of  the  open  stand  for  acoustic  studies  of 
DTRD  should  be  considered  the  development,  testing,  and  evaluation 
of  the  effectiveness  of  the  different  methods  of  noise  reduction. 


The  basic  principle  of  aviation  noise  reduction  is  the 
principle  of  the  suppression  of  its  sources.  For  DTRD,  this 
principle  is  realized,  for  example,  by  following  solutions: 


a)  an  increase  in  the  bypass  ratio  for  the  purpose  of  a 
reduction  in  the  Jet  noise; 


b)  removal  of  the  inlet  guide  ring; 


c)  an  increase  in  the  clearance  between  vane  rings,  etc. 
for  the  purpose  of  reduction  in  compressor  noise. 
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This  principle  is  most  labor-consuming,  requires  prolonged 
experimental  research  and  change  of  engine  design.  But  this 
principle  is  also  the  most  effective. 

Another  principle  obtained  wide  dissemination  -  the  principle 
of  the  energy  absorption  of  noise  on  the  path  of  its  propagation. 
The  sound-absorbing  facing  of  the  channels  at  the  inlet  and  exit 
from  the  compressor  (fan)  works  according  to  this  principle.  A 
change  in  the  form  of  the  channels  of  t,he  inlet  and  exit  gives 
the  same  result.  For  the  facing  of  the  intake  and  outlet  ducts  of 
the  fan  contour,  a  porous  sound-absorbing  material  is  used.  For 
example,  for  the  RB211  engine  installed  in  the  Lockheed  L-1011 
airbus,  18.5  m  of  facing  material  (Fig.  4)  is  required  for  the 
indicated  purposes. 

Fig.  4.  Treatment  of  exhaust 
ducts  of  a  turbofan  engine 
with  the  sound-absorbing 
facing. 


§  2.  Requirements  for  Acoustic  Tests 

1.  Assurance  of  Conditions  of  a  Distant 
Acoustic  Field 

This  requirement  should  be  fulfilled  maximally  strictly,  since 
the  accuracy  of  measurement  of  all  acoustic  characteristics  depends 
on  it . 
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Near  the  noise  sources,  the  points  of  the  arrangement  of 
microphones,  and  between  them,  there  should  not  be  any  obstructions 
(Fig.  5).  Since  the  noisy  object  (DTRD)  has  a  determinate  length 
from  air  intake  to  a  certain,  not  completely  clearly  limited,  but 
distant  from  the  nozzle  edge,  cross  section,  during  the  study  of 
the  directional  characteristic  and  sound  pressure  levels  in  the 
distant  field,  it  is  necessary  to  make  noise  measurements  at  a 
sufficient  distance  from  the  engine. 


Fig.  5.  Overall  view  of  the  open  stand. 


Fig.  6.  The  effect  of 
engine  dimensionality 
on  the  value  of  the 
necessary  distance 
between  the  engine  and 
microphones,  with  which 
the  conditions  of  distant 
field  are  provided. 
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The  value  of  the  necessary  distance  from  the  studied  object 
to  the  points  of  the  microphones  arrangement  depends  on  the  engine 
model  and  its  dimensionality  (Pig.  6).  Then  the  observer  can 
consider  the  engine  as  a  point  source  of  noise. 

It  is  necessary  also  to  eliminate  the  effect  of  the  ground 
in  order  to  avoid  the  interaction  of  the  Jet  with  soil. 

Extraneous  noise  sources  should  be  at  a  sufficient  distance 
in  order  to  not  introduce  distortion  into  the  acoustic  measurements. 

2.  Requirements  for  Stand  Equipment 

The  test  bench  should  be  equipped  with  devices  for  measuring 
thrust,  fuel  consumption,  temperature,  and  gas  pressure  in  basic 
sections  of  the  DTRD.  The  requirements  for  the  accuracy  of 
measurement  of  the  parameters  are  Just  as  rigid  as  during  the  study 
of  the  characteristics  of  the  effectiveness  and  cost-effectiveness 
of  the  engine,  since  the  latter  are  directly  connected  with  the 
characteristics  of  noise. 

Test  floor  should  be  equipped  with  equipment  for  continuous 
recording  of  the  rate  and  wind  direction,  of  the  barometric 
pressure  and  air  humidity.  These  data  substantially  affect  the 
distribution  of  sound.  At  wind  velocities  7.5-8  km/h,  measurements 
in  a  distant  acoustic  field  are  not  stable. 

It  is  necessary  to  have  a  sufficient  quantity  of  microphones 
for  simultaneous  recording  of  data  measurements  in  many  points  in 
order  to  eliminate  errors  from  a  change  in  wind  velocity, 
refraction,  conditions  of  engine  operation,  etc.  The  use  of  a 
large  number  of  microphones  will  make  it  possible  to  decrease  the 
combustion  duration  and  to  lower  consumption  for  experiment. 
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During  the  study  of  intake  and  outlet  devices  of  unsymmetric 
form,  it  is  necessary  to  produce  measurements  in  both  the 
horizontal  and  vertical  planes  passing  through  the  engine. 

3.  Requirements  for  Audio  Equipment 

When  selecting  metering  equipment  for  acoustic  studies  it  is 
advantageous  to  use  the  recommendations  of  International 
Electrotechnical  Commission  ( 1 9 6 1 ) . 

The  frequency  band  of  the  measuring  circuit  should  be  within 
the  limits  of  20  to  12,500  Hz.  Measuring  circuit  in  this  case 
should  conform  to  the  following  requirements: 

1)  nonuniformity  of  the  frequency  characteristic  of  measuring 
circuit  in  every  octave  band  should  be  not  more  than  ±3  dB; 

2)  error  for  measuring  circuit,  taking  calibration  into 
account,  should  be  not  more  than  ±2  dB; 

3)  level  of  internally-produced  noise  should  be  5  dB  lower 
than  maximum  level  of  measuring  noise; 

4)  microphone  should  be  weakly  directed  in  the  entire  range  of 
operating  frequencies; 

5)  measuring  circuit  should  be  calibrated,  and  before  beginning 
and  after  measurements  -  recalibrated. 

The  obtained  acoustic  characteristics  should  be  supplemented 
by  the  following  information  fixed  into  the  protocol  of  the  tests : 

-  date; 

-  engine  model,  its  series  and  operating  mode; 

-  atmospheric  conditions  near  the  earth  (relative  air  humidity, 
temperature,  barometric  pressure,  the  velocity  and  wind 
direction) ; 


-  the  layout  of  the  measuring  points  (radius  measurement, 
number  of  points,  etc.). 

§  3-  The  Types  of  Open  Stands  for  Full-Scale 
Acoustic  Investigations 

A.  The  Stationary  Open  Stand 

For  the  creation  of  the  stationary  open  stand  the  presence  of 
free  areas  of  sufficient  large  dimensions  is  necessary.  The  design 
of  such  a  stand  can  be  analogous  to  the  design  of  a  stand  for  usual 
engine  tests  outside  the  cell.  The  stand  is  the  experimental 
station  which  has  the  engine  and  its  maintenance  systems.  The  tool 
for  the  measurement  of  thrust  is  fastened  to  a  general  mounting 
which  is  immobilely  connected  tc  the  load-bearing  elements  of  the 
foundation.  The  stand  is  established  on  an  open  site  distant  from 
projecting  buildings  and  constructions  and  also  from  noise  sources. 

In  order  to  eliminate  the  effect  of  the  ground,  it  is  necessary 
to  raise  the  engine  by  a  height  of  1. 7-2.0  m  and  higher,  depending 
on  the  overall  dimensions  cf  the  engine  and  its  arrangement  on  the 
aircraft.  It  Is  most  advisable  to  completely  or  partially  bury 
the  control  room  and  the  fuel  tank  in  the  ground,  10-15  m  from  the 
engine.  In  this  case,  all  power-supply  and  engine  control  systems 
will  pass  in  bunkers  and  have  comparatively  low  length.  It  is 
convenient  to  locate  the  entire  control  panel  to  the  side  of  the 
engine,  armoured  with  slits  for  the  survey  of  the  engine  during 
operation  (Fig.  7). 

The  application  of  a  ground-level,  mobile,  control  panel  is 
possible.  It  must  be  removed  to  a  considerable  distance  (20-30  m) , 
which  makes  survey  worse,  complicates  throttle  circuits,  and 
maintenance . 
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^  13.  Foundation 


i  Basic  entries 
•"v  5.  Armour  plate 
j 6.  The  control  panel 

7.  Control  room  (shack) 

8.  Emergency  exit 

9.  Fuel  tank 

10.  Bunkers  for  the  power-supply  and  control 
systems 


Fig.  7.  The  layout  of  engine  and  systems  on  the  open  stand. 


During  the  rational  selection  of  all  parameters  a  stationary 
stand  possesses  the  advantages,  the  basic  of  which  are:  simplicity 
of  the  design,  maintenance  and  organization  of  the  measurements  of 
noise  at  different  points  (site  can  be  equipped  with  stationary 
all-weather  microphones). 
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B.  The  Mobile  Open  Stand 


Frequently  the  scientific  research  organizations  are  not 
furnished  with  large  free  areas  for  a  stationary  open  stand.  Then 
it  is  necessary  to  create  mobile  open  stand  (Fig.  8).  It  can  be 
mounted  on  the  base  of  a  chassis  (trailers,  etc.).  The  condition 
of  free  field  is  satisfied  in  this  case  sufficiently  simply,  since 
there  is  the  possibility  of  producing  the  selection  of  a  testing 
place  (field,  wood  clearing,  etc.)  by  the  transportation  of 
movable  stand. 

For  the  delivery  of  the  stand  into  the  area  of  tests,  the  tow 
car  or  a  special  airport  mobile  power  unit  [AMPU]  ( AIHA )  is  utilized. 
In  the  latter  case  the  problem  of  starting  the  engine  is  considerably 
simplified. 

It  is  desirable  to  transport  the  mobile  open  stand,  and  the 
equipment  available  in  the  airport  weather  service  makes  it 
possible  to  obtain  the  necessary  atmospheric  characteristics.  In 
this  case,  the  maintenance  of  the  engine  and  its  systems,  fuel 
loading,  fuel  and  lubricants,  and  special  fluids  are  simpler. 

For  test  work,  the  mobile  stand  is  established  on  a  level 
site,  the  supports  which  impart  lateral  stability  to  the  stand  are 
extended,  and  the  chassis  wheels  are  chocked.  All  systems  of  the 
fueling  and  control  of  the  engine  are  furnished  on  an  isolated 
frame,  which  can  be  moved  in  height  with  the  aid  of  rigid  adapters, 
or  on  self-raising  frame  with  a  hydraulic  cylinder. 

For  a  change  of  the  slope  angle  of  the  engine,  in  the  vertical 
plane,  the  truss-adapters  between  frames  should  have  different 
heights.  When  using  hydraulic  hoists  in  such  cases,  the  rigid, 
mechanical  locking  of  its  legs  is  necessary. 


Fig.  8.  Engine  installation 
of  a  mobile  platform. 

During  the  arrangement  of  new  position,  the  upper  frame  with 
engine,  panel,  etc.,  is  raised  by  a  crane  or  hoists,  fixed  in 
a  determined  position,  and  then  both  frames  are  rigidly  fastened. 
The  use  of  other  methods  of  a  change  in  the  height  of  the  engine 
installation  is  possible  also. 

Also,  the  control  panel,  fuel  tank,  the  system  of  thrust 
measurement,  and  part  of  the  measuring  apparatus  are  installed 
on  the  mobile  frame. 

The  control  panel  in  this  case  can  be  utilized  for  the 
limitation  of  conditions  regarding  noise  in  the  cabin  of  the 
aircraft  (during  the  use  of  similar  sound-absorbing  material, 
engine  power  rating,  etc.). 

At  the  location  of  conducting  the  research,  it  is  necessary 
to  make  the  arrangement  of  microphones  and  to  lay  cable  to  them. 
It  must  be  considered  that  stand  variations  will  occur,  and 
therefore  it  is  necessary  to  provide  the  protective  measures  for 
especially  precise  and  sensitive  equipment  on  the  control  panel. 
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The  mobile  open  stand  is  considerably  more  complex  in  design 
than  the  stationary,  but  is  also  more  universal. 

§  4.  The  Basic  Equipment  of  the  Open  Stand 

1.  For  equipping  the  stand  with  a  thrust  measurement  system, 
a  non-rigid  test  stand  design  on  connecting  rods  or  flexible  strips 
may  be  used  (Fig.  9). 
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Fig.  9.  A  closeup  of  the  engine, 
prepared  for  noise  measurements 
on  a  open  stand. 

2.  In  turbofan  engine  the  temperature  of  gases  and  pressure 
in  different  cross  sections  should  be  measured. 

3.  For  the  measurement  of  pressures  (complete  and  static) 
on  a  test  stand,  different  adapters,  manometers,  vacuum  gages 
and  differential  nanometers  can  be  applied. 

4.  For  flow  velocity  measurements  (in  Mach  numbers  of  X), 
they  use  a  combined  adapter  of  complete  and  static  pressures,  and 
also  measure  complete  or  static  temperature.  Directivity  of  the 
velocity  vector  can  be  determined  with  the  aid  of  multichannel 
adapters . 
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5.  For  the  measurements  of  fuel  consumption,  tachometer 
sensors,  flowmeter,  and  also  volumetric  and  mass  flowmeters  can  be 
applied.  Flow  meters  are  usually  used  for  more  precision 
determination  of  the  consumption  with  regard  to  engine  operating 
mode.  For  the  measurement  of  instantaneous  fuel  consumption  they 
use  tachometer  sensors. 

6.  For  the  measurement  of  r/min  electrical,  strobotacs  and 
the  total  revolution  counters  are  applied.  All  measured  values 
enumerated  above  are  available  on  control  instruments,  the  panel 
(in  the  cabin)  according  to  which  operation  of  the  engine  and  its 
systems  is  evaluated:  the  engine-control  lever,  stopcock,  the 
starting  switches,  the  thrust  indicators,  tachometer,  the  thermo¬ 
meter  of  exhaust  gases,  the  oil  thermometer,  the  fuel  and  oil 
pressure  gauges,  fuel-metering  device,  the  signaling  of  the 
position  of  the  organs  of  nozzle  control,  compressor,  etc. 

7.  The  test  stand  is  supplied  with  the  following  of  systems: 

a)  by  the  electric  power  wiring; 

b)  6  or  12  V  electric  wiring  for  operation  of  the  portable 
illumination  sources  and  power  tools; 

c)  27  V  direct  current  electric  wiring,  that  provides  operation 
of  meters  and  starting  the  engine  and  turbine  starters  from 
electric  starters; 

d)  wiring,  carrying  the  load  of  the  engine's  generators; 

e)  compressed  air  feed  system; 

f)  oil  system; 

g)  hydraulic  circuit  (in  the  presence  of  hydroaggregates); 

h)  propellant  feed  system. 
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S  5.  The  Measurement  of  Engine  Acoustic 
Characteristics 

The  study  of  engine  acoustic  characteristics  has  the  purpose 
of  measuring  and  registering  sound  pressure  levels  at  a  certain 
radius  from  the  engine.  In  connection  with  the  fact  that  the 
domestic  industry  does  not  produce  standardized  audio  equipment 
which  makes  it  possible  to  produce  complete  acoustic  measurements 
under  varied  conditions,  the  measurements  of  sound  pressures  can 
be  made  with  the  aid  of  the  following  assemblies  of  equipment: 

a)  for  the  registration  of  levels  to  135  dP  use  of  the 
equipment  manufactured  by  Danish  firm  Briiel  and  Kjaer  A/S  is  more 
expedient,  including  the  chart  recorder  (Type  2304),  analyzer 
(Type  2110)  and  capacitor  microphones  (Type  4111,  4112,  4113). 

It  is  possible  to  also  apply  domestic  capacitor  microphones  of  the 
MK-5A  and  MK-6  type,  as  well  as  electrodynamic  microphones  of  the 
type:  MD-45,  MD-59,  MD-62,  etc.; 

b)  for  registration  of  levels  greater  than  140  dB  it-  is 
possible  to  use  the  domestic  MIK-5  microphone  with  the  UPU-2 
amplifier,  the  vacuum  tube  MIK-5  voltmeter  with  the  UPU-2 
amplifier,  or  the  vacuum  tube  voltmeter  and  analyzer  and  the  chart 
recorder  of  the  firm  of  Briiel  and  Kjaer  A/S; 

c)  for  the  magnetic  recording  of  noise  on  the  ground  and  in 
flight,  a  modified  magnetic  tape  recorder  (for  example  M-30,  MEZ-28S 
MEZ-41,  MEZ-63,  Kometa,  Yauza-10,  Reporter-3,  etc.).  For  a  magnetic 
recording,  tape  type  6  is  recommended. 

The  magnetic  recording  of  noise  is  reproduced  on  tracing 
paper  in  the  form  of  the  total  noise  level  and  frequency  components 
(with  the  aid  of  the  Briiel  and  Kjaer  analyzer  and  chart  recorder). 
The  modification  of  the  magnetic  tape  recorder  consists  of  the 
conforming  of  its  input  and  output  with  the  equipment  connected 
to  it ; 
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d)  for  a  selective  control,  in  terms  of  total  noise  levels, 
in  the  experiment  there  should  be  a  reliable  and  thoroughly 
calibrated  objective  audio-noise  meter  (of  Type  Sh-60,  Sh-63  and 
Sh-52,  Sh-2,  General  Radio’s  759A  or  similar),  and  also  domestic 
octave  filters  OP-1,  1/2-octave  PF-1,  1/3-octave  AW-2M,  etc. 

For  frequency  response  analysis,  the  octave  filters  of  the 
firm  of  Cruel  and  KJaer  of  Type  1612,  1613,  1/3-octave  -  2111, 

2112,  domestic  narrow-band  analyzers  S 5—3  and  S4-7  and  an 
analyzer  of  Type  2107  from  the  firm  of  Bruel  and  KJaer  are  applied. 

For  the  registration  of  level  -  monitors  N-110,  Type  2305  from 
Bruel  and  KJaer;  for  calibration  a  pistonphone  of  Type  4220,  a 
ball  bearing  calibrator  of  Type  4240  and  an  instrument  of  Type  4142 
from  Bruel  and  KJaer.  Also  applied  are  the  audiofrequency 
oscillators.  Type  ZG-4,  ZG-12,  GAII-1  and  GNII-2,  voltmeters  V-3-13, 
VK2-6,  oscillograph  S-l-19,  etc. 

During  the  investigation  of  acoustic  characteristics  it  is 
necessary  to  maintain  certain  conditions. 

1.  Before  beginning,  and  during  the  period  of  measurements, 
the  entire  equipment  should  routinely  undergo  calibration.  These 
adjustments  are  dated  and  are  introduced  into  these  equipment. 
Equipment  from  the  firm  Bruel  and  KJaer  may  be  calibrated  only 
once  -  before  beginning  the  experiments  (if  in  the  process  of 
measurements  sharp  Jolts  and  impacts  are  not  undergone). 

Besides  standardization,  and  independently  of  it,  the 
equipment  before  every  measurement  and  after  it  should  be 
calibrated. 

2.  When  using  Bruel  and  KJaer  equipment  in  field  conditions, 
one  ought  to  transport  to  the  place  of  experiment  in  a  warm 
vehicle  and  on  a  spring  suspension.  During  measurements  the 
equipment  remains  in  the  vehicle,  and  only  lengtheded  cables  to 
the  microphones  are  Installed.  The  microphone  cannot  be  shielded. 
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3.  If  the  tape  recording  of  noise  is  not  immediately 
reproduced,  then  the  magnetic  tape  is  placed  in  iron  cover  (for 
protection  from  stray  fields).  Between  finish  of  the  recording 
and  at  the  start  of  its  reproduction  there  is  a  period  no  more 
than  days. 

4.  At  wind  velocity  more  than  7  m/s  one  should  not  make 
measurement  of  ambient  noise. 

5.  A  change  in  the  temperature  of  air  with  height  leads. 

Just  as  the  presence  of  wind,  to  refraction.  If  temperature 
decreases  with  height,  then  sonic  rays  are  bent  upwards  and  the 
distance  of  signal  reception  decreases  and  vice  versa.  With  a 
drop  in  the  temperature  of  more  than  3.5-4°C  at  heights  from  0 
to  500  m  one  should  not  make  the  measurement  of  noise. 

6.  One  must  take  into  account  that  audio  equipment  (Including 
the  firm  of  Briiel  and  KJaer)  has  only  one  channel  of  the  recording 
of  noise;  it  records  either  the  total  noise  level  or  the  level  in 
any  determined  frequency  band  (octave,  1/3-octave,  1/4-octave). 

In  the  absence  of  a  large  number  of  microphones,  measurement 
can  be  produced  as  shown  in  Fig.  10  with  the  aid  of  one  microphone 
which  is  moved  on  a  circular  arc,  in  the  center  of  which  the  noise 
source  is  located  -  the  engine.  With  a  change  in  radius,  the 
process  of  measurement  is  repeated. 


Fig.  10.  Overall  view  of  the  cell  for 
acoustic  tests. 
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Conclusions 


1.  At  the  contemporary  stage  of  studies  of  the  acoustic 
characteristics  of  turbofan  engine  with  high  bypass  ratios,  where 
the  determinate  noise  source  is  the  compressor  (fan),  experiments 
on  model  and  full-scale  objects  have  great  significance.  Model 
experiment  is  carried  out,  as  a  rule,  in  anechoic  chambers.  To 
conduct  a  full-scale  experiment,  the  creation  of  an  open  stand  is 
necessary . 

2.  The  open  stand  for  the  acoustic  tests  of  DTRD  can  be 
executed  in  stationary  and  mobile  versions.  It  should  be  equipped 
with  all  systems  necessary  for  the  measurement  of  the  parameters 
of  the  effectiveness  and  cost-effectiveness  of  the  engine  and  its 
acoustic  characteristics. 

3.  The  open  stand  which  has  as  a  basic  purpose  acoustic 
investigations,  may  obtain  multipurpose  application  -  for  tests 
of  the  characteristic  of  engine  operational  reliability. 

Clearly  formulated  requirements  for  an  open  stand  for 
acoustic  tests  of  DTRD  thus  far  do  not  exist,  and  available 
separate  recommendations  do  not  reflect  a  complete  picture  of  the 
specific  character  of  the  studies.  Thus  it  is  necessary  to  develop 
commercial  norms  and  standards  which  determine  the  conditions  of 
assuring  distant  acoustic  field,  stand  equipment  and  audio 
equipment . 
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